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The infrared spectra of HeS, HDS and D,S have been examined up to 15 uw with a mono- 
chromating prism spectrometer, and the results compared with the Raman spectra. The 
fundamental frequencies have been assigned with reasonable certainty, and a calculation has 
been made of the frequencies corresponding to infinitesimal amplitude. The cruder methods of 
the vibrational-rotational spectrum suggest a vertical angle of 110°, as compared with 92° 


obtained previously from a rotational analysis. 





(ore have been many investigations of the 

infrared absorption spectrum of hydrogen 
sulphide.'~7 The possible presence of impurities 
has led to the postulation of bands which do not 
exist,':?»4 and to the denial of others which are 
really there.’ > We have re-examined the spectra 
of this compound and its deuterium isotopes be- 
tween 1 and 15 w with a monochromating prism 
spectrometer. 

The weakness of absorption throughout the 
spectrum has increased the difficulties of observa- 
tion and interpretation; and although finality is 
not claimed, we think we have presented at least 
a concordant account of the three isotopic mole- 
cules, while we have stressed in what follows the 
uncertainties which would justify further in- 
vestigation. 





'W. W. Coblentz, Pub. Carnegie Inst., No. 35, pp. 52 
and 178 (1902). 

* A. H. Rollefson, Phys. Rev. 34, 604 (1929). 

wa H. Nielsen and E. F. Barker, Phys. Rev. 37, 727 
(A¥YS1), 

*W. Mischke, Zeits. f. Physik 67, 106 (1931). 

ms H. H. Nielsen and A. D. Sprague, Phys. Rev. 37, 1183 
(1931), 

aie D. Sprague and H. H. Nielsen, Phys. Rev. 43, 375 
(1933), 

*P. C. Cross, Phys. Rev. 47, 7 (1935). 
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EXPERIMENTAL 
1. Apparatus 


The general arrangement has been fully de- 
scribed,® the only important modifications being 
the replacement of the Downing galvanometer by 
a photo-relay using a Zernicke and a Moll gal- 
vanometer, and the fitting of the absorption cells 
and thermopile with potassium bromide end- 
plates. We would emphasize the need for great 
care in the calibration of a prism spectrometer, 
particularly in the use of the unresolved maxima 
of atmospheric CO: bands, since the position of 
these changes considerably with the percentage of 
carbon dioxide present. In our experience the use 
of two absorption cells, one evacuated and the 
other containing the gas under examination, 
seems to be essential. The omission of this pre- 
caution inevitably leads to the inclusion of false 
absorption regions, and to a good deal of imagin- 
ary resolution. 

A preliminary run was made with the rocksalt 
prism, and the located regions of absorption were 
then more closely examined with the appropriate 


' A. B. D. Cassie and C. R. Bailey, Proc. Roy. Soc. 
A132, 252 (1931). 
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Fic. 1. The infrared spectra of HS, DoS and HDS. 


prism of quartz, fluorite, or rocksalt. Readings 
were taken at intervals of approximately 0.02 u 
in order to have full advantage of the available 
resolution. The slit widths used and the spectral 
interval included are as already given.® Correc- 
tions were made for the temperature of setting 
and of observing, and for the error in setting at 
minimum deviation. The absorption cell length 
was 45 cm. 


2. Preparation of the gases 


H.S was prepared by three different methods: 
(i) by warming calcium sulphide with a saturated 
solution of magnesium chloride; the gas was then 
bubbled through water and baryta solution and 
over P,O;: it was then liquefied by solid CO2 and 
fractionated; (ii) by passing H2S into a suspen- 
sion of magnesium oxide in water, and gently 
warming the resulting Mg(HS)2: the gas so ob- 
tained was washed and dried as before; (iii) as a 
check, a sample was prepared from aluminum 
sulphide: an apparatus was devised which could 
be swept out with dry, CO2-free air and evacu- 
ated : water was dropped onto Al,S; from within 
the apparatus and the resulting HeS passed over 
baked calcium chloride and P:O;, and collected 
in a trap from which it could be fractionated ; the 
sample so prepared was found to be free from 
COs. It was found that neither KHS nor baryta 
saturated with H2S removed CO: completely. 

D.S was first obtained by method (iii) above, 
with 99.2 percent D.O. Additional Al.S; was 
introduced as a drying agent. Curve (a) (Fig. 2) 
shows that the sample contained HDS, probably 


°C. R. Bailey, A. B. D. Cassie and W. R. Angus, Proc. 
Roy. Soc. A130, 134 (1930). 
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derived from the calcium chloride and Al.S;. The 
former was therefore discarded, and the latte 
baked in vacuum before use, when curve (b) (Fig 
2) was obtained : it will be seen that although the 
amount of HDS was considerably reduced, com 
plete elimination was not effected. 

HDS was obtained by mixing equal volumes 
of HeS and D.S, and allowing the mixture to 
come to equilibrium. 


EXPERIMENTAL RESULTS 


A general comparison of the three spectra is 
given in Fig. 1, from which a rough measure of 
the relative intensities can be made. The adjusted 
percentage absorptions of Figs. 2-4 should not be 
used for that purpose. 


H.S 


Band A: 8.089 w; 1236 cm (rocksalt); Fig. 2 
There has been much discussion as to the exist- 
ence of this band, but we have found it in all 
samples of H2S, and all doubt is removed by the 
appearance of the corresponding isotopic bands. 
Although a good deal of the rotational structure 
is revealed, there is no evidence of separation into 
P and R branches, and the band is altogether 
different in appearance from the corresponding 
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fundamentals of DO and H.O. It is just possible 
that the band may be overlaid by a difference 
tone. A number of the most prominent lines show 
a separation of the order of 7 cm". 

Band B: 4.127 wp; 2423 cm (fluorite). This 
band of low intensity is probably real and 
different from that located at 4.25 u by Coblentz,! 
Rollefson,? and Mischke,* which is most likely 
due to COs. 

Bands C and D: 3.8 uw; (fluorite) Fig. 3. In this 
region there are peaks at 3.91 4 (2558 cm), 
3.70 uw (2703 cm~'), and 3.435 uw (2911 cm="). The 
separation between 3.70 and 3.43 u is 208 cm“, 
and is much too large to enable these peaks to be 
considered as a doublet: we therefore conclude 
that 3.70 u consists of the unresolved Q and R 
branches of a band, C, whose P branch is at 
3.91 uw. The line (—7) of Nielsen and Barker’s 
Curve C* may perhaps be the weak Q branch. 

Band E: 2.636 uw: 3794 cm (quartz); Fig. 4. 
The contour shows the presence of P, Q, and R 
branches in agreement with the complete resolu- 
tion obtained by Nielsen and Barker.’ 

Band F: 1.94 w; 5155 cm (quartz). The band 
was not resolved, but has been shown to be a 
doublet by Sprague and Nielsen,® whose value 
for the center is 5130 cm™. 
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Fic. 3. Bands C and D. 
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Fic. 4. Band E. 


D.S 


Band A: 11.14 w; 898 cm (rocksalt) ; Fig. 2 (b). 
The band is similar in structure to Band A of 
H.S; there is some faint evidence of rotational 
structure, but a spacing of 3.5 cm™! corresponding 
to the 7 cm~' is hardly within the limits of resolu- 
tion of our spectrometer. 

Band C: 5.155 uw; 1940 cm (fluorite); Fig. 3. 
The band was extremely weak, but has the one- 
sided appearance often associated with the pres- 
ence of a Q branch. 

Band E: 3.57 uw; 2798 cm (quartz) ; Fig. 4. The 
band resembles that of H2S, but again no resolu- 
tion was obtained. 


HDS 

Band A: 9.28 uw; 1078 cm™ (rocksalt); Fig. 2. 
The shape of the band serves to emphasize the 
difficulties in the interpretation of the spectra. 
Two peaks are observed at 1078 and 988 cm“. 
The separation between them is some 90 cm™', a 
not unlikely value for a P-R spacing in HDS; but 
the total frequency range included in the band is 
some 340 cm™! against 180 for HeS and 135 for 
D.S. We are faced with the absence or extreme 
weakness of a P branch in the latter substances, 
or alternatively the appearance of a new band in 
HDS. The determination below of the force 
constant ks makes it certain that the three peaks 
1236, 898, and 1078 cm“! correspond in the three 




























































628 BAILEY, 
substances. Since HDS is asymmetric, we for this 
molecule should possess a small component of 
change of electric moment along the least axis 
of inertia, and the band should therefore show a 
weak Q branch, and it is tempting to find in the 
sharper character of the peak at 1078 cm~™ some 
faint evidence for its existence. 

Band C at 5 u: (fluorite); Fig. 3. There is a 
much extended region of weak absorption from 
4.5 to 6.04, with peaks in the neighborhood of 
4.7 and 5.1 uw. 

Band E’: 3.5 uw; (quartz); Fig. 4. The peak at 
2721 cm is probably due to H2S, and that at 
2938 cm to DeS and HDS (see assignments). 

Band E"’: 2.7 uw; (quartz); Fig. 4. The peaks at 
3724 and 3849 cm™ are probably both due to 
HDS and to HS. 


THEORETICAL 


A very striking and complete analysis of the 
band at 9911 cm has been carried out by Cross? 
from measurements in the photographic infrared. 
The bond angle in the normal state was found to 
be 92° 20’, while for water vapor the correspond- 
ing value is 105°6’. It is difficult to predict 
qualitatively the effect upon the vertical angle of 
the exchange of sulphur for oxygen: the distance 
of the hydrogen atom from the central atom in- 
creases from 0.97 to 1.34A so that lessened inter- 
action between the extreme hydrogen atoms 
would favor closer approach to the theoretical 
valence angle for p-s binding of 90°; on the other 
hand the electronic structure of the sulphur atom 
tends to make the angular rigidity less. 

An investigation of a purely vibrational spec- 
trum will not yield results of the same high order 
as those obtained by Cross, but for the purpose of 
comparison with heavier molecules whose rota- 
tional structure cannot as yet be unravelled, it is 
interesting to enquire what accuracy may be 
expected when the number of approximations is 
reduced as far as possible.We therefore propose to 
apply the potential energy expression used by 
Van Vleck and Cross in their consideration of the 
water vapor molecule.'' For this purpose we re- 
quire an estimate of the mechanical frequencies 
v; associated with infinitesimal amplitudes; and 

10R. Mecke, Zeits. f. Physik 81, 313 (1933). 


1 J. H. Van Vleck and P. Cross, J. Chem. Phys. 1, 357 
(1933). 
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we now proceed to an attempted assignment of 
frequencies. 


ASSIGNMENTS 


We will adopt the conventional numbering of 
Dennison, in which v3 is the asymmetrical fre- 
quency, v; the short wave and v2 the long wave 
symmetrical frequencies. The limiting angle 
above which the axis of least inertia is perpen- 
dicular to the symmetry axis is 88° 16’ for HS, 
and for D2S is 86° 38’: hence we may expect 7; to 
have a more or less well developed Q branch. A, 
the smallest moment of inertia, is about 0.9B, 
whereas for water, A=0.5B; and from Denni- 
son’s curves” we see that this Q branch will be 
much weaker in H2S than in HO. This may 
account for our inability to resolve the band. For 
HS we therefore take 2632 cm™ as w3. The peak 
at 2911 cm™' gives a force constant of some 
4.9X 10° dynes per cm which seems to be too high 
to allow this band to count as an H2S fundamen- 
tal, although the possibility should not be over- 
looked. We have to confess that we can offer no 
reasonable assignment for it, and since it appears 
in all samples of the gas we do not think it can be 
ascribed to an impurity. w; is found in the Raman 
spectrum of the gas at 2615 cm™'® w. we must 
place at 8 uw, although the peculiar structure of 
the band again makes it difficult to determine the 
center exactly: a comparison of our own and 
Rollefson’s curve® suggests 1236 cm™ as the most 
reasonable value, and the data for the isotopic 
molecules are in agreement. 2432 cm is then 
2w2, and since 3794 cm™ has a Q branch it must 
be w3-+we. As far as one can say from our own 
work and that of Sprague and Nielsen® the band 
at 1.9 u has no Q branch and may be 2a, rather 
than w,;+w3. Weizel™ has shown for HO that 
even overtones of ws are not to be expected, and 
we may probably transfer his reasoning to this 
case. 9911 cm (w3+3w:) and 10,100 cm™ 
(3w3-+w1) have already been assigned by Cross.’ 

D.S. We assign 898 cm~! to we, 1940 cm“, to 
w3, and 2793 cm™ to w3+we2. The only reference 
we have been able to find to the Raman spectrum 
occurs in an authors’ abstract of a paper pre- 


sented last year to the Vienna Academy ;" as far 


2D. M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 
18S. Bhagavantam, Nature 126, 502 (1930). 

144 W. Weizel, Zeits. f. Physik 88, 214 (1934). 

1 A, Dadieu and H. Kopper, Nature 136, 235 (1935). 
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INFRARED ABSORPTION SPECTRA 


Proposed assignment of Raman and infrared 
bands of HS, HDS, and D2S. 


TABLE I. 








HDS 
(em~) 


1078 
(1910) calc. 
(2619) calc. 

2110 
3724 
2938 


H.S 
(em) 


1236 
2615 
2632 
2423 
3794 


D.S 


(cm~) 


898 
1906 
1940 





2798 


5155 
9911 
10100 











as we can say, the paper has not yet been pub- 
lished. The mixture of H2S and DS gives lines at 
1880 and 2585 cm™'. These are certainly the 
values for the liquid state, and we find for D,O 
and H,0 that the percentage change in frequency 
in passing from the liquid to the gaseous state is 
the same for both molecules. Hence we may 
assume that the change in D.S is of the same 
order as in H2S, and take the calculated gaseous 
frequency 1906 cm™ as w}. 

HDS. The peculiar double nature of w2 has 
already been described, and the additional possi- 
bilities discussed should be borne in mind. For 
the purposes of approximate calculation of w3; and 
wi, we may use the valence force expression 
derived by Lechner'® for the asymmetrical tri- 
atomic molecule. The somewhat complicated ex- 
pression is a cubic for the frequencies, and on 
inserting the values Ky,s=3.92, Kp,s=4.09, 
and d(=2k3; as determined below) =0.47 X10 
dynes/cm, together with the observed frequency 
w2=1078 cm™, we find w3=2619 and w,;=1910 
cm in agreement with the values calculated for 
the gas from Dadieu and Kopper’s results. We 
may thus assume that the weak, extended region 
of absorption at 5 u contains w; of HDS together 
with w; of DoS. 3.5 4 (band D’) contains ws; of 
HS, ws of HDS, and w;+we of HDS at 2938 cm=. 
Band D” at 2.64 4 contains 2w; and w3+we of 
HDS, with w3++ws of HS. There are indications 
of 2w2 for HDS at 4.8 u. 


The above conclusions are summarized in 
Table I. 


°F. Lechner, Wien. Ber. 141, 291 (1932). 
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DETERMINATION OF THE FORCE CONSTANTS 


The significance of the force constants will be 
seen from the potential energy expression used by 
Van Vleck and Cross :!! 

V = }[hi(Ar1)?+ke(Are)? + 2kar1” (Ay)? ] 

+ kyAriAre+2?ki3r19(ArsAy+AreAy). (1) 
On the assumption 7,°=72°, and k;=k2, the solu- 
tion of the secular equation yields 


(€11—@12)A3? — ki +ki2=0, 


[ (dir t+ai2)A1, 2 —k, — kyo |[.a331, °—k; | 
— 2[ai3\? — ki3 P=0. 


(2) 


The coefficients a;; are not given explicitly in the 
original. They have therefore been deduced and 
recorded for various values of the vertical angle, 
y (Table II). If m and M are the relative atomic 
masses of H and S, respectively, 


w=m(m+M)/(2m+M), ws=m?/(2m+M). 


Since the coefficients are all of dimensions [ M }', 
we may write \?=47°c*v?=5.863X10-*y*?, and 
thus obtain the k’s in dynes per cm. 

The frequéncy equations apply to the corrected 
frequencies, v;, and not to the observed frequen- 
cies, w;. We have therefore adopted an expression 
of the type 


Eooo=X1VitX2V2t+X2V3 
+Xi1 VP +Xe VP+X33 V3? 
+X 12 V; VotXo2VoV3s+Xa1 Vs Vi, (3) 


where V;=(n;+4), and have attempted a solu- 
tion from the observed data. We have, for ex- 
ample, w100=E100—Eooo when E is expressed in 
cm~!, and the values of X; give the required 
frequencies v;: unfortunately only eight bands 
can be located and assigned with reasonable 
certainty for H2S, and it has been necessary to 
calculate Xj. from the D2S spectrum on the 
assumption that the anharmonic constants in the 


TABLE II. 








H:S 


a; (4 
92°20" 105° 110° ‘o2r20' 105° 110° 


10.971 (0.971 |1.886 see /1.884 





Qu =a22=p—p3* sin? y/2(u—us cos y) 
ai2 =—u3 COS Y | | 

+ (us? sin? y)/2(u—us eos y) | .001 | .008 | .010 |0.008 {0.032 |0.043 
@13=a23= (1/2 y 2)u2 sin y 010 | .010 | .010 | .039 | .038 | .037 
a33 =}(u+ys cos y) 242 .240 | .240 A473 | 465 462 


10 971 





630 BAILEY, THOMPSON AND HALE 


two molecules show the same ratio (determinable 
from the reduced masses) as is found for the 
fundamentals,!” e.g., (w2/we*)? = 1.87. The results 
for H.S are X ,=2719, X= 1294, X3=2749, 
Xiu= — 37, Xo2= — 24, X33= —i1, X 12> +57, 
Xe3= —74, X31= — 114, and for DS, X,=1961, 
X2=929, X3= 2003, Anu= — 20, Xe= —13, X 33 
= —6, X 32= 30, X o3= — 40, X 31= — 62, all in 
cm~!. The high positive value for X12 is in ac- 
cordance with the other observed data: further- 
more we can now apply Teller’s relationship’ 
between the symmetrical frequency 7; for the two 
molecules H2S and D.S and thus make an esti- 
mate of the vertical angle, y. We have (v1/71*)? 
= D(S+2H sin® y/2)/H(S+2D sin? y/2) where 
H, D, and S are the atomic masses of the respec- 
tive atoms, and hence from (2719/1969)? we have 
vy =110°. This value is considerably greater than 
that of 92° 20’ deduced by Cross’ from the rota- 
tional structure of the band at 9911 cm™', but as 
will be seen below, it furnishes better agreement 
among the force constants in the vibrational 
spectrum, and on the whole one would expect the 
vertical angle for H2S to be certainly not less 
than that of the water molecule, which is 105°. 

That the mechanical frequencies so deduced 
cannot be far from true may be seen by solving 
a simplified valence force system for the two 
molecules.'? Without a knowledge of the force 
constants we can construct a cubic in the fre- 
quencies for the vertical angle: 


prp? 3 
ps 1+(2m/M) 





—(pr+pe+ps*)x 


m 
+2p2(14+—) =0; (4) 
M 


where x=1+(2m/M) sin? a, a being the semi- 
vertical angle, and p;=27v;c. On insertion of the 
frequencies deduced for infinitesimal amplitude, 
we find for H2S that 2a is 111° 26’, and for DS 
that 2a is 111° 30’. The uncorrected frequencies 
give no such measure of agreement, the value for 
H2S being 102° and for DeS, 119°. 

We may note in passing that the expression 
deduced by Rosenthal’ for the vibrational effect 

aon Teller, Det kgl. Danske Videnskaternes Selskab 
* rs R. Bailey and A. B. D. Cassie, Proc. Roy. Soc. 


A137, 622 (1932). 
19 J. E. Rosenthal, Phys. Rev. 45, 426 (1934). 


TABLE III. The force constants for H2S and DS (i 
dynes/cmX 10-*). 








ki—kie 


92° 20’ H2S 4.29 4.21 
DS 4.40 4.27 
110° HS 4.26 4.25 
D2S 4.33 4.34 


kitkie ks 























in w2, the asymmetric vibration, is incorrect. As 
given it is 


w3*/w3=([1—2ku(sin? a+ cos? a)" ], (5) 


where k=43Am(m+Am)-'!; p= M/(2m+M); and 
2a is the vertical angle. If the term on the left is 
squared, one obtains results of the right order. 
Miss Rosenthal’s formula for the ratio of the two 
symmetric frequencies is useful: we have 


(w1* /w1)(we*/we) =(1—2kw)*(1—2k)*. (6) 


Inserting v;= 2719, ve=1294, and »,*=1961, we 
have ve*=925 compared with »2*=929 cm™ ob- 
served, and so obtain a good check on our assign- 
ments and observations. 


THE FORCE CONSTANTS 


We have solved the frequency Eq. (2) using the 
calculated values of y for HeS and DS. On com- 
parison with H.O it now seems fairly certain that 
the potential energy curves of a polyatomic mole- 
cule do not correspond sufficiently closely with 
those of the isotopic molecule to enable them to 
be used to give highly accurate values of the 
vertical angle, and in fact no reasonable solution 
for this angle can be obtained on the assumption 
that the force constants are the same for the two 
molecules. We have therefore examined the Van 
Vleck and Cross relationship for certain values of 
the angle, and summarize those for 92° 20’ and 
110° in Table III. 

Although the absolute value for k; must not be 
taken too literally, nevertheless it is significant 
that the larger value for y tends to remove the 
discrepancies between H2S and D,S, and between 
ki—ky2 and ki +k. For HDS, with a33=0.632 
for 110°, the uncorrected frequency 1078 cm" 
gives ks as approximately 0.216X10° dynes per 
per cm as compared with 0.219 for HeS and 0.218 
for D2S. It is interesting to record approximate 





re —__-~= -— 


fa fF © ©. = wm “7 


DECOMPOSITION OF 
values for k; in the series H2O, HeS, and H2Se: 
with uncorrected data these are 7.8, 4.0, and 
3.2 10° dynes/cm. 

For the sake of completeness we will remark 
that Végard®® from an x-ray analysis gave a 
rectilinear structure for H2S, at least in the solid. 
It is noteworthy that the observed frequencies 
give reasonable force constants if introduced into 
equations for a rectilinear model: the peculiar 
shape of the long wave fundamental would then 
have to be explained by invoking the presence of 
an unresolved Q branch; while all observed bands 
other than we and ws; must be interpreted as 
combination tones involving an inactive fre- 
quency. The only band in H2S which cannot be 
so assigned is 2w2 at 2423 cm~, but the last point, 
together with the irregular nature of the rotational 
lines in the band at 2.6 » resolved by Nielsen and 
Barker,’ is sufficient to exclude the rectilinear 
structure. 

While this paper was in preparation, Nielsen 


20 L. Végard, Nature 126, 916 (1930). 


DEUTERIUM IODIDE 631 
and Nielsen*! published a brief note as a pre- 
liminary account of similar work. They confirm 
the existence of the long wave band at 8.0 u for 
H2S, and the corresponding bands for DeS and 
HDS. They do not appear to have located the 
fundamental w:; for D2S; in its place they give a 
band at 4.55 uw and relate it to w3 of HS. This 
wave-length corresponds to a force constant of 
some 5.2X10° dynes/cm, and hence the fre- 
quency is much too high. The strong peak in 
HDS at 988 cm™ also seems to have escaped 
observation. The following comparative table is 
given by Nielsen and Nielsen: 


HDS DS 
1.9 2.1 —_ 
2.65 3.15 3.65u 
3.70 41 4.55 
7.9 9.0 10.8 


H.2S 


We wish to express our acknowledgments to 
Professor F. G. Donnan, F.R.S., for his continued 
interest in our work. 


21 A. H. Nielsen and H. H. Nielsen, J. Chem. Phys. 4, 229 
(1936). 
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The Rate of Thermal Decomposition of Deuterium Iodide 


Joun C. L. BLAGG AND G. M. Murpny,* Chemical Laboratory of Columbia University 
(Received August 7, 1936) 


The partition functions for deuterium iodide have been 
calculated at temperatures from 300 to 1500°K. Constants 
were calculated for an equation to give the value of the 
free energy function at any temperature within this range. 
A similar calculation was made for D2 using the data of 
Johnston and Long. Combining these with the data for 
the other molecules involved, the equilibrium constant for 
the decomposition of deuterium iodide and for the exchange 
reaction was calculated. These results can be differentiated 
and integrated to give values for AH®, AC, and AS®*. The 


HE effect on the rate of several chemical 
reactions produced by substituting deu- 
terium for hydrogen has been studied by various 
investigators.! Because of the simple homo- 
geneous and bimolecular character of the thermal 
decomposition of hydrogen iodide, the corre- 


* Present address, Sterling Chemistry Laboratory, Yale 


University, New Haven, Conn. 
‘For a review of the subject, see Urey and Teal, Rev. 
Mod. Phys. 7, 34 (1935). 


equilibrium constant for the reaction H,+2DI@D.+2HI 
has been plotted against the temperature from 300 to 
1500°K. The equilibrium constant for the reaction DI 
=4D.+3I, has been calculated. The rate of thermal 
decomposition of deuterium iodide was measured between 
660 and 719°K. The experimental ratio of the rate constants 
for the reactions HI—}H2+ I, and DI—}D.+ 31; has 
been compared with the theoretical value of Wheeler, 
Topley and Eyring. Contrary to their expectations, it 
was found to be lower than the theoretical value. 


sponding reaction with deuterium appears to be 
one of the most important reactions to study at 
the present time, since theoretical calculations 
can be made for this case with some certainty. 
We have therefore studied the rate of thermal 
decomposition of deuterium iodide between 660 
and 719°K. 

Since the reverse reaction also occurs, the rate 
of decomposition of hydrogen iodide depends 
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TABLE I. TABLE II. 

DI T Q —(F°—E,)/T 
B, 3.301 300 64.550 43.828 
‘ 6.180102 500 108.585 47.399 
D. — 5.347 X10-5 800 182.750 50.768 
< 1640.307 1000 240.238 52.420 
Xedg 20.043 1500 404.875 55.471 








upon the equilibrium constant. The latter has 
been determined experimentally by Bodenstein? 
who has also made an extensive study of the 
rates of formation and decomposition. Additional 
experiments have been made by Kistiakowsky® 
and others and a critical review of both the 
experimental and theoretical aspects of the 
reaction has been given by Kassel.* The equi- 
librium constant for deuterium iodide has not 
yet been measured but Rittenberg and Urey® 
have investigated the reaction 


H.+2DI—D,.4+2HI (1) 


and have shown that the experimentally de- 
termined equilibrium constants are in agreement 
with those calculated from spectroscopic data. 
Since the energy levels and their statistical 
weights are well known for hydrogen, iodine and 
hydrogen iodide and since the corresponding 
data for the isotopic deuterium compound are 
readily obtainable,* we may proceed with confi- 
dence and use the equilibrium constant for 
deuterium iodide as calculated from the sta- 
tistical theory. This calculation has been made 
recently for hydrogen iodide by Murphy’ who 
has shown that the experimental results are in 
reasonable agreement with the theory. It 
appears, however, that the temperatures at which 
Bodenstein has worked may be in error and an 
attempt has been made to correct them. The 
other difficulty lies in the uncertainty in the heat 
of dissociation of hydrogen iodide which cannot 
be determined from spectroscopic data. We have 
therefore used Bodenstein’s data to determine 


2 Bodenstein, Zeits. f. physik. Chemie 13, 56 (1894); 22, 
1 (1897); 29, 295 (1899). 
3 Kistiakowsky, J. Am. Chem. Soc. 50, 2315 (1928). 
‘Kassel, Kinetics of Homogeneous Gas Reactions (The 
Chemical Catalog Co., New York, 1932), p. 148. 
RYe ied and Urey, J. Am. Chem. Soc. 56, 1885 
6 Jevons, Report on Band-Spectra of Diatomic Molecules 
(The University Press, Cambridge, 1932). 
7 Murphy, J. Chem. Phys. 4, 344 (1936). 








this quantity.* The calculated values of log K are 
plotted against the reciprocal of the absolute 
temperature in Fig. 1 together with the experi- 
mental points of Bodenstein and of Rittenberg 
and Urey. A more certain value of the heat of 
dissociation of hydrogen iodide would only result 
in a parallel shift of the whole curve in either 
direction, but it seems impossible to adjust the 
statistical theory in any way so as to obtain 
agreement with experiment at the highest and 
lowest temperatures studied by Bodenstein. 

The method of calculating the equilibrium 
constants for deuterium iodide is similar to that 
for hydrogen iodide, the partition functions, 


Q=E prexp(—Ee, D/kT] 


being summed by the method of Kassel.° The 
summation is taken over vibrational and rota- 
tional quantum numbers, v and J, py =(2J +1) is 
the quantum weight of the Jth rotational level 
and E is the vibrational and rotational energy of 
the molecule as given by analysis of the molecular 
spectra. The molecular constants were obtained 
from those of hydrogen iodide and are given in 
Table I. Following Giauque and his co-workers, 
the natural constants found in the International 
Critical Tables’ were used. The summation was 
made at 300, 500, 800, 1000 and 1500°K and 
the values for the partition function, Q, are given 
in Table II. The free energy may then be 
obtained by the relation 


— (F°—E,°)/T=(3/2)R1n M 
+(5/2)RiIn T+R In Q—7.267 = (3) 


8 Bodenstein has also neglected to correct for the 
decomposition of iodine molecules into atoms. This has 
been done in the points as plotted in Fig. 1, although it !s 
nearly negligible at his temperatures. The decomposition 
of Hz and D, into atoms is also negligible at these tempera- 
tures. We have therefore neglected these slight corrections 
in our rate constant. 

* Kassel, J. Chem. Phys. 1, 576 (1933); Phys. Rev. 43, 
364 (1933). 

10 Int. Crit. Tab., Vol. I, p. 16. The numerical constant of 
the free energy equation was taken as 7.267; R=1.9869; 
he/k =1.4324. ’ 
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where M is the molecular weight, E,° is the 


energy at the absolute zero and the other. 


quantities have their usual thermodynamic 
meaning. Since the normal state of deuterium 
iodide is 'Y the electronic statistical weight is 
unity. There are no other electronic states to be 
considered. The nuclear spin contributions have 
been omitted which is valid for temperatures 
above 25°C as shown by Gibson and Heitler.!! 
With the use of Eq. (3) and the values of the 
Q's as shown in Table II, we have obtained 
values for the free energy function as given in the 
same table. The atomic weights for deuterium 
and iodine were taken as 2.01363! and 126.932,” 
respectively. For convenience in interpolation, 
the values obtained from Eq. (3) have been fitted 
to a formula of the type 


—(F°—E,.°)/T=a/T+6 log T+cT+dT?+i. (4) 


The coefficients in this equation for deuterium 
iodide together with those obtained for D2 from 
the calculations of Johnston and Long" are 
shown in Table III. An equation of the form of 
(4) will reproduce the value of the free energy 
function to within at least 0.001 cal./mole 
between the temperatures of 298.1 and 1500°K 
as shown by comparison with the results obtained 
by summation at other temperatures within this 
range. It should probably not be extrapolated 
beyond these limits. The coefficients in Eq. (4) 
may now be combined with the corresponding 
constants for Iz as calculated by Murphy,’ to 
give 
—-Rin K=AF°/T=AH,/T—AYo In T 

— ZAP, 7T—23Ar2T?+J7, (5) 


providing the value for AE’ is known. This may 
be obtained from the value 971.83 14 cal./mole 














TABLE III. 
DI D, 
a — 365.512 — 6.92285 
b 11.1051 15.7301 
c 3.642 x 10-3 2.85 10-5 
d —6.78 10-7 1.09 x 1077 
i 16.5063 — 11.1506 











"Gibson and Heitler, Zeits. f. Physik 49, 465 (1928). 
® Int. Crit. Tab., Vol. I, p. 44. 
8 Johnston and Long, J. Chem. Phys. 2, 389 (1934). 
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DEUTERIUM IODIDE 
TABLE IV. 
DI—3}D.4+3lI. 
AH, 562.08 
AT» 2.8930 
AT; ° — 6.8113 x 107 
AT. 4.1889 « 10-6 
I 19.8513 








for the hydrogen iodide reaction and the zero- 
point energies of the various molecules involved, 
by using the equation 


"AE (= "AE +(3E op," 
_ Eopr® sins >Eou,°+Eont’), (6) 


where the superscripts one and two refer to 
reactions (7) and (8) respectively. 


HI-3H2+3l,, (7) 
DI4—-De2+3le (8) 
and E’ = (we. — 1 xe) -heN. (9) 


We have used constants as given by Urey and 
Teal! from the spectroscopic data of Jeppeson'4 
for He and De, obtaining for the zero-point 
energies of He, De, hydrogen iodide and deu- 
terium iodide the values 2175.02 1545.62, 
1144.84 and 815.142 cm“, respectively.'® The 
value of *AE,’ obtained from Eq. (6) is 1014.52 
cal./mole. The principal uncertainty in this 
quantity is that due to 'AE,? and therefore *AE,” 
is probably good to about 14 calories. The 
constants for Eq. (5) are given in Table IV. 
The equilibrium constant for reaction (8) is 
now obtained from Eq. (5), the constants of 


14 Jeppeson, Phys. Rev. 44, 165 (1933). 

16 The vibrational constants used by Giauque (J. Am. 
Chem. Soc. 52, 4816 (1930)) and by Johnston and Long, 
reference 13, which we have used in calculating the 
partition functions for H2 and De give somewhat different 
zero-point energies for these two molecules than the ones 
we give above. The use of their zero-point energies changes 
2AFo° by about 50 cal. but this affects the final constants 
by a negligible amount as we have discovered by calcu- 
lation. It seems better to use the more recent vibrational 
constant for Hz and D2 as determined by Jeppeson. A 
similar conclusion was reached by Davis and Johnston 
(J. Am. Chem. Soc. 56, 1045 (1934)) who report that the 
use of Jeppeson’s data in place of that used by Giauque 
changes the entropy by only 0.001 cal./mole/degree. Teal 
and MacWood (J. Chem. Phys. 3, 760 (1935)) have given 
vibrational constants for H2and Dz, obtained from Raman 
Spectra, which differ slightly from those of Jeppeson and 
which are probably better due to the method in which they 
were obtained. However, the use of these constants instead 
of those of Jeppeson would change AEF,” for reaction (8) by 
only 0.98 cal./mole which is much less than the other 
uncertainties involved. We have taken 1 cm-=2.8461 
cal./mole. 
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Table IV and this value of 7AE,°, and the results 
are shown in Fig. 1. Since the principal uncer- 
tainty in the whole procedure is that in 'AE)’, a 
better value of this would only shift both curves 
parallel to themselves in either direction by an 
equal amount. There are at present no experi- 
mental data to check the deuterium iodide 
curve. 

The well-known methods of thermodynamics'® 
may also be used to give reliable values of A.S®, 
AC,, and AH® from Eq. (4). Although the results 
obtained are not as exact as those for the free 
energy, the precision is sufficient for most 
applications and the tedious task of summing the 
partition functions a second time is unnecessary. 
We obtain in this way, from the constants of 


Table IV 


AF 093.1 = 1159+15 cal. /mole 
and 
AS o98.1 = — 1.303 cal. /mole. 


Using the value Soos.1 (12) = 62.29 as calculated by 
Giauque!? and Sog.1 (D2) = 34.36 as given by 
Johnston and Long,!* we find S293.1 (DI) = 49.62 
cal./mole. The values of AH® and AS® at this 
temperature are higher than the corresponding 


TABLE IVa. 





H.+2DI—D.+2HI 





— 306.80 

2.6378 
AT; — 1.03694 x 10% 
AT: 7.7832 10° 
I 14.3158 


AHy 
AT» 











‘6 Lewis and Randall, Thermodynamics (McGraw-Hill 
Book Co., New York, 1923), Chapter XV, p. 173. 
17 Giauque, J. Am. Chem. Soc. 53, 507 (1931). 
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quantities for hydrogen iodide by 46 and 0.1 
cal./mole, respectively. 

The equilibrium constant for reaction (1) may 
also be calculated from the constants of Table IV 
and smilar constants for reaction (7) as obtained 
by Murphy.’ The resulting coefficients in Eq. (5) 
are found in Table [Va and a plot of K against 
temperature is shown in Fig. 2. The partition 
functions of the four molecules involved vary ina 
very complicated way with the temperature and 
so it is difficult to make any exact statements as 
to the way in which the various factors involved 
contribute to the free energy change of the 
reaction. There is no reason, however, to doubt 
the statistical calculations and, if the equilibrium 
should be measured for this exchange reaction 
over a wide temperature range, this curious 
variation with temperature will probably be 
found. 


RATE OF DECOMPOSITION OF DEUTERIUM IODIDE 


The rate constant for the decomposition of 
deuterium iodide is given by the expression 


ko = (K2/(0.4343-a-1)) log ((a+s)/(a—s)), (10) 


where K¢ is the equilibrium constant for reaction 
(8), @ is the initial concentration of deuterium 
iodide in moles per cc, i is the time expressed in 
seconds, a=2K.2(1—2) and 2z is the fraction of 
deuterium iodide decomposed at time ¢. This is 
essentially the equation used by Bodenstein, 
although expressed in different units and 
rearranged to facilitate computation. As may be 
seen from Eq. (10) and as pointed out by 
Kistiakowsky,’ the errors involved in determining 
ke are smallest when the time of heating is so 
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chosen that a small but uniform decomposition 
occurs as far as possible from the equilibrium 
concentration. Our experiments have been car- 
ried out, for the most part, at the same time for 
each temperature. This is valid for the deuterium 
iodide decomposition since the reaction is known 
to be homogeneous and bimolecular from 
Bodenstein’s work on hydrogen iodide. 


EXPERIMENTAL MATERIALS 


Deuterium was prepared by decomposing 99.9 
percent deuterium oxide in an electrolytic cell so 
designed as to separate the oxygen and deuterium 
as they were generated. The iodine was resub- 
limed before using. Deuterium iodide was 
prepared by passing a mixture of deuterium and 
iodine gases over two platinum spirals, burning 
at red heat, in the manner described by Cook and 
Bates.'® An excess of iodine was used to insure 
almost complete combination of the deuterium. 
The deuterium iodide was purified by distillation 
from a trap surrounded by a dry ice-alcohol 
mixture to a trap immersed in liquid air. Three 
or four distillations were made to insure the 
purity of the deuterium iodide as shown by the 
colorless product which was obtained. The 
deuterium iodide was kept frozen by liquid air 
until the cells were ready to be filled. Hydrogen 
iodide was prepared in a similar manner using 
tank hydrogen, which had been passed over a roll 
of red hot copper screen, and was used to check 
our experimental procedure as described below. 
All chemicals used in the analysis were of high 
grade analytical purity. 


FURNACE AND REACTION VESSELS 


The furnace was a sulfur boiler in which the 
temperature was controlled by keeping the 
pressure over the sulfur constant. The device 
for controlling the pressure was developed in this 
laboratory and is described by Rittenberg and 
Urey® and independently by C. C. Coffin.” 

The temperature was measured with a thermo- 
couple made from No. 28 Chromel-P and No. 28 
Alumel wire and the voltage read with a Leeds 
and Northrup type K_ potentiometer. The 
standard cell used with the potentiometer was a 


'* Cook and Bates, J. Am. Chem. Soc. 57, 1775 (1935). 
" C. C. Coffin, J. Am. Chem. Soc. 55, 3636 (1933). 
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* Runs 1 to 8 inclusive were made with hydrogen iodide. All other 
runs were made with deuterium iodide. 


new Eppley cell which had a negligible tempera- 
ture coefficient. The thermocouple was calibrated 
against the ice point and boiling point of water, 
the boiling point of naphthalene and the boiling 
point of sulfur. These are all approved standards 
of pyrometry.” The boiling points were corrected 
to atmospheric pressure by the equations given 
in the International Critical Tables.” A least 
square equation was calculated for the tempera- 
ture. The temperature remained constant to 
better than 3° since it was impossible to notice 
fluctuations with the potentiometer. 

Quartz reaction cells were used. They were 
approximately six inches long with an outside 
diameter of about one inch. The volumes were 
between 50 and 60 cc. The volume of each cell 
was determined by one of two methods. One, by 
weighing empty and then filled with water, and 
two, by filling with water from a burette. The 
two methods checked to within 0.05 cc. The 
volume was needed in order to determine a. 


20 Int. Crit. Tab., Vol. I, p. 53. 
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EXPERIMENTAL PROCEDURE 


The cells were pumped for from twelve to 
twenty-four hours and heated to a dull red heat 
four or five times during the pumping. When the 
cells were ready to be filled, the liquid air was 
removed from the storage trap containing the 
deuterium iodide and a bath of alcohol dry ice 
was substituted. The amount of dry ice used was 
such that the temperature of the bath was 
between —50 and —45°C as measured by a 
pentane thermometer. This gave a pressure of 
deuterium iodide desirable for our experiments 
and was reproducible as can be seen from column 
four in Table V. The pressure of deuterium 
iodide was allowed to build up for about ten 
minutes and then a stopcock between the storage 
trap and the cells was opened. The gas flowed 
into the cells through capillary stopcocks which 
were closed after about fifteen minutes from the 
time the liquid air was removed from the storage 
trap. The gas in the cells was then frozen with 
liquid air and the cells sealed off from the 
vacuum system. Contact with the stopcock 
grease caused a slight amount of decomposition, 
but since the cells were filled through capillary 
tubes the diffusion of the decomposition products 
was slow and no color was shown by the deu- 
terium iodide in the cells. Bates and Lavin”! have 
shown that stopcocks lubricated with vaseline 
may be used with hydrogen iodide. In these 
experiments Lubriseal was used throughout. The 
cells for runs with hydrogen iodide were filled in 
a similar manner and, since a good check with 
Bodenstein’s results was obtained, any error 
introduced by the use of greased stopcocks is 
negligible. Before each set of runs the deuterium 
iodide was redistilled. 

After a measured length of time in which the 
cells were kept in the constant temperature 
furnace, the reaction was stopped by placing the 
cells in a blast of cold air and then freezing the 
deuterium iodide and iodine with liquid air. The 
total time elapsed between removal from the 
furnace and complete freezing never exceeded two 
minutes. Any error thus introduced was balanced 
out by the time necessary to bring the cell up to 
temperature when placed in the furnace. The 
reaction was timed by an electric clock and was 


*1 Bates and Lavin, J. Am. Chem. Soc. 55, 81 (1933). 
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measured to the nearest second from the time the 
cell was placed in the furnace until it was removed 
and placed in the blast of air. The cells were then 
opened, sealed to the vacuum system and the 
gases collected in a trap filled with a four percent 
solution of potassium iodide which was free from 
dissolved carbon dioxide and which had been 
tested for neutrality and for the presence of 
iodate.”” 

METHOD OF ANALYSIS 


The trap containing the iodine and deuterium 
iodide dissolved in potassium iodide was emptied 
into a flask and thoroughly washed. A standard 
solution of sodium thiosulfate, free from dissolved 
carbon dioxide, was used to titrate the iodine. 
The thiosulfate solution was standardized each 
day by means of a standard ceric sulfate solu- 
tion.”* This latter solution had been standardized 
against Bureau of Standards sodium oxalate, 
using orthophenanthroline ferrous complex, ‘as 
the indicator since a back titration with ferrous 
sulfate was necessary. For iodine-thiosulfate 
titrations a one percent freshly prepared potato 
starch solution was used. Weight burettes were 
used for all titrations and all titrations were 
made in a reproducible light. 

The iodate method for acid determination was 
used in analyzing for deuterium iodide.** After 
analysis for iodine had been made, two grams of 
potassium iodate were added and the solution 
allowed to stand for twenty minutes and then 
titrated with thiosulfate. The titration was 
completed within twenty-five minutes after the 
addition of the iodate. A blank was run over the 
same period of time by adding the same amount 
of iodate to twenty-five cc of the same potassium 
iodide solution in which the deuterium iodide and 
iodine were dissolved. Since the volume of each 
cell was known, the number of moles of deuterium 
iodide at the start could be calculated. The 
fraction decomposed was determined by dividing 
the number of grams of thiosulfate solution used 
to titrate the iodine by the total number of 
grams used in both titrations. 

The validity of our experimental method was 


tested by means of eight experiments on the rate 


 Fales, Inorganic Quantitative Analysis (The Century 
Co., New York, 1925), p. 353. 

*3 Suggested by Dr. G. W. Walden of this laboratory. 
po Kolthoff and Furman, Volumetric Analysis, Vol. II, p. 
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of decomposition of hydrogen iodide and the 
results obtained were in good agreement with 
those of Bodenstein. The data obtained for both 
hydrogen iodide and deuterium iodide will be 
found in Table V. The equilibrium constant, K, 
was calculated for deuterium iodide from Eq. (5) 
and for hydrogen iodide from a similar equation 
given by Murphy.’ The values obtained are 
given in the last column of Table V. Our con- 
stants refer to reactions (7) and (8) and the 
subscripts one and two refer to the hydrogen 
iodide and deuterium iodide reactions, re- 
spectively. All of our experimental rate constants 
will be found plotted in Fig. 3, together with 
those of Bodenstein and Kistiakowsky on hy- 
drogen iodide. 

We have considered all the factors involved in 
calculating our rate constants and each constant 
is good to the average deviation from the mean 
which ranges from two to four percent. This 
precision is as good as that obtained by Boden- 
stein and Kistiakowsky on the hydrogen iodide 
reaction. Any error due to K is due to the 
uncertainty in the zero-point energy and is well 
within the precision which we have claimed. 


DISCUSSION OF RESULTS 


It is now well known that the simple Arrhenius 
equation is not valid for reaction rate constants 
over a very wide temperature range. Moreover, 
there are objections to the simple collisional 
formulas as used in the past. We have, however, 
used the latter theory to derive an equation for 
our results of the form 


log ke= —9333/T+0.5 log T+11.799 (11) 


which may be compared with the similar equation 
given by Kassel‘ for hydrogen iodide 


log ki= —9555/T+0.5 log T+12.311. (12) 


Eq. (11) will reproduce log ke to +0.01 and 
should be regarded merely as an interpolation 
formula, rather than ascribing any reality to the 
heat of activation, number of collisions, tempera- 
ture dependence or cross-sectional molecular 
areas. On the basis of this theory, if we assume 
that the collision cross sections of HI and DI are 
equal and that there is no zero-point energy in 
the activated complex, the ratio of the rate 
constants for the two reactions is given by 


DEUTERIUM IODIDE 
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Ri /Ro= (M2/M;)*e4%0"/R?, (13) 


where M; and M2 are the molecular weights of 
HI and DI and AE, is the difference in their 
zero-point energies and equals 938.4 cal./mole. 
Evaluation of Eq. (13) at 700°K gives 1.97 as the 
ratio of the rate constants. Using Eq. (11) based 
upon our experimental results and the corre- 
sponding constant for HI as determined by 
Bodenstein, we find the experimentally observed 
ratio at 700°K of 


ky /Ro=1.156/0.757 = 1.528. 


The definitely slower rate for DI seems to 
indicate that some zero-point in energy is 
required in the activated state. 

Eyring” and his collaborators have calculated 
the absolute rate of reactions in terms of potential 
surfaces, which in general are those obtained 
from the Morse formula. Statistical considera- 
tions enable the probability of the activated 
state to be calculated and symmetry consider- 
ations furnish considerable information about its 
potential and vibrational properties. Later 
developments of this theory by Wheeler, Topley 
and Eyring’® have resulted in numerical results 
for the reaction rates between the halogens and 
hydrogen and deuterium. They have predicted 
for the rate of formation of hydrogen and 
deuterium iodides, the values k;'/k3’=2.6, 2.34 
and 2.26 at 575, 700 and 781°K where the primes 
on the rate constants refer to the reverse reaction 


25 Eyring, J. Chem. Phys. 3, 107 (1935). 
= Topley and Eyring, J. Chem. Phys. 4, 178 
1936). 
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of that measured by us. Using the equilibrium 
constant as calculated above we find for 700°K, 
which is in our temperature range, 


ki /ko= (K,/K2)?: ky’ /Re! =(0.818 2.34= 1.92 


according to the theory of Wheeler, Topley and 
Eyring. This result is based upon 14 percent 
additive binding in the potential surfaces, the 
zero-point energy of the activated state being 
ignored. Approximately 20 percent is required to 
give the experimental result for the HI reaction, 
if activated zero-point energy is_ included. 
However, they have used Bodenstein’s experi- 
mental results to determine the proper fre- 
quencies, activation energy and moments of 
_inertia and since the DI reaction uses the same 
potential surface, their result should give the 
correct constant for the DI reaction. Because of 
the difficulties involved in the theoretical calcu- 
lations, however, we feel the agreement with our 
experimental ratio of 1.528 is highly satisfactory. 
Unfortunately, Wheeler, Topley and Eyring 
expected the observed ratios to be somewhat 
larger than the value cited above, since they 
believe their calculated frequencies for the 
activated complex to be too stiff. Since the 
experimental ratio for this reaction has been 
highly desirable as a check of the theory, it will 
no doubt suggest to the theoretical chemists the 
proper way to adjust the potential surfaces in 


DUNCAN 


order to obtain closer agreement. 

Recently, Geib and Lendle?’ have presented 
some results on the rate of formation of DI. 
Their ratio at 700°K is 2.45 and converting 
to the decomposition reaction as before by means 
of the equilibrium constant, we find k;/k.=2.00. 
This is considerably higher than our result and 
also somewhat higher than the theoretical one of 
Wheeler, Topley and Eyring. They followed the 
course of the reaction by measuring the extinc- 
tion coefficient of iodine by a photometric 
method. Since they did not use monochromatic 
light, rather large corrections were necessary for 
the effect of H: and HI on the calibration curves 
of the iodine absorption. Their precision is 
therefore rather low and since their agreement 
with Bodenstein’s results for HI is not very 
satisfactory, the individual rate constants, as 
well as their ratio, may be in error. 
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The absorption spectrum of nitrous oxide below 2200A has been reinvestigated. Much new 
material has been found, particularly below 1550A, including some discrete systems of bands, 
and some continuous regions of absorption. No bands are found below 997A, but only con- 
tinous absorption down to about 850A, the limit of observation. A Rydberg series was found, 
converging to an ionization potential of 12.66 volts. Evidence was found for a total of ten 
electronic states including the members of the Rydberg series. The vibration frequency in the 


electronic transition around 1500A is 621 cm7. 


HE near ultraviolet absorption spectrum of 
N:O has been studied by many investi- 
gators.'~* The absorption above 2000A appears 
‘A. K. Dutta, Proc. Rov. Soc. A138, 84 (1932). 
2Q. R. Wulf and E. H. Melvin, Phys. Rev. 39, 


(1932). 
’L. Henry, Nature 134, 498 (1934). 
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to be continuous, and very weak. Below 2000A 
Leifson‘ found two broad bands, one whose center 
was apparently about 1850A and the other ex- 
tending from about 1550A to the limit of his 
observation. These regions have also been studied 


4S. Leifson, Astrophys. J. 63, 73 (1926). 
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Fic. 1. (a) source, (b) NsO0—0.001 mm, (c) —N2.O—0.01 mm, (d) —N»O—0.02 mm, (e) N,O—0.04 mm, 
(f) NxO—0.075, (g) NLO—0.09 mm, (h) NXO—0.19 mm. 


by Sen-Gupta’ who has discussed the relation of 
this absorption to the decomposition of N2O into 
NO and N. The present writer has investigated 
the absorption of this molecule from 2200A 
down to about 850A and found several new 
regions of absorption, some of them showing 
vibrational structure. 


EXPERIMENTAL 


The NO (a tank product prepared for 
anesthetic work) was purified by bubbling slowly 
through a tower filled with strong KOH solution, 
then dried by passing through P2,O;. It was 
collected by liquid air in a large bulb which was 
attached to the spectrograph. It was then further 
purified by distillation with liquid air. 

The spectra were photographed with a 120,000- 
line, 1-meter focus, glass grating used at normal 
incidence. The dispersion was about 8.5A/mm 
in the first order. The light source was the 
Lyman continuum with He used as the con- 
ducting gas in most cases. Helium, purified with 
charcoal and liquid air, was used in a few ex- 
posures in order to try to improve the Lyman 
continuum in the region 1150 to 900A where 
many emission and absorption lines occur in the 
source. The use of helium leads to no improve- 
ment, which might have been anticipated, since 
most of the lines come from the glass vaporized 
in the discharge. In spite of this the strong bands 
below 1100 could be located without difficulty, 
but weak bands in this region may have been 
missed. 


*P. K. Sen-Gupta, Nature 136, 513 (1935). 


The pressure of N2O was varied from 0.001 
to 0.53 mm in this work. The spectrograph 
served as the absorbing column, giving a two- 
meter path. At these pressures no absorption 
whatever appeared at wave-lengths longer than 
1530A. In one exposure the spectrograph was 
filled to a pressure of 7.3 mm by separating the 
spectrograph from the discharge with a fluorite 
window.® The region found by Leifson and Sen- 
Gupta about 1850 appeared at this pressure, 
extending from 1768 to 2068A. This region was 
not investigated further. 


RESULTS 


The general appearance of the spectrum 
between 1550 and 1050A is shown in Fig. 1 
which is an enlargement (about four diameters) 
of a representative plate. 

The first electronic transition (called C) below 
1550A consists of narrow bands shaded to the 
red, which widen rapidly with increase in pres- 
sure. They first appear at about 0.01 mm. 
Above 0.1 mm the bands merge in a continuum 
whose limits appear to increase slowly with 
increasing pressure. At 0.53 absorption extends 
up to 1528A. At 7.3 mm strong absorption 
extends to 1700A, and the region between 1700 
and 1768 is partly absorbed. Thus the “long 
wave limit’’ of this region given by Sen-Gupta® 
(1580A) is far from constant.’ The apparent 

6 A pressure greater than about 1 mm cannot be used in 
this spectrograph, because of diffusion through the slit 
into the discharge. 

7 The absorption limits in this paper are given solely for 


the practical use of photochemists who may wish to study 
the decomposition or other reactions of N»O by light in 



















































TABLE I. Bands of electronic transition C. 














Obs Calc.—obs. 
65939 65912 +27 
66548 66569 —21 
67135 67166 —31 
67699 67697 + 2 
68239 68224 +15 
68756 68740 +16 
69251 69232 +19 
69722 69745 —23 

70148 +22 


70170 


centers of nine bands observed fit a formula 
v= 65,939+621.2n—11.54n? (n=0,1,2---). 


The frequencies of these bands and the agree- 
ment with values calculated from the above 
equation may be seen in Table I. The bands do 
not appear to be very sharp but this does not 
imply predissociation here, since the rotational 
structure could not be resolved in any case with 
our instrument. The maximum of intensity 
occurs at the band 68,224 cm“. 

Since no discrete absorption in the ultraviolet 
had been reported for NO, we investigated the 
possibility that these bands might be due to N» 
or to NO. We had some plates of pure No» at 
pressures about 100 times that necessary to 
develop these N2O bands clearly. Nz showed no 
absorption whatever in this region. NO present 
in the original gas would have been removed, 
but some may have been formed during ex- 
posures.?~* We then tried irradiating N.O in the 
spectrograph for 30 minutes (longer than the 
longest exposure) with the Lyman continuum. 
The absorption spectrum was then taken and 
showed no additional bands, while the original 
bands were unchanged in appearance. The fre- 
quency separation is also too low for any of the 
known bands of NO in the Schumann region. 
We must conclude, therefore that these bands 
are due to N.O. 
these regions. It would appear that the correlation of long 
wave limits of continuous absorption with heats of dis- 
sociation must be carried out with great caution, for only 
in extremely fortunate circumstances would the relative 
positions of the Franck-Condon curves be such that an 
experimentally determined limit would have any _theo- 
retical significance. In addition, Boltzmann factors for the 
normal vibrations of NO, while small, do not exclude 
absorption from higher ground levels than the zero level, 
which further confuses the interpretation of a long wave 
limit. The ‘‘limit”’ for the next absorption region at longer 
wave-lengths was given as 1850A by Sen-Gupta. Here (at 
7.3 mm 200 cm column) absorption extends visibly to 


2068A. 
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We are unable to relate definitely the single 
observed excited state frequency (621) to the 
normal frequencies observed in the infrared 
It is probably a modification of 1285.4 (N N O) 


— oe = 
or 589( TNNJ| OT ). If the molecule remains 
linear in the excited state, 621 is more probably 
a lowering of 1285 than a modification of 589, 
which is not totally symmetrical. If the molecule 
becomes bent in the upper state no decision 
could be made from symmetry considerations. 
Resolution of the rotational structure of the 
bands would, of course solve the problem, since 
one of these frequencies is || and the other L. 
These remarks apply also to a less complete set 
of bands farther in the ultraviolet where the 
frequency separations have comparable values. 

The next electronic transition (D) farther in 
the ultraviolet shows no vibrational structure at 
the lowest pressures. It is first observed at a 
pressure of 0.002 mm, about 200 cm™! wide, 
with a center at about 77,900 cm~!. The region 
widens almost symmetrically until at 0.2 mm 
the limits are 82,221—73,724 cm“. 

At the end of this region of absorption, a 
region of relatively high transmission begins. 
At 0.05 mm it extends from 82,300 to 83,700 
cm™~!. At 0.2 mm it can still be seen though it is 
partly absorbed at this pressure. 

Between 83,700 and about 93,000, 
another transmission region begins, there are 
probably two electronic transitions. Correspond- 


where 


ing to one there is a strong continuous region of 
absorption very similar to the transition called D. 
The center of this region is approximately 91,200 
at the lowest pressure. This transition is called /. 
To the long wave-length side of this continuous 
region are five sharp bands which seem to be 
related to each other but not to the continuous 
absorption. There are also three other bands in 
this region which are more diffuse and weaker, 
at least in comparison with the continuous 
absorption. The first five bands are considered 
to belong to one electronic transition, called £. 
The last three bands 86,750, 87,439, and 88,241 
could not be explained in any way, and may 
belong to still other electronic transitions. The 
two frequency differences between these three 
bands are so different that a classification with 
the E bands is not permitted. These three bands 
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may be seen faintly on Fig. 1, being almost lost 
in the continuum. The second member of the E 
bands 84,992 cm is the strongest of all the 
ultraviolet absorption bands. As will be noted in 
Fig. 1 this is the only band appearing on the 
lowest pressure spectrum (6). The first band, 
84,397, is much weaker and we believe that it 
comes from the level v’’=1 of the normal elec- 
tronic state. The difference between this and 
the strong band 84,992 is 595, which is equal, 
within our experimental error, to the funda- 
mental perpendicular infrared frequency 589 
cm~'.® The other three bands of the £& transition 
are 85,473, 85,880, and 86,125. The frequency 
differences resemble those in the C transition 
around 1500A,- but are lower. The bands are 
shaded more definitely to the red than the C 
bands. 

Below these bands a short transmission region 
(seen at the left of Fig. 1) extends to about 
94,500 cm~!. There are five bands of moderate 
intensity at shorter wave-lengths beyond this 
transmission region. There is also continuous 
absorption all over this region. Beyond the last 
band measured, 100,397 cm~', continuous ab- 
sorption extends to the end of the plate, about 
850A. The centers of four of these five bands, 
together with the lowest pressure centers of D 
and F fit into a Rydberg series with the following 
formula: 


v=102,567—R/(n—0.92)? (n=3,4,5---). 


There are probably a few vibrational transitions 
accompanying these electronic transitions, but 
they are too weak to measure with certainty. 
These bands with values calculated from the 
above equation are shown in Table II. The 


TABLE II. Rydberg series in N20. 








Calc. Obs. 


~ 
= 


Calc.—obs. 


77,202 77,900 (D) — 700 
90,999 91,200 (F) — 200 
95,975 95,852 +123 
98,315 98,326 — il 
99,598 99,609 — il 
100,378 100,379 - 1 
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TABLE III. Summary of electronic states of NO. 








Lowest pressure at 
a which it appears 
Description (mm) 








center about 1900A, 7.3 mm 
discrete bands, origin 65,939 
max. intensity 68,224 
continuous, center 77,900 
discrete bands, origin 84,992 
continuous, center 91,200 
single band 94,681 
Four others in Rydberg series (Table II) 


between 0.5 and 7.3 
0.01 








ionization potential obtained by extrapolating 
this series, 12.66 ev, is a ‘‘vertical’’ (Franck- 
Condon) ionization potential, since the most 
probable transition to the upper state is measured 
in all cases. This is the center of absorption at 
the lowest pressure at which any absorption can 
be observed. This value may be compared with 
the value 12.9+0.5 ev obtained by electron 
impact® and with the vertical J of 2p7 of 
N(12.7 ev).!° 

Finally the electronic transitions which have 
been found are summarized in Table III. The 
band 94,681 cm! must represent a separate 
electronic transition, but no explanation can be 
given at present of its significance. In reference 
to Table III, the letter A is reserved for the 
absorption above 2100A'* though this may 
prove to be a continuation of B. 

The strongest objection to the Rydberg series 
above comes from the relatively low intensity of 
n=6 as compared to m=5 and n=7. This is 
however evident only when the three are com- 
pared at the lowest pressure. At pressures of 
0.01 mm the intensities appear about the same. 
The strongest band in the whole spectrum does 
not fit into this Rydberg series. It does not fit 
with any other member into a series which con- 
verges to the first ionization potential. It is 
possible however, that it is the first member of 
a series converging to a higher ionization 
potential of this molecule, which perhaps lies 
at wave-length shorter than 850A. 

The writer expresses his gratitude to the 
Warren Fund of the American Academy of Arts 
and Sciences for financial aid during the progress 
of this work. 
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(1) The order of the reaction causing the rapid dis- 
appearance of OH radicals after interrupting an electric 
discharge through H.O vapor is measured by photometry 
of absorption spectra taken in snapshots. (2) These 
spectra are taken in time intervals of a few tenths of a 
second with a 21-foot grating and rotating interrupter 
and sector. (3) Relative concentrations of OH radicals 
are derived by comparing intensities of rotational lines 
within one branch of the 0-0 band of OH, the relative 
intensities being given by the theory of Hill and Van Vleck. 
(4) The decay of OH with time for various H2O pressures 
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and addition of He is given in diagrams and a table of 
rate constants. (5) A discussion of the kinetics leads to 
triple collisions consuming OH. (6) KCI covering the 
surface of the absorption tube materially increases the 
rate of the reaction. (7) The absorption spectrum of HO, 
investigated simultaneously showed that no appreciable 
H2O2 is formed in the reaction. (8) In the processes con- 
suming OH radicals after interrupting the discharge 
through H;O vapor in clean glass tubes, the triple collision 
H+OH+M-—H:0+M plays a predominant part. 








I. PROBLEM AND PREVIOUS WoRK 


N preceding papers! observations have been 
reported of the absorption spectrum of free 
OH radicals in a reacting gas at room tempera- 
ture. In particular it was shown that OH as 
produced by the electric discharge through H2O 
vapor could be observed to last for at least § 
second after interrupting the discharge. It was 
concluded that the recombination process de- 
termining the disappearance of OH is not so fast 
that a bimolecular association has to be assumed, 
as had been suggested by other observations.’ 
The work here reported was carried out in 
order to obtain further information regarding 
this process, and consists in determining quanti- 
tatively relative values of the OH concentration 
and so obtaining the rate of disappearance of OH 
as a function of the composition of the reacting 
gas mixture. 


II. EXPERIMENTAL 


The apparatus was essentially that described 
previously, but improved in several respects in 
order to allow for more quantitative results. It 
consisted of : 


(a) Source of continuous radiation 

(b) Absorption tube _ 

(c) Rotating sector and switch 

(d) Spectrograph 
( ee J. Chem. Phys. 2, 713 (1934) and 3, 266 
1935 


? W. H. Rodebush and M. H. Wahl,"J. Chem. Phys. 1, 
696 (1933). 
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The OH 3064 band was to be observed in 
absorption. The carbon arc in nitrogen as used 
previously did not provide a sufficiently con- 
tinuous radiation to be used as a source in 
quantitative work. However, the carbon arc 
burning in flowing helium was found to be 
satisfactory in this respect and was used in a 
number of experiments. This source was dis- 
carded later in favor of a particularly intense 
hydrogen discharge tube.* The latter carried a 
current of one to two amperes through a 1X5 mm 
capillary and was less intense by a factor of two 
to four, but was much easier to operate, and 
provided a radiation which was constant over 
long periods of time. 

The absorption tube was 150 cm long and 
4.7 cm in diameter and closed at each end with 
quartz windows. Large aluminum electrodes were 
placed in side tubes at a distance of about 20 cm 
from the absorption tube proper. 

A synchronous motor operated a revolving 
sector which allowed light to pass to the spec- 
trograph during the desired interval of time. 
On the same shaft with the sector was a sliding 
contactor which served to control the discharge 
through the gas in the absorption tube. At first 
this contactor was placed in series with the dis- 
charge and the secondary of a transformer, an 
arrangement which proved unsatisfactory be- 
cause of the high potential involved. This 


3 Constructed by N. D. Smith, Thesis, Harvard Uni- 


versity. To. be published. 
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TABLE I. 


KINETICS OF OH RADICALS 








Relative Absorption 





Rotational 
Quantum Number Coefficient 
13 0.571 
7 0.725 
3 0.605 
4 0.390 
3 0.195 





2 0.081 





problem was solved by putting the contactor in 
the grid circuit of a thyratron (FG—29) which 
regulated the current through the transformer 
primary.‘ 60-cycle a.c. 110 and 220 volts was 
used in the primary, the power varying from 2 to 
4 kilowatts. One or four half-cycles were passed 
through the primary depending upon the gear 
system used with the synchronous motor. The 
shaft was rotated at 6 revolutions per second, or 
at 1% revolutions per second, different speeds 
being necessary in order to follow the reaction 
under widely differing conditions. 

A high frequency discharge was used in several 
runs and was produced with the aid of a 250- 
watt power oscillator, the high frequency po- 
tential being applied between external electrodes 
running the full length of the absorption tube. 
The advantage of this arrangement was that a 
discharge could easily be obtained with much 
higher pressures of water vapor, up to 5 mm Hg. 

For obtaining the weak band absorption, the 
second order of a 21-foot concave grating spec- 
trograph was used. The high resolving power 
was necessary in order to observe the very 
narrow absorption lines of the band. To decrease 
the exposure time as much as possible a cylin- 
drical lens was placed in front of the photographic 
plate.» This increased the intensity over the 
small spectral range to be observed by a factor of 
about 20 without affecting the resolving power. 
Eastman IV-—O plates were used. 


III. EVALUATION OF RELATIVE OH 
CONCENTRATIONS 


The individual OH absorption lines being so 
very narrow (cf. the spectrum reproduced, 
reference 1, Fig. 2) it was not found possible to 
measure the intensity of absorption with the 
available microphotometer. However, the follow- 


_ 





* Following a suggestion by Dr. F. V. Hunt. 
°O. Oldenberg, J. Opt. Soc. Am. 22, 441 (1932). 
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ing method involving only visual comparison of 

absorption lines was found to be reasonably 
satisfactory. It so happens that instead of the 
conventional intensity marks, the individual lines 
of one branch of a band represent concentration 
marks since the relative concentrations of the 
rotational quanta for a given temperature and 
their absorption coefficients are theoretically 
predictable on the basis of the theory of Hill 
and Van Vleck.® Lines of the Q; branch were 
used because of their intensity and freedom from 
overlapping with lines of other branches. Assum- 
ing that the rotation of OH is in equilibrium at 
slightly above room temperature (300°K) the 
relative values of the absorption coefficients 
(Boltzmann factor Xstatistical weight x proba- 
bility of transition) are given in Table I. 

In order to obtain relative OH concentrations 
corresponding to two different exposures, the 
exposures were placed side by side-in the field 
of view of a microscope and the absorption lines 
compared until one or more pairs of lines were 
found, the lines of each pair having nearly the 
same intensity of absorption.’ In this comparison 
care must be taken to have equal intensity of 
the continuous background in the two exposures 
and to use the photographic plate in the linear 
portion of its characteristic. One exposure is 
taken as a standard furnishing an arbitrary unit 
pressure of OH. All measurements discussed in 
the following sections are expressed in the same 
standard units of OH pressure. Suppose the mth 
line of the unknown is paired with the mth line of 
the standard and let km and k, be the respective 
values of the relative absorption coefficients 
found in Table I. Then the OH concentration for 
the unknown is k,/km units. Several such con- 
centration values may be obtained for a single 
exposure if interpolation of intensities between 
lines is made use of. Only the third, fourth, 
fifth and sixth lines and two or three gradations 
between each neighboring pair were used in the 
comparison. The available concentration range 
was about 1 : 100. A typical set of values of OH 


6 E. Hill and J. H. Van Vleck, Phys. Rev. 32, 250 (1928). 
Calculations for the special case under discussion were 
carried out by L. T. Earls, Phys. Rev. 48, 423 (1935). We 
appreciate the advice of Professor Van Vleck in the 
problem of intensities. 

7 An instrument was devised which enabled this com- 
parison to be made without the need of cutting the plates 
up into their separate exposures. 
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Fic. 1. Dependence of 1/[OH] upon time; X, HO 1.1 
mm; O, H,O 1.1 mm; He 3.7 mm; 1/[OH] in 1/standard 
pressure. 


° 


concentration obtained in this way for a single 
exposure is as follows: 0.64, 0.50, 0.50, 0.37, 
0.46; average 0.49; mean deviation 12 percent. 
The results are not very precise but they are 
accurate enough for the present purpose of 
determining the order of the reaction. 


IV. EXPERIMENTAL RESULTS 


As is shown in the next section it is expected 
that for any set of conditions there should be a 
linear relation between 1/[OH] and the time. 
That this is true within the experimental error 
is shown in Fig. 1, in which are plotted the results 
for typical cases. The slope of the straight line 
represents an experimental rate constant k for 
the particular conditions employed. Table II is a 
summary of the experimental results as expressed 
by values of k corresponding to different values 
of H,O and He pressure. 


V. CHEMICAL KINETICS 


The various possible processes leading to a 
long life of OH in the electric discharge through 
water vapor have been discussed at length in 
the preceding paper. For the present purpose we 
shall base our calculations on the following 
assumptions : 

(a) The reacting mixture consists of H2O, H 
and OH, with the addition of He in several runs. 
Ions will have recombined within a few thou- 
sandths of a second after the discharge. O atoms 


can be assumed to be present in negligibly small 
amounts if the current is not too high. The ex- 
perimental arguments for this assumption have 
been summarized by Oldenberg. From the 
molecular structure it is plausible that the 
easiest dissociation process leads to H and OH, 
not to free O atoms. 
(b) The probable reactions are: 


(1) H+H+M—H.4+M \ Termolecular Processes 
(2) H+OH+M—-H.0+M } where M is any third 
(3) OH+OH+M—H,0.+M)} molecule 


or instead of (2) and (3) 


(2a) H+OH—H:0 


(3a) OH+OH—H.0, \Bimolecula Processes 


(c) HO is only to a small part dissociated so 
that [ M] will be approximately constant during 
the reaction. 

For reactions (1), (2), (3) the rate of dis- 
appearance of OH is given by 


~ (d[OH ]/dt) = ks [OH TH ][M]+2,(0H PL] 


plus additional terms with different constants if 
more than one kind of molecule is effective as a 
third body. Since we assume the statement of 
Rodebush and Wahl that [H ] and [OH ] are the 
primary products, the reaction starts from equal 
concentrations [H ] and [OH ]. If we furthermore 
assume that H and OH disappear by triple 
collisions at not very different rates it follows 
that [H ]=[.OH ] during the time through which 











TABLE II. 
k 
Plate H,0 He (standard pressure 
number (mm Hg) (mm Hg) units sec.~) 

L3 0 24 

1 2.6 0 24 
3.2 0 48 

1.0 0 14 

2 2.0 0 17 
4.1 0 30 

3 1.0 0 13 
1.0 ta 31 

4 1.45 0 17 
1.45 2.75 27 

5 1.1 0 17 
1.1 3.7 31 

6 0.75 0 18 
0.75 4.1 27 
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KINETICS OF OH RADICALS 


the reaction is observed. In this case the two 
terms may be combined and 


— (d[OH ]/dt) =k[OH #[M], 


which may be integrated, M being constant with 
respect to the time: 


(1/[OH],) —(1/[OH}) =k}. 


The corresponding treatment for OH reacting 
according to (2a) and (3a) would result in 


(1/[OH ].) — (1/[OH ]o) =k’t. 


The same equation would hold if the wall re- 
action is predominant. In any case there should 
be a linear relationship between 1/[OH] and 
the time. This relationship has been used in the 
preceding section. 

The experimental constant k is to be compared 
with the constants k{.M] or k’ of the theory. 
If k is independent of the H,O and He concentra- 
tions the reaction proceeds according to (2a), 
(3a). On the other hand, if k increases with the 
H,O and He concentrations in such a way that 


k= ku,o[H20]+kue[ He], 


then the reaction is termolecular according to 
(2) and (3). That the latter is the case may be 
seen at least qualitatively in Fig. 1 where it is 
noticed that the addition of He increases the 
rate of the reaction. The constant ky,o may be 
calculated directly from data on the reaction in 
the absence of He. Fig. 2 is a plot of k against 
H.O pressure based on the experiments with a 
wide variation of H.O up to 5.2 mm (plates 1 
and 2 of Table II). The rate does definitely in- 
crease with increasing [H2O]. The fact that the 
best straight line does not go through the origin 
seems to indicate the simultaneous occurrence 
of a slow bimolecular reaction, presumably at 
the surface. More evidence of a surface reaction 
is given in the next section. 

The value of ku,o as given by the slope of 
the line of Fig. 2 is 


ku,o=6.8; HO in mm Hg. 


The constant ky. is most easily obtained by 
dividing the increase in the constant k due to 
the addition of He by the pressure of He added. 
The four last plates of Table II give the four 
values 2.5, 3.6, 3.8, 2.2, the average of which is 














4 
(H,oJ mm Hg 


Fic. 2. Dependence of experimental rate constant k 


upon H,O pressure; & in standard pressure units” 
Meee. 


ky.=3.0. These results for kuy.o and ky. are 
based on the same ‘‘standard units” of [OH ]. 

In order to determine whether or not these 
are reasonable values for triple collision rate 
constants the absolute concentration of OH 
must be known or estimated. If we assume that 
the maximum OH concentration obtained with 
the present apparatus, defined as a standard in 
section III, is about 0.1 mm so that the large 
majority of molecules is not dissociated, then 
the reaction constants are 


ky,0 = 68 and ky.=30 
or ku,o =2.4 x 10'° 


(mm-~? sec.~!) 


and ky. =1.1X 10" 
(cc-? mole*? sec.~'). 


These values although being very approximate 
compare favorably with rate constants for similar 
third order reactions such as 


H+H+M—-H.4+ / 
and Br+Br+ M—-Br.+ M° 


with values ranging from 10-10". 

Obviously the problem remains of determining 
the absolute concentration of OH which in this 
argument is rather arbitrarily assumed. The only 

8H. M. Smallwood, J. Am. Chem. Soc. 51, 1985 (1929) 
and 56, 1542 (1934). 

® K. Hilferding and W. Steiner, Zeits. f. physik. Chemie 


B30, 399 (1935); E. Rabinowitsch, Zeits. f. physik. Chemie 
B33, 275 (1936). 
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method to apply is a comparison of the OH 
observed at the discharge with a known concen- 
tration as it was observed by Bonhoeffer and 
Reichardt’ in water vapor partially dissociated 
at temperatures above 1150°C. An experimental 
difficulty turns up, however, in that the two 
absorption experiments to be compared are 
carried through at temperatures different by a 
factor 5 and pressures by a factor 1000, so that 
the line widths are quite different. But the total 
energy absorbed within a spectral line of a 
certain coefficient of absorption depends upon 
the width of the line—apart from the properties 
of the spectrograph. Hence the comparison 
requires a special study of the pressure broaden- 
ing of the OH lines. For this purpose even the 
second order of the 21-foot grating is not 
adequate. An attempt will be made to solve this 
problem and so determine the absolute values of 
the reaction rate constants ky,o and kue. 


VI. SURFACE REACTION 


The pressure effect reported definitely indi- 
cates a gas reaction. If the surface reaction were 
prominent He would decrease rather than in- 
crease the rate, contrary to the values of Table 
II. In the following experiments a potassium 
chloride surface was studied for various reasons. 
KCI is known not to catalyze.the recombination 
of hydrogen atoms. Hence any positive effect in 
the discharge through water vapor must be due 
to the adsorption of OH radicals. A very strong 
effect of KCl was discovered by Pease"! in the 
slow combustion of He at temperatures of 520— 
550°C. This is known to be a chain reaction. 
When the glass walls are covered with KCl the 
rate is reduced by a factor of about 1000. Since 
chain reactions are easily broken by surface 
reactions removing one chain carrier, the effect 
of a KCl surface on radicals is of interest. 
Taylor and Lavin observed the heating of 
thermometers covered with various substances 
in the flowing gases pumped out of the discharge 
through H.O; they inferred adsorption and 
recombination of OH on a KCI surface. 


10K. F. Bonhoeffer and H. Reichardt, Zeits. f. physik. 
Chemie 139, 75 (1928). 

1 R, N. Pease, J. Am. Chem. Soc. 52, 5106 (1930). 

2H. S. Taylor and G. I. Lavin, J. Am. Chem. Soc. 52, 
1910 (1930). 


The rate of disappearance of OH, observed by 
the decay of the absorption spectra, was com- 
pared for a glass surface carefully cleaned with 
cleaning solution, boiling water and distilled 
water, and on the other hand for the same 
surface covered with KCl. In both cases, of 
course, the surface was moist. The KCI surface 
increases the rate materially. While for the glass 
surface the absorption spectrum of OH is easily 
traced through 8/15 sec., it disappears com- 
pletely for the KCl surface in the same time 
interval. Furthermore, since the discharge proper 
producing OH is not instantaneous but lasts 
through 1/15 sec. it is noticeable that, right 
after interrupting it, more OH is present in 
the case of the clean glass than in the case of the 
KCI surface. 

The obvious interpretation of this experiment 
is that the KCI surface adsorbs OH radicals and 
so catalyzes their reaction, in agreement with 
the result of Taylor and Lavin. 

The noticeable increase of the surface reaction 
caused by the KCl surface indicates that for 
clean glass the*surface reaction is not pre- 
dominant. Whether or not in the experiment of 
Pease the same effect takes place, i.e., whether 
OH is one of the chain carriers, remains uncer- 
tain. An endeavor will be made to obtain further 
evidence. 


VII. ABSORPTION SPECTRUM OF H,O, 


In section 5 it was shown that OH disappears 
largely by triple collisions. It remained uncertain 
whether HO or H2Oz is formed, that is, whether 
reaction (2) or (3) takes place. The results of 
Rodebush and Wahl seemed to indicate an 
overwhelming probability for the formation of 
HO: since concentrations of nearly 50 percent 
H.O:2 were obtained in the gases pumped out of 
the discharge. It should be possible to measure a 
concentration as high as that by the absorption 
spectrum. Rodebush and Wahl failed in doing 
this. They explained that their failure was 
caused by lack of sensitivity of the experiment. 
Since the absorption spectrum of OH proved to 
be valuable for the study of kinetics, it seemed 
worth while to improve the sensitivity with 
which the absorption spectrum of HO» was to be 
observed. 


KINETICS OF OH RADICALS 


According to Urey, Dawsey and Rice,!* H2O2 
has a continuous absorption spectrum extending 
from about 3100 with increasing intensity to 
short waves beyond the range of the quartz 
spectrograph. No other molecule present in this 
experiment shows continuous absorption in the 
same range. In particular OH has the major 
intensity in the 0-0 band so that no continuous 
OH spectrum beyond a convergence limit of 
vibrational quanta is to be expected which 
might obscure the H2O2 spectrum. The difficulties 
of the two absorption experiments are entirely 
different. While for the absorption spectrum of 
OH because of the exceedingly narrow lines high 
resolving power is necessary, the absorption 
spectrum of HO. can be observed with any 
small quartz spectrograph because of its con- 
tinuous nature. In order to detect small traces of 
H.Oe2 a very accurate measurement of the light 
intensity passing through the vapor is necessary, 
hence a very steady source of light is required. 

In order to take advantage of the strong 
absorption at short wave-lengths, some experi- 
ments were made with an electric, discharge 
through Cd vapor, the resonance line of which at 
\2280 is rather isolated in its spectrum. However, 
the mercury line \2537 emitted from an ordinary 
arc proved about as effective and easier to apply. 
Instead of using a monochromator in order to 
separate the line \2537 from the intense long 
wave-length spectrum which is not absorbed by 
H,O2, a chlorine filter was applied, consisting 
of a layer of chlorine gas of 1 atmosphere, 30 cm 
long. The intensity of the light passing the HO» 
was measured by a photoelectric cell which 
indicates small changes of intensity much better 
than the photographic plate. A cadmium cell, 
argon filled, in a quartz bulb, made in this 
laboratory by Mr. H. W. Leighton, was used. 
The visible light transmitted by the chlorine 
filter has no effect on the cadmium cell since its 
long wave-length limit is \3100. The current was 
amplified with a type 32 tube and measured 
with a microammeter. This system—lamp, 
chlorine filter, cadmium cell—was calibrated for 
its power to indicate the presence of H2O2 by 
H.O2. vapor of known concentration. The 


8H. C. Urey, L. H. Dawsey and F. O. Rice, J. Am. 
Chem. Soc. 51, 1371 (1929). 

‘* Concentrated HO. was made,for other experiments 
which will be described in another paper. 


647 


calibration showed that in the same absorption 
tube of 150 cm length HO: of 1/100 mm could be 
detected. 

We are interested as to whether or not in the 
electric discharge through water vapor the OH 
radicals when recombining form H2O:s. The 
absorption experiment gave no trace of HQ». 
This result indicates that the considerable 
amount of H2Os, frozen out in the experiment of 
Rodebush and Wahl in the liquid-air trap, is 
formed in the trap, not in a gas reaction. 

The formation of HzO: in the trap was men- 
tioned as an alternative explanation by Rodebush 
and Wahl themselves. That this actually takes 
place was concluded by Campbell and Rodebush"® 
from recent experiments. Their result has recently 
been confirmed by Geib,’ who found hydrogen 
peroxide in the gases pumped out of the discharge 
tube only under conditions which according to 
his own experiments are favorable for the forma- 
tion in the trap of H2QOx>. 

In the interpretation of the experiments by 
Rodebush and Wahl the difficulty remains that 
the concentration of HzO. observed in the trap 
does not appreciably change, while the time of 
flow between discharge bulb and trap (velocity 
of flow and distance) is varied. It must be 
concluded that no appreciable combination of 
OH +H takes place during the time of flow. This 
might be largely due to the fact that they worked 
with pressures as low as 0.1 or 0.2 in the discharge 
tube, hence still somewhat lower in the flowing 
system, so that triple collisions consuming the 
OH radicals are still less probable than in the 
experiments reported above. Furthermore the 50 
percent yield of HO: reported in their experi- 
ments means the value relative to HO, not to the 
total gas content pumped out of the discharge. 
At any rate, their observation confirms the 
conclusion of section 6 that no fast combination 
of OH on a clean glass wall takes place since even 
the gas slowly flowing through long glass tubes 
shows the same yield of H,O, in the trap. 


VIII. PROCESSES IN THE DISCHARGE 
THROUGH H.O 


It remains to summarize the chemical processes 
taking place after interrupting the discharge 
16 R. W. Campbell and W. H. Rodebush, J. Chem. Phys. 


4, 293 (1936). 
16K. H. Geib, J. Chem. Phys. 4, 391 (1936). 
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through water vapor or else in the gases pumped 
out of the discharge. We are not concerned with 
the discharge proper in which reactions may be 
caused by ions as discussed in detail by Brewer.!” 

Again we assume H and OH as primary 
products, since O is produced at a comparatively 
small concentration. The main result of the 
present paper is that the disappearance of OH 
is largely due to a termolecular gas reaction. 
From our failure to observe the absorption 
spectrum of HO. and from the results of 


17 A. K. Brewer, J. Phys. Chem. 38, 1051 (1934). 


HENDRICKS AND MAXWELL 


Campbell and Rodebush and of Geib, H+OH 
+M-—H,0+M is to be assumed as the most 
plausible process. A wall reaction is certain not 
to play a predominant part for a clean glass 
surface but to become important at the KCl 
surface. 

In addition it is possible that HO» is produced 
but rapidly decomposed so that no appreciable 
concentration builds up, since H atoms and OH 
radicals are both able to decompose it.!® 


18 Cf. K. H. Geib, Zeits. f. physik. Chemie A169, 161 


(1934). Evidence for the decomposition process OH+H.20» 
will be presented in another paper. 
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Electron diffraction photographs (A=0.062A) were 
obtained from about 10-7 g of polonium that had been 
volatilized in a stream of hydrogen and condensed over 
an area of about 3 mm? on a thin collodion film. Diffraction 
patterns were also obtained from bismuth and tellurium 
since it was expected that polonium would have a similar 
crystal structure. Analysis of these patterns shows that 
the structure of polonium closely resembles that of tel- 
lurium, the lattice being pseudohexagonal with a=4.25A, 


INTRODUCTION 


OLONIUM is the superior homolog of 
tellurium in the periodic system of the 
elements and immediately follows bismuth in 
atomic number. It resembles both these elements 
in its chemical properties but, of course, differs 
from them in being radioactive. The crystal 
structure of the element could be expected to 
be similar to that of tellurium or bismuth. 
Sufficiently large amounts of polonium have 
not yet been available for x-ray diffraction 
experiments and thus no information is available 
on its crystal structure or its density. A strong 
source of polonium for radioactive work contains 
*On leave of absence from the University of Milan. I 
am particularly indebted to Professor Dr. G. Bruni, Dr. 
Charles L. Parsons, Dr. H. G. Knight and Dr. C. H. 
Kunsman for making arrangements necessary for the 


successful culmination of this work and for their many 
courtesies. 


c=7.06A, or 14.12A, and the calculated density 9.39 
assuming 3 Po in the pseudo unit of structure. The true 
lattice is probably monoclinic with a=7.42A, b=4.29A, 
c=14.10A and 6 quite close to 90°, a suggested value 
being 6=92°; the calculated density for 12 Po in the unit 
of structure is 9.24. A structure, based upon the space 
group C;3—C2, in which each polonium atom has four 
nearest neighbors gives moderate agreement between 
observed and calculated intensities of reflection. 


ca. 10~* g of material, which is quite sufficient to 
give electron diffraction photographs provided 
that it can be obtained as a thin film. 


PREPARATION OF THE POLONIUM FILM 


Pure polonium was prepared from old radon 
bulbs! under the supervision and with the 
constant help of Dr. L. R. Hafstad? who has 
developed a technique for the production of 
strong polonium sources after the method used 
by Mme. I. Curie.? According to this method 


1 We are indebted to Dr. C. F. Burnham and Dr. F. 
West of the Kelly Hospital, Baltimore, Maryland, for the 
liberal supply of old radon bulbs that they placed at our 
disposal. The purification of the polonium was carried out 
at the National Bureau of Standards. Successful completion 
of the work was due in a large part to the cooperation of 
Dr. Hafstad, of the Department of Terrestrial Magnetism, 
Carnegie Institution of Washington. 

2L. R. Hafstad, J. Frank. Inst. 221, 191 (1936). 

$1. Curie, J. chim. phys. 22, 471 (1925). 
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Bismuth 
A =0.062A 
plate dist. =24.5 cm. 


Polonium 
\ =0.0622A 
plate dist. 


POLONIUM 








Tellurium 
\ =0.059A 


24.5 cm. plate dist. =24.3 cm. 


Fic. 1. Electron diffraction photographs of bismuth, polonium and tellurium (full size). 


the polonium is separated from the solution 
containing salts of radium D and E by displace- 
ment on silver, complete separation is shown by 
absence of 6-ray activity. The polonium was 
finally displaced from solution upon a cleaned 
and outgassed nickel disk. A film was obtained 
by evaporating the polonium from the nickel 
base in a stream of pure hydrogen and con- 
densing the vapor on a very thin collodion film 
supported by a fine platinum gauze. The 
volatilizing apparatus used was similar to that 
described by Rona and Schmidt,‘ and Curie 
and Joliot.* During the evaporation process the 
nickel was maintained at a dull red heat which 
was sufficient to volatilize the polonium® without 
carrying over traces of nickel. A dark visible 
coating was obtained over an area of about 
3 mm? on the collodion film. 

Electron diffraction photographs were taken 
immediately after preparation of the sample 
since the radioactivity of the polonium soon 
destroyed the film. a-ray air ionization measure- 
ments taken after the film was destroyed showed 
the presence of about 1500 e.s.u. or 2.521077 ¢ 
of polonium remaining on the wire support. It 
is estimated that the original thickness of the 
polonium film was not less than 100A. 


ELECTRON DIFFRACTION APPARATUS 


The apparatus used for electron diffraction 
has been described elsewhere.? Electrons with a 
de Broglie wave-length of about 0.06A were 


«. E. Rona and E. A. W. Schmidt, Zeits. f. Physik 48, 784 
28). 

I. Curie and F. Joliot, J. chim. phys. 28, 201 (1931). 

‘ P. Bonét-Maury, Ann. de physique 11, 253 (1929). 
_'S. B. Hendricks, L. R. Maxwell, V. L. Mosley and M. 
E. Jefferson, J. Chem. Phys. 1, 549 (1933). 


used with a film to plate distance of 245.6 mm. 
Preliminary tests were made with collodion 
films upon which lead, bismuth and _ tellurium 
had been volatilized in the apparatus mentioned 
above. The lead films were prepared from about 
10-° g of the element that had been deposited 
upon nickel. In all these cases the diffraction 
photographs obtained were in complete agree- 
ment with the crystal structure data for these 
elements. There was only slight evidence of 
scattering from collodion on any of the photo- 
graphs. 

Five electron diffraction photographs were 
obtained from polonium films at de Broglie 
wave-lengths between 0.0623 and 0.0627A. One 
of these is reproduced in Fig. 1 together with 
photographs from bismuth and tellurium. Aver- 
age values of Q=1 d?=4sin® @ as measured 
on the polonium photographs are listed in 
Table I together with the estimated intensities 
of reflection. 


Tue LATricE AND CRYSTAL STRUCTURE 
OF POLONIUM 


Diffraction rings of polonium can roughly be 
accounted for, as shown in Table I, on the basis 
of a hexagonal lattice having a=4.254A and 
c=7.06A. These lattice dimensions are the 
expected ones if the structure of polonium is 
similar to that of tellurium or bismuth. More- 
over, as shown in Table II, the intensities are 
not greatly unlike those of tellurium. However, 
closer inspection shows that neither the above 
lattice nor the one having c=14.12A can be 
correct. It perhaps can be seen on the electron 
diffraction photograph of polonium reproduced 
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TABLE I. Electron diffraction data from polonium interpreted 
according to a hexagonal lattice. 














Q Pseudo 
Ring hexagonal Intensity* 
No. Obs. Calc indices (obs.) 
1 0.0890 to 
0.0940 0.0919 (10.1) vS 
2 0.1580 0.1520 (10.2) m 
3 0.1804 0.1805 (00.3) ms 
4 0.2175 0.2210 (11.0) w 
5 0.2481 0.2415 (11.1) mw 
6 0.2642 0.2522 (10.3) w 
7 0.3634 0.3730 (20.2) w 
8 0.370 vw 
9 0.4068 0.3925 (10.4) m 
0.4016 (11.3) 
10 0.4480 m 
11 0.46 0.4732 (20.3) vw 
12 0.5310 0.5340 (21.1) s 
0.5418 (11.4) 
13 0.7170 0.7216 (00.6) w 
0.7221 (11.5) 
14 0.8000 0.7940 (20.5) vw 
15 0.8465 0.8436 (30.3) w 
16 0.8953 0.8840 (22.0) w 
0.9040 (22.1) 
17 0.9920 0.9640 (22.2) mw 
1.014 (20.6) 
18 1,149 1.136 (31.3) vw 
1.177 (40.0) 
19 1.248 1.257 (40.2) vw 
20 1.473 1.456 (31.5) w 
1.478 (32.2) 
1.498 (40.4) 
21 1.627 1.624 (00.9) vw 
1.627 (41.2) 








*The following abbreviations are used throughout this work; vs, 
very strong; s, strong; ms, medium strong; m, medium; mw, medium 
weak; w, weak; vw, very weak. 


in Fig. 1 that the first ring is quite broad, as is 
also true for some of the other rings. The (8) 
and (10)—(11) rings cannot be accounted for by 
the hexagonal lattice and the agreement between 
observed and calculated values of Q=4 sin? 0/d” 
is extremely poor for the (2), (7) and (10)—(11) 
rings. This lack of agreement is not changed by 
doubling the value of c. 

Before continuing the direct structure argu- 
ment for polonium it is best to point out some 
of the structural characteristics of the elements 
in the 5th and 6th subgroups of the periodic 
system; the pertinent facts are summarized in 
Table III. Both the bismuth and tellurium type 
structures can be considered as distortions of 
the mercury structure which is itself a deforma- 
tion of a close-packed structure. Each mercury 
atom has six nearest neighbors; in bismuth the 
displacement is 0.15A along the c axis, leading 
to three nearest neighbors, and in tellurium 
0.25A perpendicular to that axis such as to give 
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two nearest neighbors. As the atomic number 
increases in either series the differences in 
distance between nearest and next nearest 
neighbors decreases. Now if the pseudo lattice of 
polonium is approximately correct and if the 
pseudo parameter values are those of tellurium 
then the distances to surrounding atoms are 
those shown in Table ITI. 

The true lattice of polonium is probably a 
monoclinic one and it is logical to suppose that 
the space group is the monoclinic subgroup 
common to both D;*—C3,2 and Ds,a®—R3m, 
namely C,*—C2. A strictly analogous situation 
is found for plaster of Paris, 2 CaSO,-H.O, which 
as studied by x-ray powder diffraction methods‘ 
was found to have a structure based on the 
enantiomorphic space groups D;*—C3,2 and 

3°— C322. Earlier and more thorough work on 
single crystals of 2 CaSO,-H2O,* however, shows 
that the true lattice is monoclinic with 6 =90° 38’ 

TABLE II. Diffraction results from tellurium. 





Q Intensity* 











Obs. Cale,” Indices X-ray Electron Polonium 
0.097 0.0966 (10.1) Ss Ss VS 
.184 .1840 (10.2) ms ms m 
.201 .2024 (11.0) m m w 
2ot .2315 (11.1) w mw mw 
.263 .2620 (00.3) w w ms 
.302 .2990 (20.1) m m 
o22 .3190 (11.2) vw vw 
3295 (10.3) w 
.396 .3860 (20.2) m w w 
409 4645 (11.3) m m m 
4722 (21.0) 
496 .5013 (21.1) w w s 
aoe Ooae (10.4) w w m 
(20.3) 


.5318 








* The formula used was Q =0.06746(h?2 +k? +hk) +0.02911/2 corre- 
sponding to a lattice having a=4.445 and c=5.86A, the value_of ¢ 
being slightly smaller than 5.91A as found by Bradley (Phil. Mag. 48, 
477(1924)) from x-ray diffraction. 


TABLE III. Structures of some elements of the 5th and 6th 
subgroups of the periodic table. 








Hexagonal | 


























| Lattice Neighboring Atoms 
| Dimensions A | 
| 
Space | | | | 
Element Group a c | c/a | No. | Distance | No. | Distance 
| | » 
As | Dag —R3m |3.76|10.57|2.81| 3 | 2.51 | 3 | 3.15 
Sb | D3q® —R3m_ |4.26|11.29| 2.65) 3 | 2.87 3 3.37 
Bi | Dag’ —R3m |4.56/11.85| 2.60) 3 | 3.10 3 3.47 
Se | D3t—C312_— | 4.33 9.94) 2.28) a | ae 4 | 3.46 
Te | Dst—C312 |4.44]11.82] 2.66] 2 2.86 4 | 3.74 
Pseudo polonium | | | | 
lattice 4.25 14.12} 3.32 —| 310 | — | 3.58 











® P. Gallitelli, Periodico di Mineralogia 4, 132 (1933). 
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CRYSTAL STRUCTURE OF POLONIUM 


TABLE IV. Electron diffraction data from polonium with suggested monoclinic indices. 








Pseudo 
hexagonal 
indices 


Observed 
Q 


Intensity 


Calculated 
Q Monoclinic values of 
indices J-S/Fp.? 





0.0890 vs (10.1) 
—0.0940 
0.1580 (10.2) 


(00.3) 
(11.0) 


(11.1) 


(10.3) 


(20.1) 


(20.2) 


(11.3) 


(10.4) 


(20.3) 


(21.1) 


{ (93 ; 230 
(112) 
)} (112) 
| (202) 
( (204) 
) (114) 
| (114) 
| (204) 
| (006) 
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and a/b=1.744 instead of 1.732 and that the 
space group is probably C;*— C2. This conclusion 
has recently been challenged'® but the properties 
of the crystal, particularly the biaxial character, 
gives it strong support. 

The quadratic form for a monoclinic lattice is: 
4 sin? 6 h? k? P 








2 cos B 
Al. 


? asin? B & csin?B acsin? Bs 


Q= 


Now if the lattice is pseudohexagonal about the 
¢ axis B must be near 90° and thus sin 8 must be 
approximately unity while cos 8 is very sensitive 

OW. A. Caspari, Proc. Roy. Soc. A155, 41 (1936). 


Note that the Laue photograph as reproduced does not 
have a threefold axis. 





to deviations from 90°. The pseudohexagonal 
(00.3) becomes the monoclinic (006) and thus 
csin 6B is 14.10A and c is very close to 14.10A. 
If the structure is pseudohexagonal the intensi- 
ties of the reflections from the monoclinic (hk0) 
should differ chiefly by the mvItiplicity factor. 
Thus both (020)(310)—(11.0) should have an 
appreciable intensity. If the limits of Q for the 
fourth line are taken as 0.2175+0.035 then 
b=4.29+0.04A and a=7.42+0.07, giving a/b 
= 1.763 to 1.700 and 


Q=0.0181h?+0.0544k? +-0.0050/? 
—0.0190 cos Bhl. 


In estimating the value of @ it is noted that (1) 
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Fic. 2. A projection showing the relationship of the 
polonium structure to that of bismuth and tellurium as 
obtained by deformation from the mercury structure. 
The full arrows that are parallel to the ab plane represent 
deformation towards tellurium and the curved half arrows 
toward polonium. Bismuth positions are obtained by 
displacements along the normal to the ab plane. 


is broad, that (2) does not agree well with the 
pseudohexagonal spacing, and that the agreement 
between observed and the calculated value of 
(10) for the pseudo lattice is poor (Table I). 
A suggested value of cos8 is —0.035 which 
corresponds to B=92°. Values of Q calculated 
on the basis of this lattice for planes possibly 
contributing to the first twelve lines of the 
polonium diffraction pattern are listed in Table 
IV. The agreement between observed and calcu- 
lated values of Q is quite satisfactory. 

If the space group is C,*—C2 the polonium 
atoms are probably in three sets of general 
positions : xyz; yZ;4-+x, $+y, 2;3—x, $+y, 2. 
Parameters corresponding to the hexagonal 
tellurium positions are 


x: =0.00, x2= 0.14, x3=0.36, 
n= 4 N= —.14, ¥3= me, 
21> so 22> 42, 23> .08. 


The intensity of the pseudohexagonal (10.3) 
which is present on the polonium photograph 
but absent for tellurium cannot be explained by 


11 Int, Tables for the Determination of Crystal Structures, 
(Berlin, 1935), Vol. I. 
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use of these values, nor by considerable distortion 
of them with the symmetry of C3,2 tending to 
equalize the distances of nearest and next 
nearest neighbors. 

Inspection of Table I further shows that the 
pseudohexagonal planes corresponding to the 
(13)—(21) lines either are 3n orders of (00./) or 
usually have low values of /. The probable 
explanation for the latter fact is that the spacings 
of the monoclinic planes derived from a par- 
ticular pseudohexagonal plane are grouped quite 
close together for low values of J. This is well 
illustrated by the strong (12) ring as shown in 
Table IV, there being eight possible contributing 
monoclinic forms. The high intensity of the (3) 
ring and the possible presence of reflections from 
(0012) and (0018) strongly suggests that the c 
parameters of polonium are quite close to those 
of tellurium. 

After considerable work it was found that the 
following parameter values gave moderate agree- 
ment between observed and calculated intensities 
of reflection as perhaps can be seen by inspection 
of Table IV, column 8. 


x1 =0.05, x2= 0.13, x3=0.35, 
n= .27, Yo = — .20, y2= AO, 
24> .245, 2 >= A3, 23> .07. 


It is not considered that these particular param- 
eter values are the correct ones but a deformation 
of this type from the tellurium structure is 
probably required by the data. They lead to 
the following approximate separations of the 
polonium atoms: 


X1VizZ1 tO XeVe%, 3.50A, 3.34A and 4.06A; 
X1Vi21 to X3V323, 3.30A, 3.38A and 3.95A. 


This distortion is in such a direction as to give 
four nearest neighbors at a distance of about 
3.40A about a particular polonium atom. 

A projection of the above structure on (001) 
is shown in Fig. 2 together with the indicated 
relationship to a similar projection of the bismuth 
and tellurium positions. 

The density calculated, for the atomic weight 
210, on the basis of the monoclinic lattice is 9.24 
and for the pseudohexagonal lattice 9.39. The 
calculated atomic volume is 22.7 cc per mole, 
a value to be expected from a consideration of 
the atomic volume-atomic number relationships 
of the elements of high atomic weight. 
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Attempts to produce a photovoltaic effect with reversible electrodes made of an active 
metal photoelectrically sensitive to visible light, such as sodium, and dipping into a solution 
containing ions of the metal, met with no success provided the solvent and salt which formed 
the electrolyte were carefully dried before preparing the cells. Slight photovoltaic effects of 
the order of 0.5 mv or smaller were observed in the case of two cells known to contain traces 
of water, but these effects were probably due either to the water present or to the compound 
electrode Na/Na2O rather than to the pure metal. It is concluded that the energy supplied to 
the electrons by the visible light is not sufficient to bring about the transition of the electron 


from the metal to the solution. 





CONVENTIONS AND DEFINITIONS 


UGHES and DuBridge' use the term 

“photovoltaic effect’? to include ‘‘phe- 
nomena associated with the production of an 
electromotive force in a cell consisting of two 
similar electrodes, separated by a_ suitable 
electrolyte or other substance, when the elec- 
trodes are unsymmetrically illuminated.”’ We 
shall define the photovoltaic potential as being 
positive when the potential of the illuminated 
electrode becomes more positive on illumination 
irrespective of the sign of the electrode potential 
in the dark. We believe that the use of the terms 
anode and cathode in the case of cells generating 
an e.m.f. is inadvisable and confusing if not 
entirely incorrect. The values and sign of the 
e.m.f. of the cells as given in this paper will 
always be the potential of the electrode which 
is illuminated or about to be illuminated in 
reference to the potential of the other electrode 
kept continually in the dark. 

By external photoelectric effect we mean, as 
usual, the ejection of electrons from a metal to a 
vacuum by light. By internal photoelectric effect 
we signify the ejection of electrons by the light 
from one solid phase to another solid phase, the 
ejection of electrons from a layer of cuprous 
oxide to the underlying copper, for example. 


INTRODUCTION 


The interesting question regarding the extent 
to which light may affect the potential of a 
‘A. L. Hughes and L. A. DuBridge, Photoelectric 


senate (McGraw-Hill Book Co., New. York, 1932), 
p. 352. 


metal or metal coated with an insoluble salt in 
contact with an electrolyte has been investigated 
by a number of workers from the original 
discovery of Becquerel’® in 1839 to the present 
day. Hughes and DuBridge* summarize papers 
published through 1929 while Audubert* and 
Athanasiu’ have recently published extensive 
reviews. Photovoltaic potentials have been ob- 
served for a number of compound electrodes such 
as Ag|AgI, Cu| CuO, etc., dipping into aqueous 
solutions, but these photopotentials largely dis- 
appear when the electrodes are carefully dried 
and then dipped into electrically conducting 
nonaqueous solutions, the photopotential of a 
cuprous oxide electrode, for example, drops from 
120 mv to 3 mv on elimination of water from 
the system. Audubert explains the photovoltaic 
effects in the case of aqueous solutions by 
postulating that the light ‘‘photolyzes’’ water 
molecules which are adsorbed on the surface of 
the electrodes to form either hydrogen or oxygen 
which then enter into equilibria involving ions of 
the metal, thereby affecting the potential of the 
electrode. Consider, for example, the Cu|Cu,O 
electrode whose potential becomes more positive 
on illumination ;° for this electrode the reaction 


Cut++10.+3H,0@Cut+++OH- 


is applicable. The incident light photolyzes the 
water, increasing the concentration of the oxygen 


and forcirig the above equilibrium to the right. 


2 E. Becquerel, Comptes rendus 9, 144, 561 (1839). 

’ Reference 1, Chap. IX. 

4R. Audubert, J. de phys. et rad. (7) 5, 486 (1934). 

5G, Athanasiu, Ann. de physique (11) 4, 377 (1935). 

6 At least in the case of cuprous oxide made in the dry 
way. 
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The resulting increase in the concentration of 
the cupric ion produces the positive photo- 
potential. In support of this hypothesis may be 
cited the fact that in the case of the Cu| CuO 
electrode its positive photopotential is increased 
in reducing solutions and reduced in oxidizing 
solutions, is reduced by increased concentration 
of cupric ions and is reduced by increase of pH. 

Athanasiu® has studied the influence of pH, 
temperature and wave-length on the photo- 
voltaic effect. His results may be summarized 
as follows; decrease of PH favors the positive 
photoeffect, increase of temperature lowers the 
photopotential markedly, increase of wave- 
length of the incident radiation decreases the 
positive photopotential and may produce a 
negative photopotential. In the case of Ag| Agl 
electrode immersed in a neutral solution of KCl 
after a basic treatment the positive photo- 
potential shows a maximum at 4200A and the 
negative photopotential a maximum at 4400A. 
The maximum of the spectral absorption coeffi- 
cient of silver iodide, the maximum of its 
photoconductivity and the maximum of the 
photovoltaic effect all occur in the same wave- 
length range, about 4200A, but the threshold of 
the external photoelectric effect is below 4070A. 
Athanasiu concludes that the photovoltaic effect 
begins with an internal photoelectric effect of the 
light. sensitive salt or oxide layer, followed by 
photochemical reactions in which the electrolyte 
may take part. In addition to a photochemical 
reaction there must also be a possibly reverse 
thermal reaction in order to account for the 
pronounced temperature coefficient of the photo- 
voltaic effect. 

The photopotential of pure metals has been 
investigated in the case of noble metals which 
may be assumed to be freer of oxide coating 
than more active metals like copper or zinc, 
although the presence of adsorbed layers of gases 
makes it impossible to assert that a pure metal 
surface has ever been studied for its photovoltaic 
effect. The photopotentials of the noble metals 
are much smaller than those discussed above and 
rather difficult to measure because of the 
instability of the electrode potentials in the dark ; 
the electrodes probably acting as oxygen elec- 
trodes. 
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Audubert and Roulleau’ studied the photo- 
voltaic effect for selenium on platinum and found 
large photopotentials even in the absence of 
water. They used solutions of potassium iodide 
dissolved in a number of organic solvents, and 
found in contrast to the Cu| CueO electrode that 
the observed positive photopotentials were in- 
creased in oxidizing solutions. Very few experi- 
mental details were given, so that the purity of 
the selenium surface is open to question; the 
authors believe that a true photoelectric effect 
was concerned and that electrons were ejected 
into the solution by the action of the light. 

In all previous work on the photovoltaic effect 
no investigations have been made to our know!- 
edge on metals having low work functions like 
sodium. It would seem that if a true photoelectric 
effect existed, for the boundary metal solution, 
it would be much more pronounced for active 
metals than for the relatively inert metals 
having high work functions. Accordingly we 
decided to investigate the photovoltaic effect 
for pure sodium surfaces in contact with a 
solution containing sodium ions. We hoped in 
this way to be dealing with electrodes giving 
reversible and stable potentials in the dark and 
with metal surfaces as clean as they could be 
when in contact with a liquid electrolyte. Sodium 
electrodes also have the advantage that there is 
no photosensitive oxide or insoluble salt layer 
which would make possible an internal photo- 
electric effect or photochemical reactions with 
the electrolyte. 


EXPERIMENTAL 


The general method of investigating the 
photovoltaic effect for sodium consisted in 
preparing a dry nonaqueous solvent, dissolving 
dry sodium iodide in it, pouring the solution 
onto sodium electrodes and finally measuring 
the e.m.f. between the two electrodes before and 
after the illumination of one of the electrodes. 
The cell used was Na| Nal in C;H;,NH2| Na. 

Sharples Solvent Company technical grade 
monoamylamine was selected as the solvent 
because it is conveniently a liquid at room 
temperature. It was dried by refluxing over 
liquid sodium-potassium alloy for three hours 


7R. Audubert and J. Roulleau, Comptes rendus 198, 
1489 (1934). 
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followed by repeated distillations in vacuum into 
flasks containing sodium mirrors until the solvent 
was so dry that a thin sodium mirror would no 
longer be removed. Connections between flasks 
were opened by breaking an inner bulb with a 
magnetic hammer; a large bulb was broken 
rather than a capillary tip in order to expedite 
the distillation of the solvent and later on the 
pouring of the solution from one flask to another. 
The iron magnetic hammer was encased in glass 
in order that iron rust would not catalyze the 
reaction between sodium and the amine to form 
hydrogen and an amide.* Various methods were 
tried in an attempt to dry sodium iodide re- 
crystallized from water by fusion, but invariably 
a saturated aqueous solution of the resulting 
fused iodide was alkaline to phenolphthalein so 
that the drying was usually done by heating the 
salt at a pressure of 1 micron from room temper- 
ature to 100°C and then at 300°C at a pressure 
of 0.1 micron for ten hours. This process gave a 
neutral salt which on solution in dry amyl amine 
contained insufficient water to remove a thin 
sodium mirror. The solutions were made by 
distilling the solvent onto the salt, followed by 
shaking. The concentrations of the solutions were 
estimated roughly by evaporating the solvent 
from a weighed quantity of the solution and 
weighing the dry residue or by adding a known 
volume of liquid to a known weight of salt. 
A saturated solution of the salt contains approxi- 
mately 27 percent of sodium iodide by weight. 
The sodium for the electrodes and mirrors was 
melted in vacuum, filtered through a capillary 
tube to remove oxide and finally pushed by 
means of hydrogen or allowed to run by gravity 
through a second capillary tube into a flask or 


’T. H. Vaughn, R. R. Vogt and J. A. Nieuwland, J. Am. 
Chem. Soc. 56, 2120 (1934), have discovered that ferric 
nitrate dissolved in liquid ammonia catalyzes the amide- 
forming reaction between sodium and ammonia. 
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cell. The side arms containing the sodium 
reservoir and capillary tubes were removed after 
the sodium mirror had been formed. The cells, 
described below, were at first maintained at 
constant temperature, 25°, in an air bath; but 
later it was found that the actual temperature 
of the experiment was of no apparent importance 
so that the air in the bath was merely stirred 
vigorously, the essential precaution being that 
the two sodium electrodes be at the same 
temperature. An e.m.f. of 0.5 mv or higher 
could be produced by touching one side of a cell 
with the hand. The e.m.f. was measured usually 
to 0.01 mv with a type K-2 Leeds and Northrup 
potentiometer using a Leeds and Northrup type 
R galvanometer (sensitivity 510-'° amp. per 
mm scale deflection). 

A 200-watt tungsten filament lamp with a liter 
flask filled with water or saturated copper sulfate 
solution as a condensing lens and heat filter was 
used in the first experiments, but later a mercury 
arc without any filter giving more intense light 
was constructed and used. The arc was enclosed 
in Pyrex glass so that light of wave-length below 
3000A was cut off. The monoamylamine also 
begins to absorb strongly in the ultraviolet 
region below 3000A.° 

The cells illustrated in Figs. 1 and 2 were of 


two types. Type A, Fig. 1, contained one sodium 


( 9 See J. Bielecki and V. Henri, Comptes rendus 156, 1860 
1913). 
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electrode which could be illuminated while type 
B, Fig. 2, contained three sodium electrodes 
capable of illumination. In Fig. 1 at @ and c are 
two tungsten electrodes sealed through the Pyrex 
glass, both of which were covered with sodium 
metal which was admitted through a capillary 
not shown in the figure. A sodium mirror was 
deposited on the inside walls of the cell around 
the electrode c by gently heating the sodium at c. 
Light entered through a thin Pyrex window of 
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solution, preventing measurements to be made 
to 0.01 mv. Water was purposely added to cell 
A-8 causing gas evolution at the electrodes and 
producing an erratic e.m.f. but after the water 
present had combined with the sodium leaving a 
thick layer of hydroxide on the electrode surface 
an e.m.f. could be measured, however, no 
photopotential was observed. No photopotential 
was observed in the case of pure solvent in 
contact with the sodium electrode, cell No. A-6, 































































the type described by Sonkin,'® at d while but here the measurements were not very ( 
solution was poured into the cell at 6. After the accurate due to the high internal resistance of the \ 
filling of the cell the side arm at b was sealed off. cell. Cell No. A-3 was deliberately opened to I 
In Fig. 2 is illustrated a cell of type B and the _ the air after preparation, a small photopotential V 
auxiliary capillary tube for admitting the sodium was indicated by a slight motion of the galva- f 
is also pictured. In some cells of this type no nometer on illumination but it was not definite ! 
sodium mirror was formed, the sodium electrodes enough to measure. In the case of cell No. A-5 i 
simply consisting of a pool of sodium metal at air was purposely admitted to the cell after the 
the bottom of each arm of the cell. sodium mirror had been formed and before the c 
The results of the investigation are collected solution was poured in, thereby making it e 
in Table I which is self-explanatory but to which probable that a thin layer of sodium oxide r 
the following remarks should be appended. In covered the mirror. No photopotential was n 
all the cells except No. B-1 the electrodes were observed, possibly because of the insensitivity of e 
mirror electrodes, the electrodes of cell B-1 were the measurements or because of the thickness of SI 
formed of sodium in the bulk. After a zero the oxide layer. ‘b 
photopotential had been observed for all the The only cell which showed a definite photo- tl 
electrodes of cells B-1 and B-2 sodium was voltaic effect was cell No. A-2, but this cell li 
deposited electrolytically in the form of shiny contained some water since on standing the d 
white needles on one electrode of each of the sodium covering one of the tungsten wires eC 
cells, but no photopotential resulted. The disappeared bringing the tungsten into direct sl 
electrical resistance of cells B-1, A-1, A-5, A-6, contact with the solution and causing the dark cC 
and A-7 was rather high due to the dimension of e.m.f. to rise to 0.9 v. Before this happened, the gc 
the cells or to the low concentration of the e.m.f. of the cell in the dark changed gradually m 
TABLE I. 
Dark e.m.f. of dennde 
No. of to be illuminated in Photo 
Cell Concn. of Soln. Light illuminated reference to other e.m.f. ok 
oO. in wt. percent Source electrodes electrode, millivolts millivolts it 
A-1 3.9 lamp 1 8 0 _ im 
A-2 18.2 lamp 1 25.10—16.5 +0.5+0.05 
B-1 27 arc 3 B-A 76.1 A 0+0.1 - 
C-A 16.0 BO 
B-C 61.3 C0 
B-2 27 arc 2 B-A 0.20 A 0+0.01 
BO 
A-3 27 arc, lamp 1 —0.45 trace +0.01 
A-4 27 arc, lamp 1 —0.18 0+0.01 
A-5 5 arc, lamp 1 —5.5 0+0.2 
A-6 0 arc, lamp 1 +0.1 v 0+0.1 v 
A-7 arc, lamp 1 26 0+1 
A-8 27 arc, lamp 1 


















10S. Sonkin, J. Opt. Soc. Am. 19, 65 (1929). 
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from +0.025 v to —0.0165 v while the photo- 
potential amounted to 0.5 mv irrespective of the 
dark e.m.f., the illuminated electrode becoming 
more positive. This photovoltaic production of 
an e.m.f. occurred instantaneously on illumina- 
tion as far as could be told, but it was probably 
due to a layer of oxide or hydroxide on the 
surface of the sodium since no photovoltaic 
effect was ever observed for a perfectly dry cell 
in which the sodium surface was apparently 
clean. Cu|CusO electrodes in the presence of 
water exhibit a large photovoltaic effect so it is 
reasonable to suppose that a Na| Na2O electrode 
would act similarly. Sodium iodide dried by 
fusion was used in the case of cells A-1, A-2 and 
A-7 and sodium iodide dried in a vacuum at 300° 
in the case of the remaining cells. 

In summarizing the experimental results we 
can say that no measurable photovoltaic effect 
exists for a clean surface of sodium in equilibrium 
with an anhydrous solution of sodium iodide in 
monoamylamine. Sodium surfaces formed by 
electrolysis, by condensation from the vapor 
state and by solidification of the liquid metal 
‘behaved alike. Changing the concentration of 
the solution and the source of the illumination 


likewise produced no photovoltaic effect. We 
do not believe, therefore, that a pure metal in 
equilibrium with a solution of its ions and of 
such a nature that no photochemical reaction 
can occur will show a photovoltaic effect analo- 
gous to the photoelectric effect exhibited by 
metals in contact with a gas or evacuated phase. 


DISCUSSION 


Having come to the conclusion that no 
observable photovoltaic effect exists for sodium, 
it is next of interest to consider what theoretical 
implications are inherent in this conclusion. In 
order to measure the e.m.f. to 0.01 mv a current 


FOR SODIUM 657 


of the order of magnitude 10-* amp. must flow 
through the galvanometer; the galvanometer, 
however, is sensitive to currents of 510-1!" 
amps. Currents measured usually in photoelectric 
experiments vary from 10-" to 10-'> amp., but 
may become as great as 10-* amp.,'' depending 
on the conditions, intensity of light, activity and 
dimensions of metal surface, presence of gas, etc. 
In our case we can conclude that if a photo- 
voltaic current was produced it must have been 
smaller than 10-° amp. 

The probable reasons for the absence of a 
photovoltaic effect are; first, the absence of any 
chemical reaction which might take place under 
the influence of light such as the photolysis of 
water postulated by Audubert and _ second, 
unfavorable energy relationships. When visible 
light, or light of shorter wave-length, falls on a 
sodium surface exposed to a vacuum, electrons 
are ejected into the vacuum because the energy 
of the light added to the energy of the electron 
in the metal produces electrons having high 
enough energy to escape into the vacuum. In 
our case it may be that the stable energy level, 
if indeed such a level exists, of an electron in 
monoamylamine solution is so high that electrons 
in the metal activated by the energy of the visible 
light are unable to make the transition metal to 
solution. The presence of the positive sodium 
ions in the solution complicates matters, how- 
ever, since the energy relations of the electrons 
may be modified by sodium ion transitions 
between the metal and the solution; for this 
reason quantitative calculations are at present 
impossible, but it is evident that the presence of 
the sodium ions does not change the energy 
relationships sufficiently to produce a photo- 
voltaic effect. 


1 Reference 1, p. 435. 
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Liquid Hydrogen Chloride. Electrical Conductance and Dielectric Constant 


GEORGE GLOCKLER AND R. E. PeEck,! University of Minnesota 
(Received July 24, 1936) 


Liquid hydrogen chloride has a specific electrical conductivity of 0.35 107° reciprocal ohms 
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at —85°C. The temperature coefficient is negative: —0.65X10~'° in the range from —84°C 
to —64°C. The molecular conductivity at infinite dilution is estimated to be 184 (—85°C) 
and 192 (—53°C). The ionization constant is 0.083107'® at —78°C and 0.051X1i07'® at 
— 68°C. The change in heat content of the ionization reaction is: 


(x+y+2)HCI—(H2Cl- HCl,)*+(Cl-HCI,)-, 


AFiigs =— 3500 cal. 








INTRODUCTION 


IQUID hydrogen chloride has been studied 

as an ionizing solvent by various investi- 
gators? and for this purpose its dielectric con- 
stant, electrical conductivity and viscosity must 
be known. These quantities have been deter- 
mined earlier.’ We wish to report here new values 
for the first two properties since they differ 
considerably from the figures given in the 
literature. 

The conductivity cells. Since we wished to de- 
termine the electrical properties of liquid hy- 
drogen chloride in the temperature range from 
—90° to about 15°C it was necessary to construct 
cells which would withstand pressures from 1 to 
50 atmospheres. We constructed five different 
types of cells but found the first four designs 
undesirable for various reasons.' The final design 
which satisfied all our requirements is shown in 
Fig. 1. These all-glass cells withstood the 
pressure involved satisfactorily and it can be 
seen that the platinum electrodes (nonplatinized ) 
are not affected by the pressure changes involved 
when the cells are placed in various tempera- 
ture baths. The tungsten inseals were covered 
with glass along their full length so that the 
platinum wires only were exposed to the hy- 
drogen chloride. The cells were washed with 
dilute hydrochloric acid and conductivity water. 
The cell constants were determined with 0.0125N 
hydrochloric acid and were 0.00373, 0.00321, 


! This article is based upon a thesis presented to the 
faculty of the Graduate School of the University of 
Minnesota by R. E. Peck in partial fulfillment of the 
requirements for the degree of Doctor of Philosophy. 

2P. Walden, Ostwald-Drucker Handbuch der allgem. 
mor IV (1 and 2); Akad. Verlags G. m.b.H. (Leipzig, 

3 Reference 2 and I. C. Tables. 


658 





0.00451 and 0.00400 for the four cells used. 
The cylindrical electrodes were 14.0 cm long; 
the inner cylinder had a radius of 0.60 cm and 
the outer one 0.80 cm. The distance between the 
two electrodes was 0.1 cm. The cells were baked 
out at 415°C for four days and hydrogen chloride 
was distilled into them. They were heated 
between each run, having been flushed with 
liquid hydrogen chloride. 

The cell constants should have been deter- 
mined with a standard liquid of much lower 


specific conductance than the hydrochloric acid 


used in order to avoid the Parker effect.* How- 
ever the conductivities to be measured were 
such low order that this correction appeared to 
be unnecessary. 

Preparation and drying of hydrogen chloride. 
The gas was generated by dropping concentrated 
sulphuric acid onto solid sodium chloride in a 
thoroughly evacuated and flamed out system. 
Two drying agents were employed: Phosphorus 
pentoxide and a combination of calcium chloride 
and magnesium perchlorate. 

While Fairbrother*® has shown that phosphorus 
pentoxide is a satisfactory drying agent for 
hydrogen chloride gas, if the latter is not very 
wet, we find that it is unsatisfactory when the 
gas is to be liquefied and then used for the de- 
termination of the electrical conductivity of the 
liquid. The second drying agent mentioned gave 
a much better product for the purpose of con- 
ductivity measurements. 

It should be noted that the hydrogen chloride 
did not come in contact with mercury surfaces 
after it was dried and fractionally distilled. 
Mercury vapor also was excluded from the 








4 Th. Shedlovsky, J. Am. Chem. Soc. 54, 1411 (1932). 
5 F, Fairbrother, J. Chem. Soc. 2, 1539 (1933). 
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Fic. 1. Conductivity cell for liquefied gases. Pressures up to 
30 atmospheres. 


system by appropriate means. The reason for 
these precautions was the fact that we could 
not obtain consistent results unless mercury was 
excluded. The results of the conductance meas- 
urements are shown in Table I. The measure- 
ments marked ‘“‘fresh’’ were made as soon as 
the cell had been removed from the liquid-air 
bath, where the gas was condensed, and allowed 
to warm up to the temperature of the cryostat 
(—85°C). Agitation and standing would in- 
crease the conductivity several fold. Similar 
observations have been made by others, for 
example in the case of water by Noyes.® 

The dielectric constant of liquid hydrogen 
chloride. The above-mentioned measurements 
were carried out with two different resistance 
bridges: one a slide wire bridge actuated by a 
homemade oscillator (500 and 2200 cycles) and 
a two-stage amplifier and a Jones conductivity 
bridge? No. 260418 at 500, 1000 and 2000 
cycles. With the latter instrument it was possible 
to determine the dielectric constant from the 
reading of the balancing condensers. The results 
are given in Table II. 


*A. A. Noyes, The Electrical Conductivity of Aqueous 
Solutions (Carnegie Institute, Washington, D. C., 1907), 
page 206, 

’P. H. Dike, Rev. Sci. Inst. 2, 379 (1931). 


OF LIQUID HCl 


DISCUSSION OF RESULTS 


The values of the specific conductivity of 
liquid hydrogen chloride reported here are lower 
than any found in the literature. It should be 
noted that the conductivity of freshly prepared 
hydrogen chloride is lower than that of the 
material after standing. In twenty-four hours 
the conductivity (at —85°C) changed from 
0.7 to 5.4xX10-* mho (expt 10) and from 
0.35 to 2.1X10-§ mho (expt 11). It seems that 
no better results can be obtained in glass-vessels 
and employing the usual methods of purification. 

The temperature coefficient of the specific 
conductivity has been found to be negative 
(—0.065 x 10-° mho/degree). 


TABLE I. The specific conductance of liquid hydrogen chloride 
(0.5, 1.0 and 2.0 kilocycles). 





L specific 
conductance 
at —85°C 
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TABLE II. Dielectric constant of liquid hydrogen chloride 
(1.0 kilocycle). 








Molar 
Polarization 
P (cc) 


density 
(grams/cc) 


Dielectric 
Const. 


temp. 
mr 





1.191 
1.179 
1.141 
1.046 
0.993 
0.980 


10.10 
9.80 
8.92 
7.33 
6.55 
6.32 


—85.0 
— 80.5 
— 66.5 
— 34.8 
—19.0 
—15.0 
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Average P 23.5 cc. 
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H20 0.01780 
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90 .157 
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TABLE IV. K;.p. and K;,, for the ionization of liquid hydrogen chloride. 








Temperature: °C 
LX 10° (mho) 

C (moles/liter) 
(A= 1000L/C) x 107 


Noo 

(a=i/d,,) X 109 

(K;i.p. = (H*)(CI-)) x 10% 

(K;, = (H*)(Cl-)/(HC1)) x 10° 


— 83 
3.36 
52.59 
1.032 
184 
0.562 
3.30 
0.103 


— 88 
3.70 
32.90 


Remarks 
Run 11(f) 


— 68 
2.38 
31.42 
0.757 
188 
0.402 
1.59 
0.051 


—78 
3.05 
$2.15 
0.948 
186 
0.508 
2.67 
0.083 


—73 
2.72 
31.80 
0.855 
187 
0.456 
2.03 
0.066 


Table III 








The dielectric constant decreases with in- 
creasing temperature (Table II). If we assume 
that the usual expression for molar polarization® 
holds for liquids then we find the molecular 
polarization for liquid hydrogen chloride to be 
23.5 cc. Smyth and Hitchcock® report 23.1 cc 
at —112°C. 

The equivalent conductance. It is possible to 
calculate the equivalent conductance of liquid 
hydrogen chloride from the specific conductance 
and density. The latter has been obtained by 
Maass and Wright.!® If we could obtain the 
equivalent conductance of hydrogen and chloride 
ions at infinite dilution, we could calculate the 
ionization constant and the ion product constant 
for the ionization of liquid hydrogen chloride. 
It is possible to obtain a very reasonable estimate 


for the equivalent conductance (A,,) of hydrogen 
chloride at infinite dilution in liquid hydrogen 
chloride by the use of a relation established by 
Walden? 


Noo’ Neo X= €-const., 


where 7,,=the absolute viscosity of the solvent 
in poises, x = the association factor obtained from 
the temperature variation of the molar surface 
energy and e=the dielectric constant of the 
solvent (see Tables III and IV). 

The values of the ionization constant at 
various temperatures enable us to calculate the 
temperature coefficient of its logarithm (—0.02) 


8C, P. Smyth, Dielectric Constant and Molecular Struc- 
ture. (Chem. Catalog Co. N. Y., 1931.) 

°C, P. Smyth and C. S. Hitchcock, J. Am. Chem. Soc. 
55, 1830 (1933). 

10Q, Maass and C. H. Wright, J. Am. Chem. Soc. 43, 
1098 (1921). 


and by means of the well-known relation 


(d log K;)/(dT) =AH/2.3RT? 
we find 


AH195= —4.6 X 0.02 X (195)? = —3500 cal. 


This value is the heat absorbed for the ionization 
act of HCI in liquid hydrogen chloride 


(x+y+2)HCI-(H.CI-HClI,)++(Cl-HCl1,)~, 
Al] 195 =— 3500 cal. 


Several weak acids have been shown by Harned 
and Embree"! to have negative heats of ionization 
especially at higher temperatures. These calcula- 
tions are based on the lowest values of the con- 
ductivity of freshly prepared liquid hydrogen 
chloride and since these conductivities are very 
low and hence difficult to measure, it must be 
realized that the ionization constant, the ion 
product constant and the heat of ionization are 
not determined to a high degree of accuracy. 
It is clear that the calculation just given has 
significance only in the case that the measured 
conductance is really due to the conductivity of 
liquid hydrogen chloride and not due to some 
impurity. The three following facts cause us to 
believe that we have measured the conductivity 
of liquid hydrogen chloride: (1) The method of 
purification involved several (three to seven) 
fractional distillations; (2) The fresh material 
used for the measurements was kept at low 
temperatures (—180° to —60°C), where solu- 
bility rates of impurities are slow; (3) The 
conductivity rose on agitation and standing. 


1H. S. Harned and N. D. Embree, J. Am. Chem. Soc. 
55, 1050 (1934). 





OCTOBER, 1936 


JOURNAL OF CHEMICAL 


PHYSICS VOLUME 4 


Absolute Rates of Four-Atom Reactions 
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Equipotential surfaces such as one encounters in valence problems and particularly in the 
calculation of absolute rates, are discussed in a ‘‘bond space”’ instead of configuration space, 
with much simplification resulting both in the general aspects and in the practical handling of 
such problems. The exchange part, as well as the so-called coulombic part, of the bond energies 
is seen to obey a simple vector addition rule. A complete treatment of a four-atom reaction 
without symmetry (H2,+ICI—-HI+HCI) is carried out and results are compared with the 


existing experiments. 





INTRODUCTION 


HE simultaneous interaction of more than 

two atoms with each other is a complicated 
theoretical problem. Insofar as it manifests itself 
in chemical exchange reactions it can be de- 
scribed on the following four assumptions: 


I. When an atom is joined by at least sufficiently many 
bonds to saturate its valence the lower potential energies 
will correspond to configurations which makes additional 
bonds to this atom as small as is possible without weakening 
the main bonds. 

Ila. The energy of interaction of more than two atoms 
can be expressed as a function of single bond additive and 
exchange energies. Here we are using the word ‘‘bond” in 
a somewhat wider sense than usual, meaning electron pair 
energy which two atoms at any particular distance would 
exhibit in the absence of other atoms. 

IIb. The interaction energy a of a pair of atoms consists 
of two parts: one part A arising from ‘‘classical” forces, 
which is purely additive, and a second a@ arising from 
quantum-mechanical electron exchange. This latter part 
is not simply additive and is responsible for the valence 
properties of atoms. 

III. In the approximation considered here, the addition 
law for the latter type of binding is taken to be that 
following from perturbation theory calculations. For the 
case of three and four atoms the following well-known 
formula is obtained:! 


E=A+B+C-—(}[(a—8)?+(@—v)*+(y—a@)?])#, (1) 











a 
Fig. 1. Fic. 2. Vector addition of exchange 
binding. 

'F, London, Zeits. f. Elektrochemie 35, 552 (1929); 
Eyring and Polanyi, Zeits. f. Physik B12, 19 (1931). 


where a2=A+a, b=B+8, c=C+y, are the bonds in the 

case of three atoms, while for four atoms (see Fig. 1), 
a=d\+d, b=)h\+-b2, 
A =A+Az2, etc., 


a=ai+az2, etc. 


cC=Ait+c2 


IV. The actual amounts of additive and of exchange 
energy for each pair are taken to be constant fractions of 
the total interaction energy as represented by its Morse 
curve. 

A,=)\, 
By,=)h, 
Ci= a, 


—]= (1 —r)a, 
Bi=(1—A)hi, (2) 
i= (1 —))c1. 


The assurance with which these four assump- 
tions can be accepted probably decreases in the 
order in which they are quoted. It will naturally 
be very interesting to see to what extent these 
assumptions, taken separately, can be justified 
or improved upon by comparison with experi- 
ment. Even the most specific assumption IV has 
not been shown to be in disagreement with ex- 
periment for configurations near the activated 
state though for interatomic distances less than 
the equilibrium distance of two atoms the frac- 
tion \ certainly decreases with decreasing dis- 
tance. When more experimental information is 
available we may find it desirable to replace the 
assumption of a percentage \ not varying with 
distance, by more general relations. There is 
nothing in our present procedure that would not 
be equally well adapted to such a generalization. 
We shall find that assumptions I and II, to- 
gether with very much less specific assumptions 
than III and IV, permit us to draw conclusions 
based only on the potential curves of the atom 
pairs and on geometrical relations. In view of the 
semi-empirical basis of (1) it is desirable to carry 
the more firmly established assumptions to their 
farthest logical conclusions before making more 
specific assumptions about (1). For example, the 
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conclusions that activated complexes for three 
monovalent atoms should be linear and that 
activated complexes of four monovalent atoms 
are planar or more generally that configurations 
are more stable the smaller the exchange binding 
between the atoms not paired in the ordinary 
chemical sense rest on a more general basis than 
III and IV and yet these conclusions are sufficient 
to enable us to find the “‘best reaction surface’ 
in bond space as we shall see. 

We propose to describe the potential function 
E of a three or four-particle system in bond 
space instead of the usual configuration space. 
That is, we shall represent our system by using 
the three bond energies a, b, c as orthogonal 
coordinates. That this is possible for four atoms 
depends on the fact that the three pairs of bonds 
a; and az», b; and be, c; and cz enter into (1) only 
through their sums. Apparently it is generally 
true that the interaction energy of an even 
numbered particle system, 2”, involves no more 
linearly independent combinations of bonds to 
express the energy than does the one for (2m —1) 
atoms. 

By abandoning the configuration space one 
sacrifices the foremost advantage of that treat- 
ment, namely that the configurational changes 
of the system in a reaction may be represented 
by a mass point subjected to the known potential 
E and obeying the laws of ordinary mechanics. 
But such a representation turns out to be not 
feasible for more than three particles, quite 
apart from the complications in visualizing which 
arise from the many dimensions. 

The success and the simplicity of the new 
treatment depend on the following fact: In 
general, of the six parameters specifying a con- 
figuration of four particles in space, only two, say 
a and b, need be considered as actually inde- 
pendent and as determining a surface B which 
will include the best reaction path. This surface 
will be considered in detail later. A reaction 
may be described by saying that the system 
starts from a configuration where a has its 
maximum value and )=c=0, and proceeds to a 
final state where } is a maximum and a=c=0. 
By the ‘‘best reaction path’ is meant that se- 
quence of configurations by which this transition 
is accomplished in a way most economical of 
potential energy. Thus each of these configura- 
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tions must have a minimum potential energy in 
all directions except the one in which the reac- 
tion proceeds. The point of highest potential 
energy along this path represents the activated 
configuration. It is clearly true that for this 
point—a saddle point—the energy is stationary 
in all directions, a fact of fundamental impor- 
tance in that it allows us to use the theory of 
small vibrations for setting up the partition 
function of the activated state. Mathematically, 
the activated configuration is thus determined 
as the point with the highest potential along the 
best reaction path for which dE=0. 

It is characteristic of the theory of reaction 
rates that the original and the activated con- 
figuration are the only ones whose physical 
properties determine the rate of reaction. Only 
in the case of two saddle points of nearly equal 
heights is it necessary to consider both. 

Since we propose to describe the progress of a 
reaction in bond space we must consider how 
points in bond space and in configuration space 
correspond to each other. In the case of three 
atoms, the three bonds are just sufficient to 
determine the configuration, except for a certain 
ambiguity which is readily disposed of.? In the 
case of four atoms, however, the situation is com- 
plicated by the fact that it takes six parameters 
to determine a configuration so that there exist 
in general »* configurations for a point (a, ), c). 
On the other hand there are regions in bond 
space which are not attainable by a particular 
system owing to the obvious restriction on the 
greatest absolute value a bond can take, and also 
owing to the geometrical restrictions which pre- 
vail among the six distances between four atoms. 
It so happens that in spite of this general situa- 
tion the best reaction path, because of the condi- 
tions on ¢c, always lies on the boundary surface 
between the physically attainable and the un- 
attainable region, and for this surface B there 
does exist a point to point correspondence be- 
tween bond and configuration space. 

We now show that the best reaction path lies 
completely on the boundary surface. We start 


2 This ambiguity arises in the case of four, three and even 
two-atom systems, from the fact that two different 
distances are always associated with every negative value 
of E occurring in the diatomic potential curve. No compli- 
cation arises from this fact since the reaction path never 
involves configurations where an interatomic distance is 
contracted beyond the equilibrium distance for that pair. 
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out by finding for every given pair of values, 
a, b, the smallest possible value of |c¢| which is 
attainable by our system under its geometrical 
restrictions. The surface so defined in bond space 
will clearly be a boundary surface for the attain- 
able region. Its equation can be set down as 


d(Cit¢2)a, »=0. (3) 


On the other hand, if we restrict our attention to 
all those points (a, 0, c) which for given values 
a, b, have their minimum value of £, then the 
surface so obtained must completely contain the 
best reaction path. But this condition can be 
written 


(dE) a, »= (0E/0c) (dc, +dc2) =0. (4) 


It should be borne in mind that the differentials 
dc, and dce occurring in (3) and (4) are related 
to each other in a way depending upon the 
geometry of the system, i.e., on the specific 
values assigned to @), @2, 5), b2, and also involve 
the Morse curves of the individual atom pairs. 
But from (4) it follows that a stationary value of 
E corresponds to a stationary value of c, ignor- 
ing the unimportant singular cases where 
dE/dc=0. Whether a minimum of £ corresponds 
to a minimum or a maximum of ¢, is determined 
by the sign of dE/dc, which in the particular 
case of the London formula can be written as: 


dE /dc=rA—(1—A)(c—3(a+d)) 

XL2(a—b)?+(a—c)?+(b—c)?) 4. (5) 
This is negative if |c| is the smallest of the 
bonds, as it is for the comparatively small values 
which XX assumes for nonmetallic elements. 
Thus we have shown that, to make E as low as 
possible for a given a and |, the value |c| must 
be made a minimum, and thus the best reaction 
path lies in the boundary surface. Equations 
exactly similar to (3) and (4) can be discussed 
for a and b. If |a| happens to be the largest of 
the three, then dE/da is positive so that a 
minimum value of E will now go with a maximum 
value for |a|. In some cases it is found more con- 
venient to use this latter relationship. In the 
intermediate case where |b] say, is neither the 
largest nor the smallest of the three |b| should 
be made a minimum or a maximum depending 
on whether 0£/0b=0. A mechanical device 
which enables one easily to satisfy condition (3) 
will be described in a later section. 


A VEcTOR ADDITION RULE FOR BOND ENERGIES 


After this more or less kinematical discussion 
of configurations in bond space we now turn to 
the question of the interaction energy of a system 
of atoms. In particular we shall consider the 
nature of the equipotential surfaces defined by 
(1) and (2). We now prove that the law of 
combination of exchange energies as expressed 
by the square root in (1) is equivalent to the 
following simple vector addition: 


E.z= (3[(a—B)?+(8—y)?+(y—a)?])! 
=latwit+w’y|, (6) 


where w and « are third roots of unity. This 
situation is also represented by the vector 
diagram Fig. 2 in terms of which the following 
simple proof can be given: If we plot the bond 
energies a, 8, y end-on making angles of 60° with 
each other as in Fig. 2, then the resultant has the 
horizontal and vertical components (8 — }(a+v7)) 
and 3 X3*(a—vy). Thus its length turns out to be 
just E,,. The existence of such a simple vector 
addition rule for bonds in three and four particle 
systems is significant and suggests the possibility 
of similar relations in more complicated cases. 
Its practical importance lies in the avoidance of 
practically all calculations in obtaining potential . 
energies from the individual bonds, especially 
since the 60° triangular paper on the market is 
ideally suited for such purposes. 

The equipotential surfaces in bond space like- 
wise turn out to be strikingly simple. Trans- 
forming to new orthogonal coordinates by means 
of a rotation of the a, b, c axes, 


x=1/4/2(a—b), 
y=1/s/6(a+b—2c), (7) 
g=1 ‘4/3(atb+o), 


we find the expression for the energy of the 
system using (1) and (2) becomes: 


E=v/3(dz—(1—))(3(x?+3”"))! 
=4/3[Az—((1—A)/+/2)r]. (8) 


Thus it is seen that the equipotential surfaces 
are circular cones with their axes along the 2 
axis, and with their apex angle 2 tan~ \/(1—)) 
depending only on the ratio of additive to ex- 
change binding. The relation exhibited by (8) 
that E,, is simply proportional to the distance r 
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from the z axis provides an alternative proof of 
our previously given vector addition rule since 
the three valence bonds a, 6, y when plotted 
along their axes, a, 0, c appear to make just the 
angles used in Fig. 2 to an observer looking 
down the z axis. The equipotential surfaces in 
bond space together with the boundary surface 
B determine completely the best reaction path: 
The intersection of the equipotentials with B 
traces a system of equipotential lines on B of 
much the same character as the system of 
equipotential lines found in the usual plots of a 
potential surface in configuration space for a 
three-atom problem. 

In many triatomic reactions the smallest bond, 
c, remains very small throughout so that the 
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surface B can be approximated by the plane 
c=0. In this case the equipotential curves on B 
are simply a system of confocal ellipses generated 
by the intersection of the ab plane with the 
cones. The reaction path is then the line where 
either a or b has its most negative value, c being 
zero. 

The other extreme where |c| assumes com- 
paratively large values for three-atom reactions 
is typified by the case of three H atoms. In Fig. 3 
we have given the equipotential curves plotted on 
the B surface as they appear looking along the 
y direction in bond space. This particular pro- 
jection was chosen because of the ease with which 
the construction can be carried out. 

The potential diagram for c=0 is practically 
identical with the one for 3 H atoms in the region 
lying below the best reaction path indicated by 
arrows in Fig. 3.* It has seemed desirable to give 
this brief discussion of two typical three-atom 
cases since the general features are just the same 
as for four-atom reactions. 


LOCATING THE ACTIVATED CONFIGURATION OF 
Four ATOMS 


_It is in the application to four-atom reactions 
that the use of bond space instead of configura- 
tion space proves to be of real advantage. Even 
in the presentation of a symmetric case such as 
H2+Is, the method of finding the activated 
state as well as the representation of the po- 
tentials in configuration space becomes ex- 
tremely laborious. With the present method, a 
four-atom system with no symmetry at all such 
as H.+ICl has been investigated and the re- 
action rate calculated with no unreasonable 
hardships. 

A mechanical device for selecting configura- 
tions which, for a given pair of values a, b, have 
the smallest possible |c| to go with it, is pictured 
in Fig. 4. Scales are provided on which the six 
interatomic distances between four atoms H, H, 
I, Cl are marked with their respective bond 


3 As pointed out in reference 2, we are ignoring through- 
out configurations involving ‘‘contracted” bonds. However, 
if one wants to extend the B surface into this region, it 
results in making multiple-valued the present allowed 
portion of the B surface in the case when c=0. In the case 
c0, the extended B surface bends over completely in the 
direction of the z axis, crosses it and finally departs in the 
direction of infinite positive z. 
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_ Fic. 4. Device for rapidly finding configurations which lie on the boundary surface, 
i.e., having the lowest sum, c, for the diagonal bonds for a given a and b. The activated 


configuration is the one shown. 


energies and the scales so obtained attached to 
each other so as to present a model of the plane 
configurations of these four atoms. The scales 
are hinged by means of thumb-tacks in such a 
way that, without varying four of the six dis- 
tances, the remaining two, let us say the di- 
agonal ones, can be adjusted until the sum of 
their energies (absolute values) as read directly 
from their scales is a minimum. In doing this 
one must of course try various combinations of 
a; and a2 which add up to the same a, and also 
various combinations for b, until an absolute 
minimum for |c| is obtained. By way of em- 
phasis we should like to point out that this 
method of variation leads to a unique configura- 
tion for fixed values of a and 6 thus establishing 
our previous statement that there exists a 
point-to-point correspondence between those 
points in three-dimensional bond space which lie 
on the surface B, and the points on the cor- 


responding surface in six dimensional configura- 
tion space. 

A point of the B surface can be found in a few 
minutes with the use of this device. A sufficient 
number of such combinations (a, b, c) must be 
found to map out the potential diagram on the 
B surface. A very pleasant feature of this pro- 
cedure is that it enables us to describe the es- 
sential aspects of the potential surface in two- 
dimensional diagrams in spite of the fact that the 
number of parameters of the dynamical problem 
is six. The plot for the reaction, H2+ICI—-HI 
+HCl is given in Fig. 5, showing the projection 
on the ab plane of the intersections of the bound- 
ary surface with the equipotential curves. An 
alternative way which is sometimes more satis- 
factory is to plot the xy projection. 

Since the B surface does not depend on a 
particular potential law, the points (a, b, c) de- 
termined in the way described may be used for 
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Fig. 5. Projection on the ab plane of the intersections of 
the equipotentials with the boundary surface B for the 
reaction H.+ICI—HI+HCI. 


different fractions \ of additive binding and 
indeed for different laws of interaction and the 
results compared with experiment. In this 
particular case we have located the activated 
point and determined the activation energy for 
14, 17 and 20 percent of additive binding, with 
results given later in Table II. 

The constants used in the Morse curves for the 
various bonds are given in Table I where the 
Morse curve is given by 


E=D'[e-2(r-10) _ 2e—2(r—r0) J, 


TABLE I. 








ro (A°) wo(cm=) D (kcal.) D’ (kcal.) a 


H—-H 0.74 4375 102.4 108.6 1.941 

I—Cl 2.310 384.6 49.4 49.95 1.871 
H-Cl 1.282 2840.8 101.5 105.5 1.781 
H-I 1.616 2233 70.0 73.2 1.700 
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There is nothing in our procedure which 
limits us to the constant percentage rule as ex- 
pressed by the constant \ in Eq. (2). We might 
equally well have replaced each of the energy 
scales in our model by a double scale giving the 
additive and the exchange part of the bond as 
separate functions of the distance. 

The accuracy with which the activated con- 
figuration can be found using this easily impro- 
vised mechanical model (of a size no larger than 
is convenient) is more than sufficient for all 
practical purposes. But if greater accuracy is 
required, a rapidly converging mathematical 
procedure can be devised which renders the 
question of attainable accuracy merely one of 
how much numerical computation one is willing 
to bestow on the solution. Suppose that from the 
mechanical model we have obtained a configura- 
tion close enough to the activated one so as to 
make all the first derivatives 0E/0q; small, and 
suppose that we are in possession of a complete 
set of second derivatives 0?E/dqi0qx, taken with 
respect to all five “‘plane’’ coordinates qj, for this 
configuration. Knowledge of these second deriva- 
tives is necessary anyway for the calculation of 
the characteristic vibrations, and the method 
for obtaining them will be described later in that 
connection. Denoting by Ag; the differences 
of the coordinates of the true activated con- 
figuration minus the ones for the approximate one 
which was found, we can expand the first 
derivatives 0E/dq; in the neighborhood of the 
activated point where (0£/0q;),i,=0, as follows: 


(0E/0Qi) gie—a9; = (0 E/0qidgx)(—Agx). (9) 
k 


Solving this set of inhomogeneous linear equa- 
tions for the quantities Aq, and using five arbi- 
trary quantities \; as multipliers, we can express 
our solution in a particularly condensed form as 
follows: 


—drAAq. IF/dqi 0E/dq2 0E/dqs 
1 @E/dqi0q1 3E/dq10q2 #E/dq10qs 
D= Xe =0 (10) 
As 0?E/dq;0q1 0° E/0q;09s 
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TABLE II. 
H-H I-Cl H-Cl H-I H-Cl H-I _ 

Config.* 7; rs rs 14 rs r¢ 0E/dq, AE/dq2 dE /0q3 0E/0q4 0E/0q5 
i 1 8033 2.3359 1.597 2.264 2.214 2.475 2.966 —4.022 — 3.831 — 1.519 7195 
ia 2 8048 2.354 1.597 2.200 2.195 2.460 395 — .399 .158 .200 — .513 
an 1 791 2.346 1.564 2.372 2.166 2.548 1.071 .229 — 1.050 — 4.143 _ 1415 
ies 2 792 2.345 1.564 2.330 2.190 2.508 .037 — .06 .80 48 .06 
4% «#41 «8300239 ~=S1.61si208—si2ks 








* 1 and 2 refer to the preliminary and the improved configuration. 


From (10) a particular Ag, is determined by 
putting the corresponding \, = 1 and all the other 
\’s equal to 0. For practical computation it is 
very convenient to first reduce the determinant 
D to diagonal form except for the first row and 
the first column, thus obtaining all quantities 
Ag; at once. These five quantities define a 
certain change in the originally assumed con- 
figuration, a change which can be carried out 
graphically with great ease and precision. For 
this new configuration the first derivatives 
dE/dq; may be recalculated and the improve- 
ment obtained in locating the activated point 
may be judged by the decreases simultaneously 
obtained in all five derivatives. The process has 
been applied to the activated configurations of 
H.+ICl both for \=17 percent and 20 percent, 
with results tabulated in Table II. The simul- 
taneous decrease of all five derivatives by ap- 
proximately a factor 0.10 must be very striking 
to anyone who tries to improve the approxi- 
mate configuration in a less systematic fashion. 
The activation energies E, (see Fig. 6) for 14 
percent and for 20 percent were found to be 41.2 
and 24.9 kcal., respectively. By linear inter- 
polation a percentage of 17 was indicated as 
giving E.=33. The experimental value of E is 
approximately 34 kcal. The activated point for 
17 percent was then actually determined and 
gave an activation energy E.=32.8 kcal. The 
activated configuration for 14 percent was not 
determined with the same accuracy as for the 
other two cases. Sherman and Li*® have calcu- 
lated a value of 24.6 kcal. using 14 percent 


ee 


asa Gore and Yost, J. Am. Chem. Soc. 57, 2723 
5). 


* Sherman and Li, J. Am. Chem. Soc. 58, 690 (1936). 


Distances are given in A units. 


additive binding in agreement with the value 
given here. 


THE VIBRATIONS OF THE ACTIVATED STATE 
(VIBRATIONS OF PLANE CONFIGURATIONS 
OF Four UNLIKE ATOMs) 


Six of the twelve modes of motion of four 
atoms arranged in a plane involve a change in the 
relative positions of the atoms while the other 
six amount to rigid translations or rotations of 
the whole molecule, and proceed without change 
in the potential energy. Of the ‘“‘internal” 
motions, five preserve the plane arrangement of 
the four atoms whereas the sixth involves dis- 
placements different in sign for adjacent atoms 
and at right angles to the plane of equilibrium. 
Since this plane is obviously a plane of symmetry 
for the molecule, this last mode belongs to a 
symmetry class different from the rest and the 
secular equation will split into two factors of 
first and fifth degree, respectively. There are no 
other symmetry conditions imposed on the 
activated configuration as there is in reactions 
of the type H2+I.—2HI where the secular 
equation decomposes into 1+2+3 degrees of 
freedom. ® 

We consider first the vibrations where the 
plane arrangement is preserved but ignore 
translational motions in this plane by keeping 
the center of gravity fixed at the origin. Tenta- 
tively we introduce 6 parameters the meaning 
of which is apparent from Fig. 7 and the follow- 
ing remarks: q3 is the distance between the two 
centers of gravity of m,; with m2 and m; with my. 
One of these centers of gravity separates g; into 


6 Wheeler, Topley and Eyring, J. Chem. Phys. 4, 178 


(1936). 
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Average energy of activated complex 
at fem rature of reaction] 
at absosolute zero 
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Fic. 6. Activation energies: E is the usual difference in 
average energies at temperature of reaction; Ey=activation 
energy at absolute zero; E,=‘‘classical’”’ activation energy 
at absolute zero. All are useful quantities in the theoretical 
treatment of reaction rates. 


two parts in the ratio p : (1—p) and the other 
separates ge in the ratio o : (1—«a). These six 
coordinates not only describe the relative posi- 
tions of the four particles but also the angular 
position yg of the whole system relative to a 
reference line R fixed in the plane. Thus the 
kinetic energy, when expressed in all six param- 
eters, includes a part arising from rotations in the 
plane which can be separated from the rest 
because such a motion involves no changes in the 
potential energy. 
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Fic. 7. Notation used in treating four atoms. 


The final coordinates were chosen in such a 
manner as to give a sum of square terms for the 
kinetic energy. This determines the coordinates, 
of course, only to within an orthogonal trans- 
formation, but the rotational coordinate ® is 
uniquely determined by the requirement that it 
be normal to all the other coordinates none of 
which can involve the parameter ¢. 

The expression for the kinetic energy in terms 
of the parameters defined in Fig. 7 is easily 
calculated to be: 


2K =mi2g 2+ mage? + mirage tag?+b(2¢a+a*) +c(298+6%), (11) 


where m2, ™34 and m2 34 are reduced masses of m, with mo, m3 with m4, mi+me with m3+m, 


respectively, and where a, 0, c stand for 


4 


a= > m:R?, b=qPmw, c= Q2?mga. 


1 


The R’s measure distances from the center of gravity of the system. The last three terms of (11) are 
easily transformed into a sum of squares if we substitute 


w= g(a—b—c)!, b=b"a+¢), w=c(B+¢), 


which property is preserved under any orthogonal transformation, in particular the following one: 


U=a-?((a—b—c)'u+b4v+ciw), 


V=[(a—b—c)/(a(b+c)) JEL —(b+c)/(a—b—c)'u+b+clw], (12) 


W=(b+c)—(clvo—biw). 


Substituting for these variables in terms of the original parameters ¢, a, 8 and also for a, b, ¢ we 


obtain: 


9 





2K= M291? -+M34g2? + m aude +| tina bos _ 


sieht inal _, Gr gem i2mss - 
ee ee 


qs? ++———, (13) 


> m:R? qm 12+4G2"M34 mR? 
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where a= kat (1 —k)pg (qi?2m 120+ 42M348) (q:?m12+G2?m4)—, qs=a —£, 
=> m:R29+g2myat+qemsf. (14) 


The newly defined coordinates g, and q; do not involve the angle ¢ and are thus orthogonal to the 
rotations of the system as a whole. The rotational coordinate is thus ®; since it does not enter into 
the potential function for the system, it may be omitted from further consideration. 

As was stated previously there exists but one internal degree of freedom which carries the four 
particles out of their plane of equilibrium z=0. Indeed, of the four independent motions involving 
displacements 2; at right angles to the equilibrium plane we may choose two as the rotations 7, and 
r,, and a third as the translation ¢., and represent them by the following coordinates: 


r2= OM BiXi= iM12 Sin a(Z1—Z2) + Goma, Sin B(Z3—24), 
ry= DL misiyi=(qs/ domi) { — (msg) (121+ MeZ2) + (m1+me) (mMs323+ M424) }, 
{.= > mizi. 


The internal motion is then completely determined by the requirement that its coordinate g¢» be 
orthogonal to 7,, ry, fz. Let gs be a linear function of the z displacements gg= )/;2;, then the following 
relations must hold among the o’s: ‘ 

4 


4 4 
2. m;o;x,;=0, = m;oyi:=0, :. m,o,=0. 


=1 i=1 i=1 


Thus we find by a simple elimination process, that gs is proportional to: 




















21 Z2 23 24 
qim,2 Sin @ —gim2 sina q2M34 sin B —g2m;34sin B 
q~™ m3+m, ms+m, mM +My m,+my2}- 
iM 12 COS a@— M1g3 — iM 12 COS @— Moz G2M34 cos B+m3q3 —Q2M34 cos B+ m4g4 
ys > 
ams am; “mj 2m; 
my, Me ms m4 


One possible way of describing the internal motion connected with g¢ is the following : 


























cot a—cot B 1 : | —cot a—cot B 1 7 
Az: = Ade - ii Az3= Ade _ . ‘ 
i q3(mitme) gim,sina | g3(m3-+ms) g2m3 sin 6 1s) 
cot a—cot B 1 7 r cot a—cot B 1 
Az2= Ade _ ,» Azys=Agqd — 4 
L gs(mit+ms) gqumesin a] L = g3(ms-+m4) gam, sin B! 


The displacements (15) comply with the three conditions: 
> m,x;Az;=0, > myy:Az;=0, > m,Az;=0, 
as can be readily established. Thus they are orthogonal to the three rigid motions (Ar,, Ar,, Afz), and 


of course also to the other modes which do not involve z,. Using (15), the kinetic energy Kg associated 
with gs is then 


9 


cot a—cot B\? 1 1 1 21 1 2 1 
nef Y(t al 
qs M,;+M. M3+m, gisina/ my. gzsinB/ mg, 














The next problem is the computation of the force constants of the activated configuration. Since 
the activated configuration, like a stable configuration of atoms, has a stationary value of the inter- 
action energy, it is actually in stable equilibrium in all degrees of freedom except the one along which 
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the reaction proceeds. That means that all first derivatives 0E/dq_ are 0 and that for small displace- 
ments from the activated configuration the energy function can be expanded : 


OE 





dE=>, gad gs. (16) 
aB 0409p 


The second derivatives in (16) are force constants for small displacements from equilibrium. When 
transformed into a sum of squares, only one force constant in (16) will exhibit a negative sign. 

It has been preferable here to follow an explicit analytical procedure in deriving numerical values 
for the force constants rather than the previously followed methods which though applicable in 
simpler cases become excessively laborious for more degrees of freedom. One other fact should be 
borne in mind. The present method of second derivatives gives the force constants for infinitesimal 
vibrations so that in case of considerable anharmonicities they may not give the correct frequencies of 
finite vibrations. 

Formulas necessary for the calculation of the second derivatives are presented in Tables III, 
IV and V. Formulas in IV and V are of course perfectly general and depend only on the geometry of 
the quadrilateral (Fig. 7). They give only first derivatives since tabulation of all second derivatives 
would require too much space ; however, calculation of the second derivatives 0?7;/0qa0qg from IV, V 
requires no further consideration of Fig. 7. 

The second derivative 0?E/dq,? remains to be calculated. We have 














CE @E dr; dr; |= = dr; 0°3i| 
— + , (17) 
dz; age! 


a + panama 
ge i Or,Or, Ags Og, * Or; | d22\ dq¢ 


which reduces to 


—(Zi2— (18) 


CE OE da; 1f d P 
rt own 
dqe 


0g6" i Oa; Or; 7; 


since all the derivatives 0r;/0q5 are zero on account of the symmetry character of this mode of motion 
(all the r’s are minima in the plane arrangement) and since 


02; 0?r; 1 
=0 and =—, 
0g¢ 027 7; 








Here the 2; in (17) refer to the z projections of the distances r; and the quantities 2;; and 2;2 in (18) 
refer to the two end-points of the interatomic distances 7;. Thus the force constant is given by (18) 
and all the mixed derivatives involving gs vanish identically. 


THE VIBRATION FREQUENCIES FOR THE REACTION H2.+ICI—HI+HCl 


In evaluating the terms for the kinetic and the potential energy using 17 percent additive binding 
we are led to the following secular equations for the five normal vibrations in a plane: 


(402.82 —.25X) 5.29 93.30 107.33 87 
5.29 (246.10 —13.86X) 5.72 34.54 7.78 
93.30 5.72 (—62.70—.9879X) — 78.64 — 20.12 =0. 
107.33 34.54 — 78.64 (112.37—3.99.X) 39.50 
87 7.78 — 20.12 39.50 (29.86—.161X) 


The twisting vibration out of the plane is given by the equation 272.88 — 1.5788X =0. Here the force 
constants are in units of kcal. per mole A~? and in the kinetic energy mass is expressed in molecular 
weights and distance in angstroms (where distance occurs). The six characteristic frequencies »; for 
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TABLE III. 
E=x(a+b+c)—(1—A)R. 
Here R=[3((a—6)?+(b—c)?+ (c—a)?) } and a@=a\+a2, b=bi+b2, c=c1+C2 
ee OE 2 i 7 ee i - % Jj 3(1- —_ 
PE _ydE Or | y PE ari an dE_dE_, _3( ill a4b+) 
0ga0Qg i OF; OGaOQg ik OF: OTK AQa OB 0a, dae 2 & 2 R 
dE _ dE da; _ dE da; #E_@E_ ®E __31=2dy_ 
Or; Oa; dr; da; dr, 0a;2 ax? Oa,da. 4 R® 
PE _ PE daida, , dE Fai &E =-ja-»| “a -2]- 3 (1—)) =1.2 
Or;dr, Oa:0ax dr; dr, da; Arr: dajoh; 4 R3 R| oac2?'4 R sine. 
TABLE IV. Formulas dr;/dq; valid for any configuration of four atoms in a plane. Notation as in Fig. 7. 
| 71 q2 q3 4 | 5 
rn | 1 0 0 0 | 0 
rs | 0 1 0 0 | 0 
rs —p Sin €; o sin ¢3 sin (a—e;) ds Cos (a—€;) | —p(1—k)r; cos e;—akr2 cos £3 
14 | —(1—p) sine, | (1—c) sin sin (a+e4) ds Cos (a+e,) | (1—p)(1—k)r; cos eg + (1—a)kr2 cos &4 
rs|  (1—p) sin €5 o sin ¢5 sin (a—e;) gs Cos (a—e;) | (1—p)ri(1—k) cos e;—okr2 cos £5 
16 psine,| (1—o) sin f sin (ate) g3 cos (a+es) | —p(1—k)r: cos eg +(1—a)kre cos &% 
a | 
TABLE V. Derivatives of ¢:,, a and B with respect to the coordinates qx. 
| qi q2 | q3 4 q5 
és | —»- cos «; ? cos 3 os cos (a—e;) 2 sin (a—€3) : [eri(1—k) sin es+orek sin ¢;]+1—k 
r3 rs 13 13 rs 
—p —7o 1 93. " 1 : a 
€;| ———cosg cos ¢%4 — cos (a+e) | —= sin (a+e.) | ——[(1—p)ri(1—k) sin eg—(1—o)rok sin £4] 
"4 "4 "4 14 "4 
—1+k 
| —p a 1 qs. 1 , ; 
€5 —— COS €5 — cos ¢5 —— cos (a—€;) =" sin (a—es) | ——[(1—p)ri(1—k) sin est+orek sin f;]+1—k 
’s ’5 5 ’s ’s 
p —o 1 q3_. 1 ‘ , 
€6 — COS €% cos ¢¢ — cos (ates) | —= sin (a+e¢) [pri(1—k) sin eg —(1—o) rok sin f6 ]—1+k 
’6 6 ’6 6 ’6 
| 0 0 0 1 1—k 
8 | 0 0 0 1 1—k 
The corresponding derivatives for the angle ¢; follow from the same table as the «; by use of the identities: ej —{ij = +(a —8) i=3, 4, 5, 6, where 
the plus sign goes with i =3, 5 and the minus sign goes with i =4, 6. The geometrical significance of the symbols is given in Fig. 7. 
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TABLE VI. 
Moments of Inertia Frequencies of char. 

Activated (gram atom, A?) vibrations Act. 
configura- (cm) Energy F, 
tion for A 3 B; C3 V} ve V3 V4 V5 Vn (keal.) 
17% 152.95 5.81 158.76 4391. 1580. 1426. 481. 435. 41128 32.8 
20% 158.54 6.00 164.54 4696. 1522. 1344, 550. 442. 1933 24.9 

additive 


binding 








17 percent additive binding are given in Table VI in wave numbers. They are obtained from the six 
roots X; of the secular equations by the simple relationship v;= 108.46X 4. 
The corresponding secular equations for 20 percent additive binding are: 


(465.07 —.25X) — 2.48 76.26 81.46 —2.11 
—2.48 (252.68 — 13.86X) 11.65 38.50 6.60 
76.26 11.65 (—47.14—.9878X) —46.41 — 19.00 =0 
81.46 38.50 —46.41 (148.57—3.99X) 34.62 
—2.11 6.60 — 19.00 34.62 (26.96 — .1565X) 


and 250.08 — 1.631X =0. 


THE RATE OF REACTION 


The general expression for the rate of a reaction of this type has been given.”: ® Thus we have the 
equation 








m,'m-! 03 T I,'T2! 5 hv; 
II (2 sinh —) 
2kT 


1 hv» 2 E. 
x(1-—(*) ) exp (- ). (19) 
24\ kT kT 
Here the subscripts 1 and 2 refer to the hydrogen and iodine chloride molecules, respectively, and 
3 to the activated complex. «x, m;’, ¢; and v; refer to the transmission coefficient which we take equal 
to unity, the mass of thezth molecule, the rotational symmetry number and the frequency respectively. 
A;B;C; are the three principal moments of inertia of the activated complex and I,’ and I,’ the mo- 
ments of inertia of the hydrogen and iodine chloride molecules respectively. v, is the imaginary mode 
of vibration normal to the barrier. The primes indicate that masses are to be taken in molecular 
weight units and lengths in angstroms in which case k; is given in the units cc moles“ sec.—. 
If; as is usual we define the energy of activation as 


E=kT"d log (k/T*)/dT 
5 hv; hy hv 3 


we have E=Eot+> _ - —-kT 
i=1exp (hv;/kT)—1 exp (hvi/kT)—1 exp (hve/kT)—1 2 





hy, hve 
(2 sinh —) (2 sinh ) 
b ( m3, ) 7102 5.785 X10" (A3’B3'C3')3 2kT 2kT 
=< 








5 
and Eop=E.4+>, shv;—4(hyy+hre). 
i=1 


7 Eyring, J. Chem. Phys. 3, 107 (1935), Eq. (7). 














1X 


he 


nd 


ual 


no- 
yde 
lar 





We thus obtain for 17 percent additive binding 
E,=32.8; Eo=37.85 and E=36.62 at 503.2°. 
For 20 percent additive binding E,=24.9; 
Ey= 30.29; E=29.04. Now defining three A’s 
by the following equations 


kR=Ae-EelkT = A ye EBlkT = A e-EIkT 


we obtain Table VII. 

These values are to be compared with the 
experimentally obtained specific reaction rate 
of Bonner, Gore and Yost.‘ We change k 
from their units to cc mole sec.-' and so 
obtain k=.52=1.64 K 105-33 ,900/k478.2. 2.8 = 1.49 
4 10!5e—38 ,900/k503 .2 and 5.8 = 1.60 4 10 !5e—38 ,900 /k513 .2 

By inspection of Table VII we see that A is 
independent of whether 20 percent or 17 percent 
additive binding is used. This is important since 
this means it is uniquely fixed by the potential 
function and depends on no arbitrary constant. 
If the percentage additive binding were raised to 
a fraction over 20 percent E would be lowered by 
a calorie to 28 kcal. and k would then agree 
rather well with the experimental values. The 
theoretical k would, however, not increase with 
temperature as fast as the experimental value 
does. Thus, if Bonner, Gore and Yost are, as they 


DIAMAGNETIC ANISOTROPY 


TABLE VII. 























% additive a k Ac X 1019 | AoX1013 | A XK 1012 
17 478.2| 7.441075 7.36 1.50 4.11 
17 503.2 | 5.12 10-4 9.06 1.46 4.16 
17 513.2 | 1.03 x 10-8 9.64 1.37 4.09 
20 478.2 | 2.2810 5.51 1.60 4.30 
20 503.2. 1.01 6.61 1.46 4.17 
20 513.2 1.79 7.26 1.43 4.21 








believe, actually measuring the reaction rate here 
calculated and if their values of k are dependable 
then our calculated A is too small by a factor of 
=~ 350. A would be increased by this amount if 
the three lowest frequencies are actually around 
50 wave numbers instead of the calculated values. 
Such a change would also reduce the excess 
(E,—£.) of the zero-point energy in the acti- 
vated state over that in the initial state to 
practically zero. More experiments seem neces- 
sary before any of these quantities can be as- 
signed from experiment with assurance. The 
corresponding isotopic reactions will also be 
useful in calculating E»)—E,. Extensive accurate 
data on reactions of this type will be particularly 
valuable in determining the accuracy and in sug- 
gesting methods for improving the calculation of 
potential functions of activated complexes. 
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The Diamagnetic Anisotropy of Aromatic Molecules 


Linus PAULING, Gates Chemical Laboratory, California Institute of Technology, Pasadena* 
(Received July 11, 1936) 


Values of the diamagnetic anisotropy of benzene and other aromatic hydrocarbon molecules 
are calculated on the basis of the assumption that the p, electrons (one per aromatic carbon 
atom) are free to move from carbon atom to adjacent carbon atom under the influence of the 
impressed fields. When combined with the assumed values for the contributions of the other 
electrons (—2.0X 10~® for hydrogen, —4.5 X 10-6 for aromatic carbon, —6.0X 10° for aliphatic 
carbon) these lead to principal diamagnetic susceptibilities of molecules in approximate 
agreement with the available experimental data. The diamagnetic anisotropy of graphite is 


also discussed. 


Pres diamagnetic crystals, such as 
bismuth and graphite, show pronounced 
diamagnetic anisotropy, the susceptibility being 
of normal magnitude in the basal plane and 
abnormally large in the direction perpendicular 
to this plane. A similar phenomenon is shown by 


* Contribution No. 555. 





crystals of aromatic substances. It has been 
pointed out by Krishnan and his collaborators! 
that the magnitudes and orientations of the 
principal axes of the susceptibility ellipsoid of 

1(a) K. S. Krishnan, Proc. Ind. Sci. Cong., Madras 
Session, January, 1929; (b) K. S. Krishnan, B. C. Guha 


and S. Banerjee, Phil. Trans. Roy. Soc. A231, 235 (1933); 
(c) K. S. Krishnan and S. Banerjee, ibid. A234, 265 (1935). 
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such a crystal are determined by the values of 
the principal susceptibilities of the molecules of 
the aromatic substance and by the orientations 
of the molecules relative to the crystal axes. 
The susceptibility ellipsoids of the aromatic 
molecules are found to be approximately prolate 
ellipsoids of revolution, with the long axis 
normal to the plane of the ring system. As a 
result of the valuable investigations of Krishnan 
and his collaborators this phenomenon has been 
made the basis of an important auxiliary method 
for determining the structure of crystals.': ? 

A qualitative explanation of these abnormally 
large diamagnetic susceptibilities as arising from 
the Larmor precession of electrons in orbits in- 
cluding many nuclei’ has come to be generally 
accepted. With the aid of simple assumptions, 
I have now developed this idea into an approxi- 
mate quantitative treatment, described below. 


Benzene 


The Hamiltonian function for an electron in a 
constant magnetic field of strength H parallel to 
the z axis is‘ 


K= H°+5C/ +50", (1) 


in which 3° is the Hamiltonian function for zero 
field strength, 


5’ = — (He/2mc)(xpy—ypz), (2) 


and 5C”’ = (He? /8mc*) (x*+”). (3) 


The contribution of the electron to the diamag- 
netic susceptibility of the system can be calcu- 
lated by the methods of quantum-mechanical 
perturbation theory, a second-order perturba- 
tion treatment being needed for the term in 3’ 
and a first-order treatment for that in 3¢”. In 
case that the potential function in 3C° is cylin- 
drically symmetrical about the z axis, the effect 
of %’ vanishes, and the contribution of the 
electron to the susceptibility (per mole) is given 


2S. Bhagavantam, Proc. Roy. Soc. A124, 545 (1929); S. 
Banerjee, Ind. J. Phys. 4, 557 (1930); K. S. Krishnan, 
Nature 130, 313, 698 (1932); K. S. Krishnan and’ S. 
Banerjee, ibid. 131, 653 (1933); 132, 968 (1933); Zeits. f. 
Krist. 91, 173 (1935) ; K. Lonsdale, Nature 137, 826 (1936); 
N. Ganguli, Zeits. f. Krist. 93, 42 (1936). 

3 P, Ehrenfest, Physica 5, 388 (1925); Zeits. f. Physik 58, 
719 (1929); C. V. Raman, Nature 123, 945 (1929); 124, 412 
(1929); C. V. Raman and K. S. Krishnan, Proc. Roy. Soc. 
A113, 511 (1927). 

4 See, for example, J. H. Van Vleck, Electric and Magnetic 
Susceptibilities (Oxford University Press, 1932). 
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by the Pauli expression, resulting from 3”, 
+ edie (Ne*/4mc*) (p") avy (4) 


in which (p?),, is the mean square of the distance 
of the electron from the z axis. This expression is 
valid for atoms and monatomic ions, the 2 axis 
being passed through the nucleus. It is not ap- 
plicable (without correction for the term in 5’) 
to molecules. For most molecules a satisfactory 
approximate treatment can be used _ which 
consists in assigning electrons to the nearest 
nucleus, Eq. (4) being applied to the electrons 
about each nucleus. 

Of the forty-two electrons of the benzene 
molecule, twelve (the carbon K electrons) make 
no appreciable contribution to the magnetic 
susceptibility. Twenty-four of the LZ electrons 
(which may be considered to be involved in elec- 
tron-pair bonds between adjacent atoms) are in 
orbitals which are symmetric with respect to 
reflection in the plane of the nuclei. It is probable 
that the contribution of these to the magnetic 
susceptibility is normal. The remaining six L 
electrons, which give to benzene its characteristic 
electronic structure and properties,> occupy 
orbital wave functions which are antisymmetric 
with respect to reflection in the plane of the 
nuclei. The probability distribution function for 
these electrons is large only in two ring-shaped 
regions, one above and one below the carbon 
hexagon. We may well expect that in these 
regions the potential function representing the 
interaction of an electron with the nuclei and 
other electrons in the molecule would be ap- 
proximately cylindrically symmetrical with re- 
spect to the hexagonal axis of the molecule, the 
electron, some distance above or below the plane 
of the nuclei, passing almost imperceptibly from 
the field of one carbon atom to that of the next. 
Under these conditions the contribution of each 
of these six electrons to the susceptibility with 
the magnetic field normal to the plane of the 
nuclei would be given by Eq. (4), with p meas- 
ured from the hexagonal axis of the molecule, 
and it would hence be very large. On the other 
hand, because of the nodal plane of the wave 
functions these electrons would make only 2 


very small contribution to the susceptibility 


5 E. Hiickel, Zeits. f. Physik 70, 204; 72, 310 (1931); L. 
Pauling and G. W. Wheland, J. Chem. Phys. 1, 362 (1933). 
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with the field in the nuclear plane. There would 
thus arise a large diamagnetic anisotropy of the 
molecule. If we use for (p?),, the value R?, in 
which R is the distance from the axis to the car- 
bon nuclei, 1.39A, the anisotropy is calculated to 
have the value AK = —49.210-*, which is in 
reasonably good agreement with the experi- 
mental value for benzene,'» AK = —54x10-. 

A more detailed consideration may be given 
this argument by evaluating the terms cor- 
responding to KX’ and 3’’ with the use of as- 
sumed expressions for the wave functions. 
Using hydrogen-like 2, wave functions, the 
term corresponding to ,3¢’”’ (Eq. (4)) has the 
value —49.210-* for large values of the effec- 
tive nuclear charge, the absolute value increasing 
with decrease in the effective nuclear charge, as 
a result of the greater contribution of large radii 
to (p”)a,. The term in 3’ can be represented by 
means of the approximate second-order per- 
turbation theory® as 


x’=([(xpy—ypz)? Joo/((Wm)ar— Wo), (5) 


in which the numerator represents the diagonal 
matrix element for the normal state under 
consideration, W» the energy of this state, and 
(W.)av an average energy for significant excited 
states. It is found on evaluating this expression 
that it is small for reasonable values of the 
effective nuclear charge. However, this ex- 
pression as well as that of Eq. (1) are rather 
sensitive to changes in the assumed wave 
functions, and there is little prospect at present 
of carrying out an accurate theoretical calcula- 
tion of the susceptibility. 


Condensed aromatic ring systems 


We may now calculate the magnetic anisotropy 
of other aromatic hydrocarbons in the same way 
as for benzene, making the assumption that the 
2p. electrons are free to move under the influence 
of the impressed fields from carbon atom to 
adjacent carbon atom. Our problem is equivalent 
to that of calculating the magnetic effect of the 
currents induced in a conducting network. In 
benzene we have assumed the induced current to 
flow from atom to atom along circular arcs. In 
large condensed ring systems, such as graphite, 
the current would be assumed to flow along 


*L. Pauling and E. B. Wilson, Jr., Introduction to 
Quantum Mechanics, Sec. 27e. 
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Fic. 1. Fic. 2. In- Fic. 3. As- 
duced currents sumed path of 
in anthracene. induced currents 

in anthracene. 


straight lines between adjacent atoms. In order 
to simplify the calculation we shall first consider 
rectilinear hexagonal nets, and then make a 
small correction by replacing the outermost lines 
by circular arcs. We assume in all cases that the 
electrons are uniformly distributed along the 
segments of the network, the linear density of 
electrons being co. 

Let us consider anthracene as an example. One 
solution of Maxwell’s equations corresponding 
to a magnetic field parallel to the z axis, with 
strength increasing linearly with the time, 
involves an electric field acting in a direction in 
the xy plane at right angles to the radius vector p 
measured from the z axis and with strength 
proportional to p. The integrated ponderomotive 
force for a segment AB of a conductor is ac- 
cordingly proportional to the area swept out by 
the radius vector p, in this case OAB. Using as 
our unit the area of one-sixth of a carbon hexa- 
gon, the integrated ponderomotive forces for 
anthracene (Fig. 2) are P;=1 (for the segment 
AF), P2=2 (for FGH),and P;=7 (for ABCDEF), 
with values for the other segments given by the 
center of symmetry of the molecule and field. In 
addition we must introduce an auxiliary potential 
at the junction F (relative to A) represented by 
the symbol p. The effective resistances of the 
circuits, arising from the inertia of the electrons, 
are proportional to the lengths of the segments, 
and have the values R:=1, Ro=2, and R;=5. 
We thus obtain the equations 


4:=1- 4, 
2i2=2+>, 
Sis=7—p 


for the currents 7), 72, and 73, together with the 
conservation equation 
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(a) 


(b) 
(c) 


Fic. 4. The aromatic hydrocarbons (a) CosHi2, (b) CssHis, 
and (c) CogH a4. 


12=1; +73. 

The values found on solution of these equations 
are 4;=3/17, i2=24/17, and 73=21/17. These 
correspond to a completed circuit 7; about one 
hexagon and another completed circuit 73; about 
three hexagons, the strength of the equivalent 
magnetic dipole being 66/17. A similar calcula- 
tion for a benzene hexagon gives the value 1; 
the quantity k=66/17 accordingly represents 
the relative strengths of the induced magnetic 
dipoles for the anthracene network of three 
hexagons and the benzene circuit of one hexagon, 
difference in the linear electron density o being 
neglected. 

To obtain the susceptibility difference AK we 
make the following calculation. For naphthalene 
and higher aromatic hydrocarbons the average 
C—C distance’ is 1.41A, corresponding to 
AK = —38.0X10-* for a six-electron hexagon. 
In general we write for AK the expression 


AK = —38.0X10-°X kof, (6) 


in which g is the electron density per C—C bond 
and f a correction factor described below. Values 
of o as well as of k& are included in Table I; for 


7L. Pauling, L. O. Brockway and J. Y. Beach, J. Am. 
Chem. Soc. 57, 2705 (1935). 
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anthracene, for example, o has the value 14/16, 
there being 14 , electrons and 16 C—C segments 
in the network. The factor f is designed to 
correct for the assumption that the electrons 
move along rectilinear paths between the carbon 
atoms. In benzene we assumed the transit from 
atom to atom to be along circular arcs. For the 
condensed ring systems we may correspondingly 
replace the outermost rectilinear segments by 
circular arcs, as shown in Fig. 3 for anthracene. 
This introduces a correction® of 23 percent 
(f=1.23), leading to AK = —159X10-* for an- 
thracene. Table I includes the results of a similar 
treatment of several hydrocarbons. 

In order to predict approximate values of the 
principal susceptibilities K, (=Ke) and K; 
(=K,+AK) of these substances we may adopt 
the following empirical procedure. The molal 
susceptibility of He. is —4.0010-*; we accept 
—2.0010-* for H bonded to C. The atomic 
susceptibility of diamond is —6.00X10-*; we 
accept this for aliphatic carbon. Now in the 
aliphatic carbon atom there are four L electrons 
contributing to the susceptibility in all directions, 
whereas in the aromatic carbon atom only three 
remain to be considered, the fourth being 
involved in the calculation of AK. For aromatic 
carbon, aside from the , electron, we ac- 
cordingly assume the atomic susceptibility 
—4.510-*. The principal susceptibilities of a 
hydrocarbon are thus given by the approximate 
expressions 
K,=Kgz 

=> (2.0 ny+6.0 No ait4.5 no ae x 10-* (7) 
and 
K;=K,+AK, (8) 


8 The error introduced by using the values of relative 
currents for rectilinear segments is small. 











TABLE I. 
o k —AKxXx10® —K,x10® —K;x10 
Benzene, C.He 1 1 49 39 88 
Naphthalene, C\,Hs , 10/11 12/5 104 61 165 
Anthracene, phenanthrene, Cy,Hio 7/8 66/17 159 83 242 
Pyrene, CisHio s 16/19 84/13 237 92 329 
Naphthacene, chrysene, etc., CisHiz 6/7 156/29 214 105 319 
Triphenylene, CisH:2 6/7 60/11 217 105 322 
Perylene, CooHi2 5/6 111/13 314 114 428 
Picene, dinaphthanthracene, etc., CoH, 11/13 227/33 268 127 395 
CuHie 4/5 §2/3 572 132 704 
CsaHig 3/4 2667/29 2770 279 3050 
CocHx 8/11 67333/227 8200 480 8680 
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TABLE II. 
Observed Calculated 
— K,-108 — K,-10° — K;-10° — K,(K2)xX10® —K;-108 

Benzene, CsHe 37 37 91 39 88 
Durene, CioHy 77 86 141 79 128 
Hexamethylbenzene, Ci2His 101 103 164 99 148 
Hexaethylbenzene, CisH3o 165 165 231 159 208 
Biphenyl, CyHio 67 67 175 74 172 
Terphenyl, CisHi4 98 98 260 109 256 
Quaterphenyl, CosHis 129 129 345 144 340 
1,3,5-Triphenylbenzene, C2His 131 131 347 144 340 
Fluorene, Ci3Hio0 73 73 194 80 178 
Naphthalene, CioHs 39 43 187 61 165 
Anthracene, Cy4H10 46 53 273 83 242 
Phenanthrene, CisHi0 74 74 240 83 242 
Pyrene, CisHio 81 81 303 92 329 
Chrysene, CisHie 88 88 306 105 319 

‘ 110 358 127 395 





1,2; 5,6-Dibenzanthracene, CoH 14 110 





in which my, Nc ai, and Nc ar represent the num- 
bers of atoms of these kinds in the molecule.°® 
Values calculated in this way are included in 
Tables I and II. 

The hydrocarbons CosH12, C5sHis, and CosHog 
included in the table have the hexagonal sym- 
metry D¢,, as shown in Fig. 4. It is interesting to 
note that the shape of the molecule has a large 
effect on the susceptibility; AK for the nearly 
circular molecule CosHie is over twice that for 
the long molecules picene, etc. (C22H14), although 
the molecules differ only a little in composition. 
Pyrene (compared with naphthacene) shows the 
same effect. Aside from this, the configuration 
of the molecule is of little significance, the 
susceptibilities for isomers such as anthracene 
and phenanthrene being the same. Branching of 
the chain of benzene residues has a small effect 
(compare naphthacene and triphenylene). 


Comparison with experiment 


In Table II there are given the available 
experimental values of the principal suscepti- 
bilities of aromatic hydrocarbon molecules, from 
the papers of Krishnan and his collabora- 
tors.!>: le. 10 Tt should be mentioned that in some 
cases the values have been obtained from crystal 
susceptibilities in a straightforward manner, with 





* The expressions could probably be made more accurate 
by the introduction of various constitutive correction 
terms (such as corresponding to Pascal’s system). At 
present it does not seem worth while to attempt to refine 
the simple and straightforward treatment given above. 

The values for hexaethylbenzene are those given by 
N. Ganguli, Zeits. f. Krist. 93, 42 (1936), and for 1,2; 
5,6-dibenzanthracene those given by K. S. Krishnan and 
S. Banerjee, Zeits. f. Krist. 91, 173 (1935). 


the use of the results of x-ray investigations of 
molecular orientations, but that in other cases 
the values given depend on various assumptions 
(such as that K;— XK, for phenyl and phenylene 
groups has the value found for benzene). It is 
seen that with certain exceptions mentioned 
below the observed and calculated values agree 
to within about 10 percent, the disagreement 
probably being due in part to experimental error 
and in part to theoretical inaccuracy. 

The largest differences between calculated and 
observed values are shown by naphthalene and 
anthracene. I feel confident that the differences 
are due to errors in the reported experimental 
values, especially since good agreement between 
experiment and theory is found for phenanthrene. 

The principal susceptibilities of graphite, per 
mole of carbon atoms, are reported!! as K,;= —5 
x10-* and Ks= —275X10~. It is to be noted 
that the value of K; is less than for diamond 
(—6X10-*), and agrees with our assumed value 
—4.5X<10~ to within the accuracy of its determi- 
nation. The very large magnetic anisotropy of 
graphite shows that the p, electrons are free to 
move over large areas of the hexagonal network 
of carbon atoms in the crystal. Extrapolating 
from the values of Table I, we find that molecules 
about 40A in diameter would have an anisotropy 
equal to that observed for graphite. It is inter- 
esting that the anisotropy and the average 
susceptibility of graphite are decreased by 
decreasing the size of the crystal grains.” 

1K. S. Krishnan, Nature 133, 174 (1934). 


1225. Paramasivan, Ind. J. Phys. 4, 139 (1929); K. S. 
Krishnan and N. Ganguli, Current Science 3, 472 (1935). 
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The Free Energy of Electron Gas 


A. R. Gorpon, Chemistry Department, University of Toronto 
(Received July 14, 1936) 


The energy and free energy of a semi-degenerate gas obeying the Fermi statistics are com- 
puted as functions of temperature and concentration. The significance of the deviation of the 
free energy from the limiting high temperature value is illustrated by calculating the degree 
of thermal ionization of potassium vapor under conditions of high electron concentration. 





T has long been known that the translational 

partition function for a gas obeying either 
the Bose-Einstein or the Fermi statistics ap- 
proaches asymptotically with increasing temper- 
ature the limiting value given by the Sackur- 
Tetrode equation. In general this approach is to 
all intents complete within a few degrees of 
absolute zero, but in the case of an electron gas 
(owing to the small mass of the particle) the 
deviation from the classical value persists to 
much higher temperatures. Sommerfeld! in a 
well-known paper has derived expressions de- 
fining the translational energy, free energy, 
entropy etc. of a gas obeying the Fermi sta- 
tistics, with special consideration of the case in 
which, owing to high concentration and moderate 
temperature, the deviations from the classical 
values are particularly great (“‘complete de- 
generacy’’); and Mott? has recently computed 
the energy, heat capacity and magnetic suscepti- 
bility for a few points in the ‘‘semi-degenerate”’ 
range. In dealing with the equilibria involved in 
problems of thermal ionization, however, it is 
the free energy that is of primary importance, 
and the present paper deals with the calculation 
of this and related thermodynamic quantities 
for an electron gas when temperature and 
concentration are not such as to permit the 
conventional Sackur-Tetrode expression to be 
used. 

Sommerfeld has shown that for a gas contain- 
ing N=6.06 X10” electrons in V cc at T°K, 


(1) 


(2mkT)*? 2 u}-du 
N=22VG:-———_ f 
“0 


h’ et-e"41 


where G (the spin weight of the electron) =2, 


1 Sommerfeld, Zeits. f. Physik 47, 1 (1928). 
2 Mott, Proc. Camb. Phil. Soc. 32, 108 (1936). 


a=-—F/RT, F is the free energy of one “mole” 
of the gas, m is the mass of the electron, and the 
other symbols have their conventional meanings. 
Eq. (1) can be written 


Q* =e*/Tyj2(a), (2) 


where Q*, the classical expression for the parti- 
tion function, is equal to 


(24 MkT)3!2- VG/N®2h3 


M=Nm=1/1838, and In(a) is defined by 


o u™-du 
In(a)= f —— /r(m+1). (3) 
0 e“t+e-* 


Similarly, the energy of one mole of the gas is 
given by 


E=(3/2)NkT : I3;2(a)/I12(a), (4) 


while pV/NRT =Is3;2(a)/Iijo(a). For large values 
of a, i.e., large T or V, it is easily verified that 
both Ji(a) and I3(@) approach unity, i.e., 
a(= — F/RT)=In Q%, etc. 

Table I gives the values of Ji2(a@) and of 
I3;2(a) for —1<a<6, and Table II gives a as a 
function of InQ* for —0.4<In Q*<8.0. For 
a>S0.3, In(a) was computed by the obvious 
relation 


n(x) = (—1)"-e-"#/ (nA), (5) 


For a<0.3, the range of integration was split; 
for 0<u<0.3, e“ was expanded and the integrand 
reduced to a power series in which could be 
integrated term by term. For 0.3<u<2.0, the 
integrals were evaluated by tabular integration, 
while for 2.0<u< ©, the contribution to Ih(a) 
was computed by the rapidly converging series 
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TABLE I. TABLE II. 
2 7? ub-du 4. 7°? wh-du In Q* a In Q* a | In Q* a 
I, a)=—- [ ’ Iona) =—- | . I] 
Vm do evte* 3Vr Jo ette —0.4 | —0.9169|| +0.7 | +0.5256|) +2.6 | +2.5738 
0.3 | 0.7686 0.8 | 0.6421||) 2.8 | 2.7785 
| ae } | ae 0.2 0.6247 0.9 0.7571 | 3.0 2.9824 
oat wes | Seeshend > usta) | taste) —0.1 | 0.4849/| 1.0 | 0.8706|| 3.5 | 3.4893 
—1.0 | 0.57965 | 0.73659 || +0.8 | 0.87138 | 0.93126 0 0.3488 1.2 1.0940 4.0 | 3.9935 
0.9 | .59976 | .75206 0.9 | .88154 | .93702 +0.1 | 0.2160 14 | 1.3131 4.5 | 4.4961 
0.8 | .61968 | .76704 1.0 | .89104 | .94236 0.2 | —0.0862 1.6 | 1.5288 5.0 | 4.9976 
0.7 | .63935 | .78152 1.2 | .90815 | .95185 0.3 | +0.0408 1.8 1.7417 5.5 | 5.4986 
0.6 | .65871 | .79547 1.4 | .92292 | .95990 0.4 | 0.1652 2.0 | 1.9522 6.0 | 5.9991 
0.5 | .67769 | .80886 1.6 | .93556 | .96669 0.5 | 0.2874 2.2 2.1609 7.0 6.9997 
0.4 | .69624| .82169 || 1.8 | .94631 | .97240 +0.6 | +0.4074 || +2.4 | +2.3680)| +8.0 | +7.9999 
0.3 | .71431 | .83395 2.0 | .95539| 977180 ———————————— a 
0.2 | .73184 | .84562 2.2 | .96303 | .98116 
—0.1 | .74880 | .85670 2.4 | .96942 | .98447 ’ ; : 
0 76515 | .86720 2.6 | .97475 | .98721 without introducing an error greater than 
+0.1 | .78085 | .87712 2.8 | .97919 | .98948 ; : : _ pe 
02 | ‘70588 | ‘98646 30 | ‘98286 | ‘90135 0.01RT; Table II shows that in this case P P 
0.3 | .81022 | .89524 3.5 | .98950| .99472 within one-half of one percent or less. This 
= a. pre re ped ‘oan condition implies that 7 must be greater than 
0.6 | .84901 | .91835 5.0 | .99763 | .99881 1250° for P*=1 atmosphere, and greater than 
+0.7 | 0.86054 | 0.92504 5.5 | .99856 | .99928 ° * 
460 | 0.99913 | 0.99956 20,000° for P* = 1000 wey geen 
For smaller values of In Q*, however, the use 




















2 r@ u)-du oo 
—f ——— =e") (-1)""! 


Vr, etpem® na 


2\! e?" P(2n) 
OSE] 
7 n ni 


where P(@) =(2//7) f ev’ -dy, 
t 





with a similar relation for J3/. 

In dealing with electron concentrations, it is 
convenient to define a fictive pressure P* such 
that P*V = NRT; in terms of this new quantity, 
In Q* = —14.239+1n T*/?/P* where P* is meas- 
ured in atmospheres. It is evident from Table II, 
that for In Q*>3.6, ie., (5/2) log T—log 
P*>7.75, F may be replaced by —RT In Q* 


of such approximations may introduce appreci- 
able errors. To take an extreme case (selected 
from Russell’s tables* of metal ionization in 
stellar atmospheres) consider the ionization of 
potassium at 5000°K under an electron concen- 
tration such that P* = 1000 atmospheres; In Q* is 
then +0.1462, a= — 0.1557, and P= (J3/2, ‘Tiy2) P* 
=1150 atmospheres. From spectroscopic data,‘ 
it follows that if x be the fraction of potassium 
ionized, log x = —4.76; had the Sackur-Tetrode 
approximation to the free energy of the electrons 
been used, the result would have been log x 
= —4.63. 


3 Russell, Astrophys. J. 75, 337 (1932). 

4 Bacher and Goudsmit, Atomic Energy States (1932), 
p. 245. In computing the State Sum for neutral potassium, 
only the low lying states up to and including ?F were taken 
into account. 
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The Exchange of Energy Between Organic Molecules and Solid Surfaces 


Part I. Accommodation Coefficients and Specific Heats of Hydrocarbon Molecules* 
Otto BEECK, Shell Development Company, Emeryville, California 
(Received August 17, 1936) 


A molecular beam method primarily devised for the 
investigation of heterogeneous thermal decomposition of 
hydrocarbon-molecules (to be published later) was used to 
determine the accommodation coefficients of various 
hydrocarbon molecules at 20°C on an ordinarily clean nickel 
surface at 140°C and at pressures in the region of 10-* mm 
Hg. The absolute values of the accommodation coefficients 
were obtained. by relative measurements with the help of 
the following simple relation that was found to exist for 
paraffin hydrocarbons between the accommodation coeffi- 
cient a, the specific heat per mole at constant pressure for 
the ideal gas state C,,, and the molecular weight M, 


a(C,,+R/2)/M=constant =0.352, 


where R is the gas constant per mole. This equation 
supports a further relation which can be set up for the 
mean specific heat between 20° and 140°C for the 
paraffins by using a few reliable C,, values: C,,=0.352M 
+1.80. A combination of both relations gives the simple 
equation for a of the paraffins: 


a=OM/(aM+e), 


where c,=0.352 and c.=1.80+R/2. The accommodation 
coefficient for hydrogen was found to be 0.310, in excellent 


agreement with earlier experimental values. The accom- 
modation coefficients of ethane and ethylene 
specific heats are known from normal frequencies were 
found to be equal within a few tenths of a percent; acety- 
lene was found to be 3 percent lower. Specific heats have 
been calculated for paraffin gases and vapors up to n- 
heptane by a combination of the results of this work and 
the specific heat of ethane from the normal frequencies. 
The specific heats are given to 300°C. Under the assump- 
tion of equal accommodation coefficients for paraffins and 
corresponding olefines based on the accommodation 
coefficient measurements of ethane and ethylene and on 
their specific heats, specific heats are given for propylene 
and n- and isobutylene to 300°C. It was found that the 
absolute difference between the specific heats of a paraffin 
and the corresponding olefine decreases with the chain 
length. The C,, values of the following pairs were found 
to be equal within the experimental error of one percent: 
n- and isobutane, - and isobutylene, methyl acetylene 
and allene. Approximate values for the mean specific 
heats between 20° and 140 C are given for methyl acetylene 
and allene. A discussion is given of the relationship of the 
present work to spectroscopic work in the field. 


whose 





INTRODUCTION 


HIS paper is part of a broader investigation 

that is chiefly concerned with an examina- 
tion of the heterogeneous thermal decomposition 
of organic molecules by a molecular beam 
method. During the course of the investigation 
it was necessary, in order to facilitate the 
analysis of hydrocarbon mixtures and mixtures 
of hydrocarbons with hydrogen at pressures of 
some 10-* mm Hg by the hot wire method, to 
study the relative accommodation coefficients 
of various hydrocarbon molecules. This study 
has disclosed an interesting relationship between 
the accommodation coefficient a, the molar heat 
at constant volume for the ideal gas (zero 
pressure) C,,, and the molecular weight J for 
paraffin hydrocarbons which makes it possible 
to determine the true accommodation coefficient 
of any gas from relative measurements only. 


* Presented in part before the American Physical 
Society at Seattle (Washington), June, 1936. 
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THEORY OF HEAT CONDUCTION AT 
Low PRESSURES 


According to M. Knudsen,! in agreement with 
M. v. Smoluchowski,? that part of heat transfer 
between two surfaces which is due to conduction 
by the gas molecules at low pressures, provided 
the mean free path is large compared with the 
relevant dimensions of the apparatus, may be 
written for monatomic gases 





O- F6n)! 





n@R(T2—T;), (1) 


67)! 


2 
noC,(T2—T;) a” 


where Q is the heat transfer per unit area per 
second, 2 the number of mols per cm’, @ the 
mean molecular velocity, C, = (3/2)R, the specific 
heat per mole at constant volume and (72.—7;) 
the temperature difference of the two surfaces. 
Analogously we may write for polyatomic 


molecules according to Knudsen 


1M. Knudsen, Ann. d. Physik (4) 34, 593 ff (1911). 
2 M. v. Smoluchowski, Ann. d. Physik (4) 35, 990 (1911). 
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where the expression in brackets is equal to the 
sum of the translational part of the specific heat 
plus ? of the internal specific heat. The factor ? 
expresses the fact that, under the assumption of 
independent Maxwell distributions of trans- 
lational and internal energy, a molecule with 
high translational energy contributes more to 
the energy transfer than a slow molecule. 

Eqs. (1) and (2) hold only if there is no 
temperature jump between the surfaces and the 
mean temperature of the gas molecules leaving 
the surface, i.e., if the temperature difference 
(T.’—T,) of the gas before and after leaving the 
surface is equal to the difference of temperatures 
(T:—T,) of the surface and the gas before the 
collision. Otherwise we have to write, more 
generally, 

(T2’ — T)) =a(T2—T)). 


The coefficient a was called by Knudsen the 
accommodation coefficient. Since the present 
measurements, like most of the earlier work, 
have been done with a heated wire the di- 
mensions of which are very small compared 
with the surrounding colder glass tube, the gas 
acquires the temperature of the outer tube and 
we may write Eqs. (1) and (2), by introducing 
the accommodation coefficient a between the 
heated wire and the gas, 





Q= noRa(T2,—T}) (3) 


1 

T 2 

for monatomic gases and 
1 


sal C4) (Tr—T 4 
O= inl +5) “— 





for polyatomic gases. 


EXPERIMENTAL 


As stated in the introduction, the apparatus 
was not primarily designed for the investigation 
of accommodation coefficients. The molecular 
beam method proved, however, to be so appro- 
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priate for the measurement of accommodation 
coefficients at extremely low pressures that it 
may very well be adopted for measurements of 
this type in the future. 

The beam was produced in agreement with 
the specifications outlined in three important 
papers on the molecular beam method by 
O. Stern and collaborators.* The source aperture 
was a short channel 0.50 mm long and 0.528 
mm in diameter. The image orifice, 26.2 mm 
away from the source, was a channel ~1 mm 
long and 0.900 mm in diameter. In line with 
these two orifices was the receiver canal, 62.2 mm 
away from the entrance to the image orifice. 
The canal had a conical shape for better defini- 
tion of the beam although this is not essential 
for the present part of the investigation. The 
entrance orifice had a diameter of 0.795 mm, 
the end of the canal a diameter of 1.60 mm. 
The canal was 11.96 mm long. Its Knudsen 
resistance was 5.19 times larger than the re- 
sistance of an ideal aperture of the diameter of 
the entrance orifice, if the latter resistance is 
taken as 


(27)? 3/2\3 dl 
R,= and the former as R.-(-) = 


’ 
rr? 8i\r 0 8 





. r being the radius and L the length of the canal.4 


The maximum gas pressure in the reservoir 
from which the gas effused was of the order of 
0.1 mm Hg and was measured with an accuracy 
of 0.3 percent by a special manometer previously 
described.’ The two spaces between the three 
apertures were evacuated by two high speed 
high vacuum oil diffusion pumps charged with 
Apiezon oil B. The fore vacuum for these two 
pumps was provided by another oil diffusion 
pump and finally by a “Cenco Megavac’”’ 
rotating oil pump. With a reservoir pressure of 
0.1 mm Hg, a vacuum of about 2-10-' mm Hg 
could be maintained in the first space and better 
than 1-10-° mm Hg in the second space. A 
molecular beam entering an otherwise closed 


3Q. Stern, Zeits. f. Physik 39, 751 (1926); Knauer and 
Stern, ibid., p. 764; Knauer and Stern, ibid. 53, 766 (1929). 
See also the excellent monography on Molecular Rays by 
Ronald G. J. Fraser (The Macmillan Company, 1931), 
The Cambridge series of Physical Chemistry. 

4M. Knudsen, Ann. d. Physik 28, 76 (1909); see also v. 
Smoloct owski, ibid. 33, 1559 (1910). 

5 QO. Beeck, Rev. Sci. Inst. 6, 399 (1935). 
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vessel behind the receiver canal will set up a 
maximum stationary pressure: 


p.. = (pa/mr*)k (5) 


if p is the reservoir pressure, a the area of the 
source aperture, 7 the distance of the source 
aperture from the entrance of the receiver 
channel and k=(R.+R,)/R;, where R., R, are, 
respectively, the Knudsen resistances of the canal, 
and of an ideal aperture of the same cross section 
as the entrance to the canal. Eq. (5) is correct 
only if the source aperture is ideal (infinitesmally 
thin), if source and receiver are at the same tem- 
perature, and if the molecules in the beam suffer 
no collisions in the two spaces between the aper- 
tures. The first and the last assumptions cannot 
be realized in practice. We see from Eq. (5) that, 
under ideal conditions, the pressure p,, is inde- 
pendent of the molecular weight of the gas and 
of the cross section of the canal. A careful and 
critical study of all the errors involved led us to 
the conclusion that a pressure p=0.1 mm Hg in 
the reservoir produced on the average a maxi- 
mum stationary pressure p,,=1-10~° mm Hg in 
the receiver vessel, which is about 10 times lower 
than would be expected under ideal conditions. 
This relationship also varies from gas to gas due 
to several factors, such as the resistance of the 


tube leading from the reservoir to the source: 


orifice, the pumping speed, etc., all depending 
essentially on the mean free path of the gas 
under consideration. The pressure set up in the 
receiver vessel by the molecular beam was 
measured by a highly sensitive Pirani manometer 
of the type used by Knauer and Stern.* An 
exact counterpart of this manometer was placed 
beside the measuring manometer in the same 
constant temperature bath. The opening of this 
manometer led into the same space as the beam 
manometer. Thus it was subject to the same 
temperature and pressure fluctuation as the 
beam manometer. The manometers formed two 
arms in a balanced sensitive Wheatstone Bridge. 
The bridge is unbalanced when the molecular 
beam enters the beam manometer, and the 
observed galvanometer deflection is for small 
pressures proportional to the pressure change. 
The sensitivity of the arrangement used was 
approximately 2-10-' mm Hg pressure change 
per mm deflection of the galvanometer for the 
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monatomic gas argon. The sensitivity was of 
course very much higher for large polyatomic 
molecules. 

Let us, for purposes of simplification, introduce 
the relative sensitivity S which is the galva- 
nometer deflection in cm produced by a reservoir 
pressure of 0.1 mm Hg. The sensitivity is 
specific for a given gas. However, it also depends 
on the reservoir pressure since, as has been 
pointed out above, the conditions are not ideal. 
These difficulties have been overcome by meas- 
uring the sensitivity S as a function of the reser- 
voir pressure for all gases under consideration. 
By extrapolation of these sensitivity curves the 
sensitivity Sp for the reservoir pressure zero may 
be obtained. Except for the possibility of a 
small error due to vapor pressures of the sealing 
materials used in the apparatus, which might 
scatter the various gases in the beam by different 
amounts, the values for Sp should obey Eq. (5), 
i.e., the values for So are independent of the 
pressure p in the reservoir, and the pressure , 
is independent of the molecular weight of the gas. 

Thus the values of Sp may be used directly 
to find the relative accommodation coefficients 
of the gases investigated by combining them 
with Eq. (3) or (4) in which @ is inversely 
proportional to 1/M. One may write, for 
instance, 


(So) gas Qgas(M 4) L (Cy) eas + R/2] 
(So) o4(Mens)!2R 





where the subscript A refers to the monatomic 
gas argon. Then 








Qgas (So) eas( Megas) 12R 
=aq’/= " (6) 
QA (So) a(Ma)3L (Cv) gast R/2 ] 
a’ being the accommodation coefficient relative 
to argon. 


Fig. 1 shows the sensitivity S as a function of 
the reservoir pressure p for all gases and vapors 
investigated (for hydrogen see Fig. 2). It is 
evident that the extrapolation of the sensitivity 
curves to p=0 is permissible only if the indi- 
vidual accommodation coefficients do not change 
in the pressure interval in which the curve has 
been established. Reasonable as this assumption 
may seem because the pressure range is not 
more than twenty-fold and because the pressure 
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Fic. 1. Sensitivity S (i.e., galvanometer deflection per 0.1 
mm Hg gas pressure in the reservoir) plotted against the 
reservoir pressure. 


independence of the accommodation coefficient 
has been well established in higher pressure 
regions for various gases by several investi- 
gators,® the sensitivity curve of ethane shows a 
distinct break at a reservoir pressure of about 
0.1 mm Hg which corresponds to a pressure rise 
of about 1-10-* mm Hg in the manometer vessel, 
caused by the entering beam. To this must be 





*See for instance: Rowley and Bonhoeffer, Zeits. f. 
physik. Chemie B21, 84 (1933). 





added an initial pressure which varies also with 
the reservoir pressure and is probably of the 
same order of magnitude. The sudden change in 
the accommodation coefficient of ethane is 
undoubtedly due to a sudden change of the 
surface conditions (type of adsorption layer) of 
the manometer wire. This change is evidently a 
function of the pressure of the surrounding gas. 

Fig. 2 shows the sensitivity curve of ethane on 
a larger scale together with the sensitivity curves 
of ethylene and acetylene and two points for 
hydrogen. While there is no sharp break in the 
two latter curves, a discontinuity is readily 
observed in the case of ethylene, and acetylene 
also follows the same trend. This discontinuity 
is unique for these three gases. No other gas 
has shown a similar behavior. It seems that the 
pressure region was particularly favorable for 
such an effect with ethane, ethylene and acety- 
lene. The same effect might be expected to occur 
for larger molecules at still lower pressures or 
higher temperatures of the manometer wire since 
the boiling points of these substances are higher. 
Consequently we have used the extrapolation of 
that part of the ethylene curve which belongs to 
the higher reservoir pressures. This extrapolation 
is in excellent agreement with the general 
relations derived from our measurements and 
discussed below. 

The influence of the pumping speed on the 
sensitivity curves is illustrated in the case of 
methane (Fig. 1). Here two sensitivity curves 
have been taken with different pumping speeds; 
curve I with a lower speed than II. Both fall 
together as expected at low reservoir pressures. 

It was found that for paraffins with the 
exception of methane the values for Sp are, as 
shown in Table I, proportional to the square 
root of the molecular weight to a surprising 
degree of accuracy. 

In combining this result with Eq. (4) and in 
setting Q proportional to So, i.e., proportional 
to 4/M, and @ inversely proportional to \/M, 
we arrive at the simple equation: 


a(C,+R/2)/M=const. (7) 


In order to evaluate Eqs. (6) and (7) the 
specific heat C, of the polyatomic gases under 
investigation must be known. To obtain high 
sensitivity, the temperature of the gauge wire 
































was chosen considerably higher than the temper- 
ature of the gauge wall. We therefore have to 
consider the mean specific heat of the gas mole- 
cules between these two temperatures. The mean 
temperature of the wire itself was measured by 
its resistance, the temperature dependence of 
which had been established by a special calibra- 
tion. This calibration was made by placing the 
gauge in constant temperature baths of various 
temperatures and measuring its resistance in a 
high-sensitive bridge with the lowest possible 
current passing through the bridge. By varying 
the current slightly, a curve was established 
whose extrapolation to zero current gave the 
correct resistance for the temperature of the bath. 

Throughout the investigation the gauge wires 
were operated at a mean temperature of 140°C 
and the gauge walls were kept at 20°C. The 
gauges were Pyrex glass tubes about 12 cm long 
and 1 cm internal diameter, each containing a 
nickel ribbon 20 cm long, 50u wide and 3y thick. 
The nickel was shown by conductivity measure- 
ments to be almost pure. Its surface was brilliant 
like that of a platinum wire. The wire was not 
handled, and was cleaned only with alcohol and 
distilled water to remove the flux used in solder- 
ing the wire to its terminals. The wire was not 
heated above the operating temperature of 
140°C. We may call the nickel surface ‘‘ordinarily 
clean.’’ The mean specific heat between 20° and 
140° was calculated from the normal frequencies’ 
for methane, ethane, ethylene and acetylene, by 
the Planck-Einstein function. A value for butane 
at 80°C was taken from a recent investigation by 
Sage and Lacey.® A value for propane by the 


7 The following normal frequencies expressed in wave 
numbers were used: 





CHa CoHe CoHs CoHe 
1304(3x) 600(2x) 730 750(2x) 
1520(2x) 729(2x) 950 827 (2x) 
2915 975 1097 993 
3020(3x) 3277 1110 1370 

3370 1160 1460(3x) 
1342 1480(2x) 
1444 2890(2x) 
1623 2899(2x) 
2990 2955 
3019 3020 
3110 
3240 


For references on this subject see: A. Eucken and A. 
Bertram, Zeits. f. physik. Chemie B13, 361 (1936), and 
E. Teller and B. Topley, J. Chem. Soc. London, 885 (1935). 

8B. H. Sage and W. N. Lacey, Ind. Eng. Chem. 27, 
1484 (1935). The authors give a straight line between 60° 
and 220°F. Since in case of ethane the mean specific heat 
between 20° and 140°C is very close to the true specific 
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Fic. 2. Sensitivity S (i.e., galvanometer deflection per 
0.1 mm Hg gas pressure in the reservoir) plotted against 
the reservoir pressure. 


same authors was discarded because of its 
uncertainty. 

Remembering that Soddy and Berry,’ Knud- 
sen! and others have shown previously that on 
ordinarily clean surfaces the accommodation 
coefficient rapidly approaches unity with in- 
creasing molecular weight, we see from Eq. (7) 
that the increase of C, must be proportional to 
M for high molecular weight paraffins. This 
proportionality, which is demanded by Eq. (7) 
is in fact approximately true for the specific 
heat of the three available lower paraffins, 
methane, ethane and butane and has also been 
observed by Lewis and McAdams." It is needless 
to say that since we are concerned with specific 
heats at extremely low pressures, we are dealing 
essentially with C,, values, i.e., with specific 
heats of the ideal gases. Calculation from normal 











TABLE I. 
So VM So/v¥ M 
CH, 13.9 4.00 3.48 
C:H¢ 19.5 5.48 3.56 
C;Hs 23.5 6.64 3.54 
CiHio 27.1 7.63 3.55 
CrHie 35.6 10.01 3.56 








heat at 80°C, we have chosen for butane the true specific 
heat value at 80°C for reference in order not to depend on 
the slope of the line which fits the experimental points far 
better in the neighborhood of 80°C than at lower tem- 
peratures. e 

°F. Soddy and A. J. Berry, Proc. Roy. Soc. A83, 254 
(1910); ibid. A84, 576 (1911). 

10 W. K. Lewis and W. H. McAdams, Chem. and Met. 
Eng. 36, 336 (1929). See also p. 689 of this paper. 
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frequencies give C,, values directly. The C.,, 
value for butane was calculated from Sage and 
Lacey’s measurements by means of Berthelot’s 
equation of state. The almost perfect propor- 
tionality of the increase with M of those three 
values for C,, is at least a good support for the 
deductions made from Eq. (7). We feel, however, 
that there are good reasons for believing the 
methane value to be slightly lower (see discussion 
below). The equation for C,, was set up, there- 
fore, by use of the ethane and butane values 
only. Thus we arrive at an equation representing 
the mean specific heat of the paraffins, between 
20° and 140°C, in terms of the molecular weight. 


C,, =0.352M+1.80. (8) 


Next a’ has been calculated according to Eq. (6) 
for all cases in which the specific heat is known 
or can be calculated from Eq. (8). The values 
for a’ are given in Table IT. 

By introducing the values for a’ and for C,, 
into Eq. (7) we find for paraffins the relation 


(C.,+R/2)a’ =0.388M. (9) 
and by combining Eqs. (8) and (9) 
a’ =0.388M/(0.352M+1.80+R/2). 
For large M, we obtain 


lim a" _ 9 388/0.352 


and since this value must be identical with the 
true accommodation coefficient a=1, we obtain 
a=(0.352/0.388)a’, 














0.352 0.3881 
ae ) 
0.388 \0.352M+1.80+R/2 
0.352M 
0.352M+1.80+R/2 
TABLE II. 

Gas a’ _ “ ; 
Argon 1 0.90; 
Hydrogen 0.342 0.31 
Methane 0.737 0.665 
Ethane 0.87; 0.79 
Ethylene 0.87 0.79; 
Acetylene 0.84, 0.765 
Propane 0.93, 0.84; 
n-Butane ” 

Isobutane } seth — 


n-Heptane 1.02; 0.925 
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Fic. 3. The accommodation coefficients of the paraffins and 
of hydrogen plotted against the molecular weight. 


The values for a are given in Table II. In Fig. 3, 
a for the paraffins is plotted against the molecular 
weight. a for hydrogen is added for comparison. 


SpEcIFIC HEATS 


In examining Table II we find equal accom- 
modation coefficients for iso- and normal butane 
and for ethane and ethylene. The former was 
derived from the fact that the sensitivity curve 
(Fig. 1) is the same for both gases. Strictly 
speaking, this only means that in both cases 
the product a-(C,,+R/2)=constant. If, there- 
fore, we assume a to be the same for iso- and 
n-butane, we must conclude that C,, is also the 
same. In view of the fact that ethane and 
ethylene have been shown to have the same 
accommodation coefficient, equality of accom- 
modation coefficients of the isomers of butane 
seems almost certain. Furthermore, the assump- 
tion seems justified that corresponding paraffins 
and olefines have the same accommodation 
coefficient. Again, since isobutylene and normal 
butylenes" lie on the same sensitivity curve, 
their accommodation coefficients and specific 
heats are presumably equal. The same sensitivity 
curve has been observed for allene and for 
methyl acetylene. Their absolute value, however, 


11 A mixture of 25 percent a-butylene and 75 percent 
8-butylene was used. 
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may not be obtained with the same certainty as 
the absolute values of the olefines, because we do 
not feel justified in setting their a equal to that 
of propane and propylene due to the fact that 
the accommodation coefficient of acetylene is 
about 3 percent lower than that of ethane and 
ethylene. The discrepancy will probably be less 
in the case of higher molecular weights and an 
approximate value of the specific heat of methyl 
acetylene and allene may be obtained by setting 
their accommodation coefficient equal to that of 
propane and propylene. In Table III are given 
the mean specific heats between 20° and 140°C 
at constant volume for the ideal gas state of all 
gases and vapors investigated. 

This, however, is not the only information on 
specific heats that may be gained from these 
observations. In case of the paraffins, Eqs. (7) 
and (8) have been established for mean specific 
heats between 20° and 140°C. It follows therefore 
that the increase of the true specific heats at any 
temperature in this temperature interval must 
also be proportional to the molecular weight. 
Knowing accurately from spectroscopic data, in 
the case of ethane the relation of specific heat to 
temperature we are able to calculate true specific 
heats of all other paraffins at any temperature. 
Assuming that the laws established in the tem- 
perature region between 20° and 140°C will hold 
likewise at higher temperatures, we have calcu- 
lated the specific heats of the paraffins up to 
300°C, using ethane as reference. Furthermore, 
the specific heat of ethylene being known from 
spectroscopic data, we find by comparison of the 
curves for ethane and ethylene, that their ratio 
stays almost constant. At 300°C the ratio differs 
only 2 percent from its value at 0°C. Such a dis- 
crepancy necessarily will be smaller for larger 
paraffin molecules and their respective olefines. 

Using Eq. (5) we arrive at the specific heat of 
the paraffins at any temperature by means of the 
relation 


C.= (Ci — 1.80)/(C’20 to 1407 1.80)) 


X (C20 to 140— 1.80) +1.80, (10) 


where C; is the specific heat per mole of the gas 
considered, at zero pressure and at the tempera- 
ture ¢, and C2o_149 is the mean specific heat be- 
tween 20° and 140°C of the same gas; C’ refers to 
ethane, the reference gas, and 1.80 is an experi- 
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mental constant. The values for the olefines are 
obtained by using the experimental ratio of 
paraffins and their respective olefines found for 
the mean specific heats between 20° and 140°C. 

In Table IV are given the C,, values for 
C,,Hen42 up to m=7, and also ethylene, propylene 
and acetylene. Further values for paraffins may 
be calculated from Eqs. (8) and (10) using the 
corresponding ethane value as standard. 

Table V gives the specific heat C, at constant 
pressure (one atmosphere) calculated with the 
help of Berthelot’s equation of state. This equa- 
tion is known to hold well for hydrocarbon vapors 
at pressures sufficiently below the critical pres- 
sure, which is true for all vapors under considera- 
tion. From Berthelot’s equation we derive : 


C»p=C.,+R+(81/32)R(T./T)*p/p., 


where 7. and , are the critical temperature and 
the critical pressure. 


DISCUSSION 
Paraffins 


A discussion of the theoretical aspects of the 
law governing the accommodation coefficients of 
the paraffins will be postponed until data for 
various gas temperatures and various surface 
temperatures have been completed. 


Hydrogen 


The value for the accommodation coefficient a 
of hydrogen is in excellent agreement with values 
recently obtained by other investigators on vari- 
ous surfaces, especially with the value obtained 











TABLE III. 
Cory 
mean 20° to 
Gas 140°C Remarks 
CH, 7.17 from normal frequencies 
CH 7.25 “ rae (Eqs. (7) and 
8 
C2He 12.35 from normal frequencies (refer- 
ence value) 
C3Hs 17.25 this investigation (Eq. (8)) 
n-C4Hi0 22.2 Sage and Lacey® (reference 
value) ; 
iso-C4Hio 22.2 assuming equal accommodation 
coefficient for n- and iso-CsH1o 
n-C7Hie 37.0 this investigation (Eq. (8)) 
CHa 9.55 from normal frequencies ; 
C3He 15.09 assuming equal accommodation 
coefficient for CsHs and CsHs 
n- ‘ae 20.5 assuming equal accommodation 
4s0- coefficient for m-CsHio and 1s0- 
and n-C4Hs 
C2He 9.50 from normal frequencies ; 
allene \ouHy 13.5 assuming equal accommodation 


coefficient for CsHs and Csli« 


methyl acetylene / 
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TABLE IV. 








C,, in cal./mole °C 
100 


200 300 
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TABLE V. 








C, (at 1 atmos.) in cal./mole °C 
0 100 20 





8.76 
16.06 
12.32 
10.89 
22.70 
19.82 
29.34 
27.09 
35.98 
42.62 
49.27 


10.19 
18.78 
14.32 
11.73 
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by Rowley and Evans” on an ordinarily clean 
iron surface at the same temperature as was used 
in this investigation. Their value is 0.31. Rowley 
and Evans by comparing their results with others 
came to the conclusions that they had measured 
the accommodation coefficient on an iron surface 
covered with adsorbed hydrogen. This conclusion 
was based on the fact that approximately the 
same value was found by other investigators on 
various materials. It seems most likely that at a 
surface temperature of 140°C a gas layer was 
present on the nickel used in this investigation 
or on the iron used by Rowley and Evans. Since 
this layer should not be affected by the hydrogen 
pressure at a moderately high temperature of the 
surface, it is not astonishing that the value found 
in this investigation with hydrogen pressures of 
some 10-* mm Hg is in agreement with the value 
of Rowley and Evans at approximately 0.025 
mm Hg. Earlier work of Rowley and Bonhoeffer, 
for instance, has shown the accommodation co- 
efficient of hydrogen on platinum to be inde- 
pendent of the pressure between 0.03 and 0.3 
mm Hg. 

It must be emphasized that a direct measure- 
ment of the accommodation coefficient at pres- 
sures of 10-®° mm Hg would lead to great experi- 
mental difficulties since the energy radiated from 
the surface at 140°C is about 10‘ times larger than 
that resulting from the collision of the gas mole- 
cules. By means of the relations found for the 
paraffin series, the true accommodation coefficient 
for hydrogen or any other gas at 10-* mm Hg was 
found by relative measurement only. 


” H. H. Rowley and W. V. Evans, J. Am. Chem. Soc. 
57, 2059 (1935). 


Specific heats and spectroscopy 

The observed proportionality between the in- 
crease of C,, and M for the paraffins suggests at 
once a simple relationship of infrared and Raman 
spectra in the paraffin series. A glance at the 
graphic presentation of the Raman lines from 
n-butane to n-dodecane reveals indeed striking 
similarity and comparative simplicity. This 
agreement is not only true for frequencies as- 
signed to C—H bonds and C—H deformations 
which lie above 1100 cm, but is surprisingly 
good in the range of the chain oscillations proper. 
Propane and ethane show a much simpler Raman 
spectrum. The intimate relation between ethane, 
propane and -butane, however, may be readily 
shown by an empirical summation rule of Kohl- 
rausch and Barnés,'* as applied to the C—C fre- 
quencies. This rule is applicable also to cyclo- 
propane, isobutane and neopentane (tetra- 
methylmethane). The rule is given by the equa- 


tion: - 
v=((1/(s—1))dv2)}, 


where y is a frequency characteristic for C—C 
chains. In making the summation, frequencies of 
2 or 3 fold degeneracy must be taken 2 or 3 
times. s—1 is the number of similar valency bonds. 
Kohlrausch and Képpl found that the gases 
named above satisfy the rule with an error of 
one percent. 

A few remarks are necessary in respect to 
methane. In spite of the fact that Eqs. (8) and 
(9) seem to hold for methane within a few per- 
cent, there are good reasons for believing that the 
methane value is slightly lower. 


13K, W. F. Kohlrausch and D. Barnés, Ann. Soc. 


Espafi. Fis. Quim. 30, 733 see also: K. W. F. Kohlrausch 
and F. Képpl, Zeits. f. physik. Chemie B26, 209 (1934), 
especially p. 218 ff. 

























































688 


At 0°C the contribution of the oscillatory de- 
grees of freedom to the specific heat is very small 
in the case of methane, and although this contri- 
bution still decreases towards lower tempera- 
tures, it is almost constant in the close neighbor- 
hood of 6R/2. Towards lower temperatures the 
methane value is therefore too high and Eq. (8) 
must break down completely in the case of 
methane. At higher temperatures the value for 
the specific heat of methane which satisfies Eq. 
(8) is slightly greater than the value calculated 
from normal frequencies. In Fig. 4 are given 
curves for ethane and methane as calculated from 
normal frequencies and the curve for methane 
given by Eq. (8). It is readily seen that below 
0°C. Eq. (8) breaks down entirely. As to the 
discrepancy at higher temperatures the following 
explanation may be offered without going into 
the details of this complex problem. Let us con- 
sider only the C—H frequencies and the C—H 
deformation frequencies. If the specific heat of 
the paraffins is proportional to M, we may as- 
sume that the contribution of the C—H fre- 
quencies and the C—H deformation frequencies 
is additive within the limit of error of our in- 
vestigation. This can only mean that these 
oscillations are well localized and practically 
independent of the chain length, or in other 
words that we have a constant binding force. The 
group of C—H valence frequencies and deforma- 
tion frequencies in the Raman spectra is indeed 
remarkably constant, especially from butane 
upwards. 

There are two frequency regions of C—H de- 
formation frequencies for paraffin molecules from 
ethane upwards, one region at about 750 cm 
and the other about 1400 cm-. The frequency 
750 cm seems definitely to be associated with a 
bond angle that only occurs when two or more 
carbon atoms are present, probably the angle 
C—C-H. 

Any attempt, therefore, to extrapolate to 
methane the measurements with higher paraffins 
would include the consideration of the deforma- 
tion frequency at about 750 cm-!. The number 
of frequencies depending on the bond angle 
C—C-—H decreases according to: 10 for butane, 
7 for propane, 4 for ethane, which would give the 
extrapolated value 1 for methane. The true value 
for methane is, however, zero. The contribution 
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Fic. 4. The specific heat in cal. per mole and °C for the 
ideal gas state calculated for ethane and methane (curve a) 
from the normal frequencies and the values for methane 
(curve 6) calculated from Eq. (8). 


of this excess frequency to the specific heat is in 
rough agreement with the discrepancy between 
the values from normal frequencies and the value 
obtained by Eq. (8), the latter value being too 
high. This shows clearly that Eq. (8) must not be 
expected to hold for methane. 

It is noteworthy that the absolute value of the 
difference between the specific heat of a paraffin 
and the respective olefine becomes smaller with 
increasing molecular weight. This indicates a 
definite effect of the chain length on the double 
bond. The only (and an apparently safe assump- 
tion), in obtaining the C,, value for the olefines 
is that of equality of the accommodation co- 
efficients for paraffins and respective olefines; no 
other errors are possible. A similarly pronounced 
effect cannot be observed in Raman and infrared 
spectra, but this might be due to the funda- 
mentally different method of excitation. 

Special mention should also be made to the 
equality of C,, for the following pairs: iso and 
normal butane, iso and n-butylene, methy! 
acetylene and allene. In the latter case two 
double bonds are substituted for one single and 
one triple bond without altering the specific heat 
within the limits of the experimental error of 
one percent. 

While it is realized that this discussion is in- 
complete from the spectroscopic standpoint, it is, 
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nevertheless, hoped that the specific heats and 
their simple relationships now presented will 
contribute considerably to the proper assignment 
of fundamental frequencies to the higher paraf- 
fins and by so doing will provide the foundation 
for a correct direct calculation from the normal 
frequencies of the specific heat of many higher 
hydrocarbons at higher temperatures. 


Sources of error 


Apart from the fact that the limit of error of the 
molecular beam measurements has been found to 
be one percent, thus establishing the proportion- 
ality of a(C,,+.R/2) with M to within one per- 
cent, the actual equation for C,,, Eq. (8), has 
been set up by the use of the specific heat value of 
ethane (from normal frequencies)’ and n-butane 
(Sage and Lacy).§ Inasmuch as the value of Sage 
and Lacey may be slightly in error, all other 
derived values may show about the same relative 
error. We have, however, confidence in their 
value for n-butane since it was obtained by an 
adiabatic expansion method which Eucken and 
collaborators!’ have used previously with much 
success. We furthermore call attention to specific 
heat values for the paraffins obtained by Lewis 
and McAdams" by a flow method not designated 
in detail. Inasmuch as the flow method is known 
to introduce systematic errors easily, it is not 
astonishing that they find C, values at 50°C for 
methane and ethane which are about 1.2 calories 

't A. Eucken and K. von Liide, Zeits. f. physik. Chemie 


assy (1929); A. Eucken and A. Parts, ibid. B20, 184 
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larger than those from spectroscopic data. How- 
ever, they find in agreement with this work an 
increase in specific heat of the paraffins from 
methane to butane, which is, as they prefer to 
state, proportional to the number of carbon 
atoms. If their error is indeed systematic, as may 
be expected, their values should give the same 
slope as the values obtained from this investiga- 
tion. This is actually the case within one percent 
accuracy, which speaks for the correctness of 
Sage and Lacey’s value for n-butane. 

As may be seen above, the accuracy of the 
values for the olefines depends on the accuracy of 
the values for the respective paraffins and at 
higher temperatures it depends also on the as- 
sumption that at any temperature the ratio of the 
specific heats of an olefine and the respective 
paraffin is constant, an assumption which is 
practically fulfilled in the case of ethane and 
ethylene as has been stated above in more detail. 
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1. INTRODUCTION 


[\ the thermodynamics of surfaces there is no 

theorem more important than that known as 
“Gibbs’ adsorption theorem,”’ which determines 
how the surface tension of a solution varies when 
the composition is varied at constant tempera- 
ture. Its derivation by Gibbs has all the logical 
conviction and elegance characteristic of Gibbs. 


But the abstract character of his proof makes it 
difficult to follow by anyone not thoroughly 
conversant with Gibbs’ methods. Owing to the 
numerous applications of the theorem to all 
branches of surface chemistry and particularly 
colloid chemistry it is very desirable to have a 
more elementary proof. 

Various attempts in this direction have already 
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been made. Of these the place of honor must with- 
out question be given to a proof by Koenig and 
Swain.' This proof is accurate, clear and con- 
vincing to anyone who is able to follow it. It 
however postulates a familiarity with the thermo- 
dynamic functions known as free energy and 
chemical potentials. It therefore seems to the 
author for some purposes insufficiently ele- 
mentary. The proof is in fact almost exactly 
Gibbs’ proof translated into less abstract and 
more comprehensible terms. 

Of a very different class is the alleged proof of 
Milner,’ which is reproduced in numerous ele- 
mentary textbooks, and a somewhat similar 
proof of Harlow and Willows.’ In neither of these 
treatments is the formula derived an accurate 
expression of Gibbs’ adsorption theorem in its 
most general form. Both these treatments assume 
firstly that the solution contains only two com- 
ponents, secondly that the solution is ideal and 
thirdly that one of the components, regarded as 
solute, is present in very small quantity com- 
pared with the other, regarded as the solvent. 
Actually many of the systems to which the 
theorem is most widely applied have none of 
these properties. These treatments also have the 
drawback of relating the surface tension to the 
osmotic pressure instead of to the partial vapor 
pressures, a drawback common to numerous 
thermodynamic discussions of a past generation. 
We shall return to a more detailed criticism of 
Milner’s formulae in a later section. 

There is also a proof by Porter,* which does 
avoid the introduction of the quite irrevelant 
osmotic pressure, but there is a step in the argu- 
ment which we are unable to follow.’ Moreover 
the proof involves unnecessary approximations, 
in particular an inexact form of Raoult’s law is 
assumed. Finally the proof seems in no respect 
easier to follow than that of Gibbs. 

The proof which we shall give makes no refer- 
ence to any thermodynamic potential functions 
and yet is free from restriction to ideal solutions 
or solutions of only two components. We con- 
sider that an elementary proof of such an im- 


1 Koenig and Swain, J. Chem. Phys. 1, 723 (1933). 

2 Milner, Phil. Mag. 13, 96 (1907). 

* Harlow and Willows, Trans. Faraday Soc. 11, 53 (1915). 

‘ Porter, Trans. Faraday Soc. 11, 51 (1915). 

® The step in question is worded: ‘‘since dN and dn must 
be in the same proportion as N and n.” 
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portant thermodynamic theorem is highly il- 
luminating and instructive, provided it is logical 
and accurate. Alleged elementary proofs in which 
simplicity is obtained at the expense of clear and 
exact reasoning are to be deprecated. As there is 
so much loose reasoning prevalent in elementary 
thermodynamics due to confusion between the 
total reversible work of a process and the net or 
useful reversible work (in the terminology of 
Lewis and Randall) we shall commence with a 
brief discussion of these two expressions. We 
shall next discuss an important theorem equiva- 
lent to the formula called Gibbs-Duhem-Mar- 
gules. An illuminating elementary proof of this 
theorem has been given by Bjerrum.® The Gibbs- 
Duhem-Margules relation has nothing to do with 
surfaces but its form is very similar to that of the 
Gibbs adsorption formula. It is therefore perhaps 
not altogether surprising that we are able to de- 
rive the latter from the former in a brief and 
simple manner. 


2. REVERSIBLE WORK AND NET (OR USEFUL) 
REVERSIBLE WORK 


If we consider two completely defined states I 
and II of a system, the two states being of the 
same temperature, then the work done by the 
system in passing isothermally from the state | 
to the state II has a maximum value if the path 
of the change is a reversible one; this maximum 
value is the same for all reversible paths and is 
called the (total) reversible work. It will be 
denoted by w. 

In the particular case that in states I and II the 
whole system is at the same external pressure P, 
then one part of the reversible work may be 
regarded as done against this external pressure. 


If we denote by V! and V™ the initial and final 


volumes of the system, this part of the work is 
evidently P(V!'— V‘). If we subtract this quan- 
tity from the (total) reversible work, we obtain 
the residue w— P(V™— V!) which is available for 
doing electrical work or other useful work. It has 
therefore been called by Lewis and Randall’ the 
net (reversible) work of the process. It will be 
denoted by w’. In the more general case that the 
initial and final states are not at the same pressure 


6 Bjerrum, Zeits. f. physik. Chemie 104, 410 (1923). | 
7 Lewis and Randall, Thermodynamics (McGraw-Hill 
Book Company, 1923), p. 157. 
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or that there are phases at different pressures, the 
net (reversible) work is defined by 


w =w—DPUyN+y Pry, 


the summation sign >> extending over the 
various phases. 

Let us now compare the values of w and w’ for 
two important processes. The first is the evapora- 
tion of one mole of a component from a con- 
densed phase (crystal, single liquid, or solution) 
through a membrane permeable only to the vapor 
of this particular component to the vapor phase 
containing only this component. If P denotes the 
pressure on the condensed phase, p that of the 
vapor, v the partial molar volume of the particu- 
lar component in the condensed phase and u its 
molar volume in the vapor, then evidently 


(2.1) 


(2.2) 
(2.3) 


w= pu—P2, 
w’ =0. 


The second, and more interesting, process that we 
shall consider is the transfer of one mole of a 
component from one condensed phase denoted 
by the superscript a to another condensed phase 
denoted by the superscript 8. The simplest 
method of reversible transfer is the familiar 
distillation in three stages: 


(1) evaporate from a through membrane permeable only to 
vapor 
w= pru*— P%y*, 


(2.4) 

(2) expand from p* the partial vapor pressure over a to p® 
the partial vapor pressure over 8 
p=pB 

pdu, 


p=p® 


(2.5) 


(3) condense into 6 through membrane permeable only to 
vapor 


w= — Pub + PHP, (2.6) 


For the complete transfer, we have by addition 
p= p8 
w= prur — Pry*+ [ pdu — piu + Pby8 (2.7) 
e p=p= 
or integrating by parts 
pB 
a f udp—Pv*+P%P. (2.8) 
p® 


Hence, according to the definition (2.1) of w’, 
we obtain the extremely simple yet exact 
formula 
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pb 
w=—f udp. (2.9) 


Pp 


If the vapor behaves as a perfect gas, then 
u= RT p, (2.10) 


where J is the absolute temperature and R the 
gas constant. Hence by substituting into (2.9), 
we obtain 


»B 
w'=— [" (RT/p)dp, =RT log (p"/p4). (2.11) 


If on the other hand the vapor does not behave as 
a perfect gas, then following G. N. Lewis* we 
define the fugacity p* by the relation 


f udp=RT log (p*/p"), (2.12) 


p° 


where p° is some pressure sufficiently low for the 
vapor to behave as a perfect gas. Substituting 
into (2.9) we obtain 

w' =RT log (p**/p**). 


For the sake of simplicity we shall write for the 
general case 


(2.13) 


w’ = RT log (p*/p?) (2.14) 


on the understanding that if the vapor is not a 
perfect gas then p*, p*® denote the fugacities and 
not the partial vapor pressures. With this reserva- 
tion we shall usually refer to p*, p* as the partial 
vapor pressures, though strictly they denote 
fugacities if the vapors do not behave as perfect 
gases. 

Suppose that two phases a and £6 at the same 
temperature, but not necessarily at the same 
external pressure, are in equilibrium with respect 
to a particular component 7. Then there is no 
tendency for this component to distil over from 
the one phase to the other and so its partial vapor 
pressure must have the same value over both 
phases : 


pt=pe. (2.15) 


By substitution of (2.15) into (2.14), we obtain 

w’ =0. (2.16) 
That is to say, the condition for the equilibrium 
with respect to a given component of two phases 


8 Lewis, Proc. Am. Acad. 37, 49 (1901); Zeits. f. physik. 
Chemie 38, 205 (1901). 
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at the same temperature, but not necessarily at 
the same pressure, is that the net work of transfer 
of the component from the one phase to the other 
be zero. It is with respect to this condition that 
the net reversible work w’ is more important than 
the reversible work w. In particular the net work 
of transfer of a component from the interior of a 
phase to its surface is zero, provided the surface 
is in equilibrium with the interior. 


3. GiBpBs-DUHEM-MARGULES THEOREM 


This theorem can be expressed in a very simple 
form by making use of the conception of net 
work. Consider two solutions at the same temper- 
ature and differing only infinitesimally in compo- 
sition; then the net work of transfer of a finite 
quantity, say one mole, of any single component 
will evidently be infinitesimal. But the net work 
of transfer of a finite quantity of either solution 
into the other solution is an infinitesimal of a 
higher order or effectively zero. This last state- 
ment is equivalent to the Gibbs-Duhem-Mar- 
gules relation. It can be expressed mathemati- 
cally as follows: consider a mixture composed of 
n, moles of component 1, 72 moles of component 2 
and so on, and imagine its composition to be 
varied at constant temperature and pressure, 
then 


mRTd log pi1+n2RTd log po+-+--=0 (3.1) 
or dividing throughout by RT 


nid log pitned log poet +++ =). (3.2) 


As an elementary proof of this important 
theorem has been given by Bjerrum® we shall 
merely illustrate it by a simple example. In a 
mixture of 0.40 mole of alcohol to 0.60 mole of 
water the partial vapor pressure of the alcohol is 
34.2 mm and that of the water is 18.35 mm. Ina 
mixture of 0.41 mole of alcohol to 0.59 mole of 
water the partial vapor pressure of the alcohol is 
34.6 mm and that of the water is 18.25 mm. The 
net work of transfer of one mole of alcohol 
from the first mixture to the second is RT log 
(34.2/34.6) and the net work of transfer of one 
mole of water isRT log (18.35/18.25). These quan- 
tities are of opposite sign and both are of course 
small compared to RT because there is only a 
small difference in the composition of the two 


® Bjerrum, Zeits. f. physik. Chemie 104, 410 (1923). 
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mixtures. But the net work of transfer of one 
mole of the first mixture, made up of 0.40 mole 
alcohol+0.60 mole water, into the second mix- 
ture will be 


0.40RT log (34.2/34.6) +0.60RT log (18.35/18.2) 
= RT{0.40 log (34.2/34.6) +0.60 log (18.35/18.2)} 
= RT(—0.0046+0.0049), 


which is numerically very much smaller than 


‘either of the two terms. 


The special form of the theorem which we shall 
use is the following: if to a large, effectively 
infinite, volume of a phase we add a finite 
quantity of another phase of the same tempera- 
ture, the same pressure and infinitesimally differ- 
ent composition, the net work is zero (strictly an 
infinitesimal of a higher order). As we shall apply 
this theorem to a system containing free surfaces, 
which may not be ignored, it is necessary to be 
clear as to the conditions of validity of the 
theorem in the presence of such surfaces. These 
conditions are simple. We may apply the theorem 
in the form stated above even in the presence of 
free surfaces, provided these surfaces are in 
equilibrium with the volume phases which bound 
them and that the process of mixing does not 
create or destroy any free surfaces. 


4. GiBBs’ ADSORPTION THEOREM 


After these preliminaries the actual proof of 
the adsorption theorem is extraordinarily simple. 
In the diagram KABCDEFG/J is the cross section 
of the wall of a container. HC is a partition 
separating the small volume ABCH, from the 
much larger, effectively infinite, volume HDEF. 
The partition HC is hinged at H and can be 
rotated slightly so as to allow the contents of the 
two compartments ABCH and HDEF to mix. 
ABCH is filled with a solution of composition 
such that the components 1, 2, --- have partial 
vapor pressures (strictly fugacities) ~1, po, --:- 
HDEF is filled with a solution of slightly different 
composition such that the components 1, 2, --: 
have partial vapor pressures (strictly fugacities) 
pitdpi, p2tdpe, ---. AH isa cover which can be 
slid from the position AH to the position HG, 
each supposed to be of area unity. This lid is 
imagined to have the following properties: 


(1) The composition of the solution in contact with the 
lid remains absolutely uniform right up the surface of the lid. 
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(2) The interfacial tension between the lid and the solu- 
tion in contact with it is independent of the composition of 
the solution. 


Actually according to the Gibbs adsorption 
theorem, which it is our object to prove, these 
two properties are not independent but thermo- 
dynamically equivalent to each other. 

Suppose that we rotate the partition HC so 
that the two solutions mix. Then, since the mix- 
ing takes place without the formation or dis- 
appearance of any free surfaces, we may apply 
the Gibbs-Duhem-Margules relation and so the 
net reversible work of this process is zero. Now 
let us imagine the system in its original state and 
let us bring it to the same final state as before by 
a different path. Let us commence by gradually 
sliding the lid across from the position AH to the 
position HG. As we do so we gradually build up a 
free surface A// and destroy the free surface HG. 
In so doing we wish to keep unaltered the com- 
positions of the interiors of the mixtures. This we 
can do by distilling over from the right to the left 
such amounts (positive or negative) of each com- 
ponent as are present in excess at the free surface 
on the left. If we do this then the surface tension 
of the two free surfaces, the one expanding, the 
other contracting will each remain constant 
throughout. After the lid has been moved to HG 
we rotate the partition HC so that the solutions 
mix. Since the mixing takes place again without 
formation or disappearance of any free surfaces, 
we may again apply the Gibbs-Duhem-Margules 
theorem and therefore the net work of mixing is 
again zero. But after mixing we can evidently 
slide the lid back to its old position AH without 
any work. But we have now reached the same 
final position as by the first and simpler procedure. 
Hence the net work is zero for the series of proc- 
esses: sliding the lid from AH to HG, mixing 
solutions, sliding lid back to AH. But we have 
already pointed out that the net reversible work 
of the last two processes is zero; hence it must 
also be zero for the first process. All that remains 
to be done is to write down the net work of the 
process of moving the lid across from AH to HG 
while at the same time distilling over such 
amounts of the components as to keep the com- 
position of each interior and therefore each 
surface tension unaltered and then equate this 
net work to zero. 
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-++ denote the number of moles of 
- in excess (positive or 


Let 1), Ts, 
the components 1, 2, 
negative) per unit area of free surface over what 
would be present if the composition remained 
uniform right up the surface. Then T;, Ts, ---, 
are the quantities of the components 1, 2, --- 
that are distilled over from the right to the left. 
The net reversible work of these distillations is 


T,RTd log pit TeRTd log pot+--+-. (4.1) 


If y denotes the surface tension of the mixture on 
the left and y+dy that of the mixture on the 
right, then the work done by the surface tension 
when the lid is moved from AH to HG is evi- 
dently dy. Adding this to the net work of the 
distillations and equating to zero, we obtain 


dy+TRTd log 1+ T2RTd log po+---=0 (4.2) 
or 
—dy=T,RTd log p1+- T2RTd log po+- - - =0. (4.3) 


This is the accurate and general form of the 
Gibbs adsorption theorem. 


5. PuHysicAL MEANING OF IT AND INVARIANT 
FORM OF ADSORPTION FORMULA 


Actually formula (4.3) is applicable to the 
surface separating any two phases, but we wish 
to confine our attention to surfaces separating a 
liquid mixture from a vapor phase. We may then 
assume that every component of the solution has 
a volume concentration in the vapor phase 
negligible compared with that in the liquid phase. 
We may therefore with sufficient accuracy define 
the I’; as follows.'® If we compare a portion of the 
liquid containing unit area of free surface with 
another very nearly equal portion in the interior 
of the liquid then the former will contain I’; moles 
of each component 7 more than the latter. This 
definition is sufficient for the derivation given 
above of formula (4.3). It is however not a com- 
plete definition of the ['; because it leaves un- 
specified the relation between the quantity of 
liquid with free surface and the quantity of liquid 
in the interior. For instance we might choose the 
two quantities of liquid to be of identical mass, or 
to contain exactly the same total number of 
molecules, or to have exactly the same volume or 
to satisfy any one of an infinite number of alter- 


10 Cf. Guggenheim and Adam, Proc. Roy. Soc. A139, 
223 (1933). 

































694 E. A. 





native conditions. The choice is immaterial as 
far as the validity of (4.3) is concerned. This 
means that the actual values of the individual I; 
occurring in (4.3) are of little physical sig- 
nificance. It therefore seems desirable to trans- 
form formula (4.3) into another formula involving 
quantities with values independent of the quan- 
tity of liquid in the interior with which we 
compare the liquid having a free surface. We can 
do this and at the same time remove another 
defect in formula (4.3) namely, the following. For 
variations of composition at constant tempera- 
ture and pressure, the differentials d log p; are not 
all mutually independent, since they are related 
by the Gibbs-Duhem-Margules formula 


mRTd log pit+neRTd log po+---=0. (5.1) 


Let us now multiply (5.1) by T'1/m; and subtract 
from (4.3). We obtain 


Ne 
—dy= (=r. Rra log pe 


ny 


nN3 

+(1-“1r) ara log pst+---. (5.2) 
nN, 

We have thus eliminated p; and have obtained a 

formula involving only the independent variables 

po, ps: + :. We may write (5.2) in the abbreviated 

form 


—dy =T.RTd log pe 


+T3%RTd log ps+---, (5.3) 
where we define T,“, [3--- by 
To) =P2—(ne/m1)Ti, 
(5.4) 


T36)) = T3 -_ (3/m)T). 


We may thus regard F.“™, T'3, --- as measures 
of the quantities of the components 2, 3, --- ad- 
sorbed per unit area, when we take component 1 
as reference substance with I, zero by defini- 
tion. Now it is easy to prove that the quantities 
ro T3, --+ thus defined have values independ- 
ent of the quantity of homogeneous liquid with 
which we compare the liquid with unit free 
surface. For suppose we increase the quantity of 
homogeneous liquid in the ratio 1+x : 1. Then 
since the total number 7;+TI;; of moles of each 
component 7 in the system with unit free surface 
has a definite value, independent of any conven- 
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tion, each T; will have to be decreased by the 
amount n,x. If we denote by the symbol — the 
variation due to increasing the quantity of homo- 
geneous liquid used as reference to define the I, 
we thus have 


N;—Nn;+NXx, (5.5) 
TT —nx, (5.6) 
l,-T;— mx. (5.7) 
Hence we deduce 
T;—(n;/my)T iT; —(ni/ny)T 1 (5.8) 


or in other words the [;™ remain invariant. Thus 
formula (5.3) has two advantages over formula 
(4.3); the first that the differentials are inde- 
pendently variable and the second that the quan- 
tities [; have values completely defined by the 
system and independent of any arbitrary choice. 
Readers familiar with Gibbs’ treatment will 
recognize the quantities [T; to be the same as 
those introduced by Gibbs" and designated TI’ ,1) 
by him, but our definition is quite different from 
and much less abstract than that of Gibbs. 


6. RELATION BETWEEN SURFACE TENSION 
AND OSMOTIC PRESSURE 


We shall show that in the special case of a 
system of only two components there is a differ- 
ential relation between the surface tension and 
the osmotic pressure of the solution. This rela- 
tion is of small intrinsic interest and, in the 
author’s opinion, of little value. The reason for 
not ignoring it is that it is this relation, not the 
Gibbs adsorption formula, that was derived by 
Milner by an elementary proof. 

In a two component system if the component 1! 
is regarded as solvent and the component 2 is 
regarded as the solute, the osmotic pressure 7 
is related to the vapor pressure of the solvent by 
the formula 


dr= — (RT /v,)d log Pi, (6.1) 


where 7 is the molar volume of the solvent. This 
formula is thermodynamically accurate subject 
to the assumption that v; is independent of the 
composition (mixing without volume change) 
and independent of the pressure (compressibility 
negligible); the formula is independent of any 
assumption of ideality of the solution. But for a 


11 Gibbs, Collected Works, Vol. I, p. 234. 
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two component system the Gibbs-Duhem- 


Margules relation is 
nid log pi+Med log po=0 
and the Gibbs adsorption equation is 


—dy=T,RTd log p1+ T2RTd log po. 


(6.2) 


(6.3) 


If we eliminate ; and p2 between the three Eqs. 
(6.1), (6.2), (6.3) we obtain 


—dy= {T2—(n2 1) 1} (myv1/n2)dr. (6.4) 


If we introduce the quantity ['. defined above, 


we have 
—dy=T2™(ny1/n2)dr (6.5) 


or —dy=(T2/m)dr, (6.6) 


where m denotes the number of moles of solute 
per unit volume of solvent. Instead of using the 
quantity 2“, we might fix the values of T;, T's 
according to the arbitrary convention that the 
solution having a free surface be compared with a 
portion of homogeneous liquid having an equal 
volume (but a different mass, different total 
number of moles, and different numbers of moles 
of each component). We shall for convenience 
refer to this convention as the (v) convention. If 
we ignore any expansion or contraction when a 
molecule is taken from the interior of the liquid 
to the surface, the (v) convention is represented 
by the equation 


Tv, +T2v.=0, (6.7) 


the superscript (v) referring to the particular 
convention used to fix the values of I; and I>.!° 
If we eliminate [', between (6.7) and (6.4) we 
obtain 

—dy = Te (nv; + noWv2/n2)dr (6.8) 


—dy=(T2/c)drx, (6.9) 


where c denotes the number of moles of solute per 
unit volume of solution. Eq. (6.9) is the accurate 
expression of the formula deduced by Milner. 


7. MILNER’S FORMULA 


We shall now review briefly what Milner does 
in his alleged deduction of the Gibbs adsorption 
theorem. He tacitly assumes that there are only 
two components and uses the (v) convention to 
fix the values of I’, and I's. Subject to this arbi- 
trary convention he correctly derives the formula 
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—dy=(T2/c)dr. (7.1) 


He then assumes for z the formula 
x=RTc, (7.2) 
which is inaccurate even for an ideal solution. By 


combining (7.2) with (7.1) he obtains the in- 
accurate formula 


—dy=T2™RTd log c, (7.3) 


which he incorrectly calls the Gibbs adsorption 
theorem. 

The accurate formula for the osmotic pressure 
for an ideal solution of negligible compressi- 
bility is 


°=-— (RT/v:) log {21/(m1+M2)}. (7.4) 


If this is substituted into (7.1), we obtain 
NyV1+NoV2 


—dy=—-T.™ RTd log 
Ni+N» 


(7.5) 


N21 
or by a simple transformation 


NW +N Ne . 
—RTd log ———. (7.6) 


Ny+N2 


—dy = re™ 
Ny) 
This is the formula which Milner would have 
obtained if he had used the accurate formula for 
the osmotic pressure of an ideal solution. But it is 
to be noted that this formula is valid only for 
ideal solutions, whereas the Gibbs adsorption 
theorem is applicable and often applied to solu- 
tions that are by no means ideal. A further ob- 
jection to Milner’s formula (7.1) is that it holds 
only for T'2 defined according to the arbitrary 
(v) convention. The formula (6.5) for T2™ is 
different from Milner’s formula even for ideal 
solutions except at infinite dilution. Finally 
Milner’s formula is applicable only to binary 
solutions, whereas the Gibbs adsorption equation 
is, of course, applicable to mixtures of any 
number of components. 
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The Rate of Oxidation of Carbon Monoxide in Presence 
of Nitrogen Dioxide 


The effect of nitrogen dioxide has been studied by Crist 
and Roehling.! They suggested that the reaction had a 
complicated chain mechanism for low concentrations of 
nitrogen dioxide and that it was homogeneous for larger 
amounts of the catalyst. 

We have studied the reaction in Pyrex, quartz, alu- 
minum, and Pyrex treated with a potassium chloride 
solution. It was found that the reaction for low nitrogen 
dioxide pressures (0.06—5 mm) was greatly affected by 
the surface volume ratio and by the material of the reaction 
vessel. The rate was comparatively small in the Pyrex 
vessel washed with a potassium chloride solution. At 
higher pressures of nitrogen dioxide the reaction was unaf- 
fected by the nature of the vessel. 

The homogeneous portion appears to consist of the 
reactions: 

2 NO+0;=2 NOz, (1) 


2 NO.=2 NO+0:z, (2) 
NO.+CO=CO.+NO. (3) 


The experimental method was varied so as to isolate these 
reactions. In a mixture with the approximate concentra- 
tions, total nitrogen dioxide=10 mm, oxygen =600 mm, 
and carbon monoxide =50 mm, the reaction was found to 
be first order with respect to both carbon monoxide and 
nitrogen dioxide and proportional to the square root of the 
oxygen pressure. When the oxygen concentration is large 
it would be practically constant and the rate would be 
expressed by 


d 
-— =k:[NO2][CO]. 


Then, in order to obtain k; we calculate the nitrogen 
dioxide on the assumption that it is constant at the 
equilibrium concentration. The values of kz; at 385°C and 
527°C are, respectively, 1.12 10?(+0.07) sec. and 8.50 
X 108(+0.20) sec.-1. 

In order that reaction (1) should determine the rate, 
the mixture taken was approximately as follows: nitrogen 
dioxide=10 mm, carbon monoxide=600mm and oxy- 
gen = 20 mm. Since the oxygen pressure is low the nitrogen 
dioxide concentration will also be small and the amount of 
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nitrogen dioxide lost by reaction (2) will be relatively small. 
The results show that the reaction under these conditions 
is first order with respect to oxygen, independent of the 
carbon monoxide, and proportional to the square of the 
nitric oxide when its concentration is not too large. This is 
in accordance with the equation 


dO 
— =k, [NO}(0.]. 


The results are in Table I. 


TABLE I. Rate constants for the oxidation of nitric oxide. 








Temperature =385°C Temperature =527°C 





ki X10 
cc? mole? sec.~! 


Pressure in mm 
of NO (Total) 


ki X 1079 
cc? mole sec. 


Pressure in mm 
of NO (Total) 





7.37 

7.65 
10.8 
13.4 
15.7 
15.8 
16.0 
19.0 
28.0 


6.80 
10.0 
115 
15.4 
17.4 
25.2 


2.02 
2.05 
1.80 
1.48 
1.47 
0.79 


4.21 
4.37 
4.20 
3.92 
3.38 
3.37 
3.45 
3.68 
2.47 

















The decrease of the rate constant at higher nitrogen 
dioxide pressures could be due to increased nitrogen dioxide 
decomposition. The best value for k;: would then be that 
obtained at the lower concentration of nitrogen dioxide. 
The value of the constant at 385°C is in fairly good agree- 
ment with the work of Bodenstein.2 The constant at 
527°C is in accordance with the predictions that the 
temperature coefficient should become positive at tem- 
peratures above 400°C. 

R. H. Crist 
G. M. CaLHoun 


Department of Chemistry, 
Columbia University, 
New York, N. Y., 

August 31, 1936. 


1 Crist and Roehling, J. Am. Chem. Soc. 57, 2196 (1935). ; 
2 Kassel, Kinetics of Homogeneous Gas Reactions (Chemical Catalog 
Co., 1932), p. 157. 
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